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UNSTEADY  AND  SEPARATED  FLOWS 
AFOSR  Contract  F49620-86-C-0134 


Abstract 

This  is  the  final  report  of  research  on  Unsteady  and  Separated  Flows  (AFOSR 
Contract  F49620-86-C-0134)  conducted  at  the  California  Institute  of  Technology. 
The  overall  objectives  of  this  research  were  to 

i.  extend  understanding  of  the  physics  and  mechanisms  of  unsteady  and  sep¬ 
arated  flows,  particularly  of  the  role  of  large  vortical  structures,  for  various 
boundary  conditions  (steady  and  time  varying),  and  for  three-dimensional, 
as  well  as  two-dimensional  geometries; 

ii.  develop  techniques  to  control  such  flows  for  desired  effects;  and 

iii.  seek  ways  of  improving  methods  of  computation  for  unsteady  separated 
flows. 
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1.  Introduction 

Our  work  under  the  sponsorship  of  this  contract  on  unsteady  and  separated 
flows  built  on  the  study  of  large  scale  vortical  structures  that  began  at  GALCIT. 
When  studying  these  flows,  we  note  that  they  are  dominated  by  large  vortical 
structures.  These  structures  result  from  the  basic  underlying  instabilities  of  shear  in 
flows.  As  a  consequence,  the  interference  or  removal  of  instabilities  might  provide  an 
effective  means  of  controlling  the  flow  field.  On  the  other  hand,  vortical  structures 
are  useful  in  flows  that  rely  on  mixing  and  combustion.  Control  of  large  scale 
structures  in  these  flows  can  yield  an  unsteady  flow  that  is  mere  desirable  than  the 
naturally  occurring  flow. 

In  the  present  investigations,  simple  geometries  were  studied.  Experiments 
and  simulations  examined  the  addition  of  unsteady  motions,  of  flat  plates  inclined 
normal  to  the  flow,  and  circular  cylinders,  to  the  conventional  steady  motions  of 
these  bodies.  The  simulations  also  investigated  the  (three-dimensional)  interactions 
of  vortex  rings  and  the  startiq)  flow  of  a  square  flat  plate  oriented  normal  to  the 
direction  of  motion  (bluff'  body  configuration). 

Laboratory  experiments  and  measurements  on  flow  past  flat  plates  nominally 
normal  to  the  direction  of  motion  were  performed.  The  objectives  of  this  research 
were  to  elucidate  the  dynamics  of  large  vortical  structure  in  globally  steady  and 
unsteady  separated  flows,  and  to  develop  techniques  to  modify  tlieir  beliavior  and 
control  them  for  desired  effect.  In-line  accelerations  and  transverse  oscillations 
of  flat  plates  were  evaluated  to  modify  the  flow  structure.  These  investigations 
examined  the  behavior  of  the  large  vortical  structures  in  the  resulting  unsteady  and 
separated  flow,  and  the  associated  forces  on  the  flat  plate. 

Vortex  methods  were  used  to  study  two-  and  three-dimensional  flows.  The 
objectives  of  this  research  were  to  develop  methods  to  examine  the  large-scale  sep¬ 
arated  vortices  in  the  wake  of  bluff  bodies,  study  the  influence  of  transient  motions 
of  the  body  on  these  dynamics,  and  investigate  control  strategies  to  produce  de¬ 
sired  effects.  Two-dimensional  simulations  of  the  flow  past  a  flat  plate  executing 
unsteady  motions  were  coordinated  with  laboratory  experiments  for  the  purpose 
of  elucidating  the  physics  of  natural  and  forced  vortex  shedding.  Lateral,  in-line, 
and  rotary  motions  of  the  plate  were  studied.  Three-dimensional  simulations  of  the 
interaction  of  vortex  rings  and  the  flow  past  an  impulsively  started  square  plate 
were  also  studied. 

Investigations  on  the  flow  j)ast  a  cylinder  executing  rotary  motions  were  made. 
The  objective  of  this  part  of  the  research  was  to  implement  a  control  scheme  to 
take  a  basically  unstable  flow  situation  and  either  render  it  effectively  stable,  at 
least  for  some  time  interval,  or  to  exploit  the  nature  of  the  instability  for  a  par¬ 
ticular  desired  purpose.  We  initially  intended  to  use  adaptive  feedback  control  to 
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remove  instabilities  in  these  flows.  Informal  investigations  had  suggested  that  an 
artificial  feedback  path  introduced  into  the  flow  could  provide  a  measure  of  control 
with  rather  modest  control  perturbations.  In  the  course  of  the  present  investiga¬ 
tions,  however,  we  discovered  that  the  natural  feedback  paths  in  wake  flows,  when 
combined  with  open  loop  forcing  (albeit  requiring  larger  control  inputs),  provide 
considerable  control  authority  over  the  flow  structure. 

The  technical  discussion  for  this  report  is  roughly  divided  into  three  parts  which 
describe,  experiments  with  flat  plates  (fixed  separation  points)  in  Sec.  2;  experi¬ 
ments  with  circular  cylinders  (variable  separation  points)  in  Sec.  4;  and  additional 
two-  and  three-dimensional  numerical  simulations  in  Sec.  3. 
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2.  Investigations  of  flow  past  a  flat  plate  (experimental  and  numerical) 

The  Laboratory  experiments  and  measurements  on  flow  past  flat  plates  nomi¬ 
nally  at  90°  angle  of  attack  were  performed  in  the  X-Y  Tow  Tank  in  close  interaction 
with  corresponding  numerical  simulations.  The  findings  may  be  briefly  summarized 
as  follows: 

a.  In  accelerated  starts,  the  flat-plate  motion,  the  forces  and  flow  field  mea¬ 
surements  in  the  Tow  Tank  during  the  start-up  (acceleration)  phase  of  the 
plate  motion  agree  well  with  those  computed  in  an  exactly  two  dimensional 
simulation  using  vortex  methods.  This  indicates  that  the  laboratory  flow 
is  nearly  two  dimensional  during  the  initial  acceleration  phase  and  that 
the  numerical  simulation  represents  the  flow  well. 

b.  After  the  acceleration  period,  the  flow  fields  and  forces  in  the  still  develop¬ 
ing  (transient)  flows  develop  differently  in  the  laboratory  and  in  simulation. 

c.  After  about  fifty  body  chords  of  travel  at  constant  velocity  the  flow  reaches 
the  steady,  vortex  shedding  state.  The  value  of  the  mean  drag  coefficient 
in  our  experiments  {Cd  =  2),  agrees  well  with  previously  known  results. 
The  corresponding  value  from  the  simulation  differ  considerably  {Cd  =  3.5 
and  Strouhal  number,  S  =  0.13,  respectively).  These  differences  may  be 
attributable  to  the  effects  of  three  dimensionality  in  the  laboratory  flows. 
These  include  effects  of  spanwise  non-uniformities  in  the  shedding  which 
develop  from  instabilities  and  imperfect  boundary  (end)  conditions. 

d.  Experiments  and  numerical  simulations  of  the  flat  plate  in  forced,  periodic 
plunging  motion*  at  the  natural  Strouhal  frequency  have  shown  that,  when 
steady  state  was  reached,  the  difference  between  the  estimated  Cd  in  the 
laboratory  experiments  {Cd  =  2.2),  and  simulations  {Cd  =  2.9),  was 
narrowed.  The  decrease  seen  in  the  simulations  must  be  connected  with 
changing  phase  relations  in  the  vortex  interactions  with  the  plate;  in  the 
experiment  similar  changes  should  be  expected  but  the  overriding  effect  is 
probably  from  a  reduction  of  three-dimensionality. 

e.  These  investigations  suggest  that,  in  accelerating  motion,  it  is  the  vorticity 
(of  both  signs)  which  is  generated  is  well  correlated  along  the  entire  span 
of  the  flat  plate.  This  may  reduce  the  relative  importance  of  small  three 
dimensional  effects  when  compared  with  the  larger  two  dimensional  motion 
introduced  by  the  plate  acceleration. 


Note  that  this  motion  involves  unsteady  lateral  acceleration. 


3.  Numerical  simulation  of  separated  flows 

In  this  work,  numerous  techniques  have  been  developed  which  bring  the  vortex 
method  paradigm  closer  to  modeling  real  separated  flows  in  both  two  and  three 
dimensions.  See  Chua  (1990). 

In  two  dimensions,  flow  past  a  flat  plate  nominally  at  90°  angle  of  attack  was 
investigated.  The  objective  of  this  study  was  to  determine  the  effect  of  forced 
oscillations  in  the  motion  of  the  plate  on  its  time-dependent  forces.  Both  lateral 
(see  Sec.  2)  and  rotational  motions  of  the  plate  were  also  investigated. 

In  three  dimensions,  two  problems  involving  strong  interactions  of  vortex  tubes 
were  carried  out.  In  both  cases,  the  models  we  developed  successfully  characterized 
the  complex  motions  and  structure  of  the  vortex  tubes,  and  the  reconnection  of 
vortex  lines.  The  two  prototypical  three  dimensional  flows  which  were  investigated 
are: 

a.  The  interaction  of  vortex  rings  was  studied.  These  results  compared  favor¬ 
ably  with  the  experimental  data  of  Schatzle  (1978),  and  other  numerical 
results. 

b.  The  early  stages  of  flow  past  an  impulsively  started  square  plate  were  also 
studied. 
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4.  Active  control  of  separated  flows 

The  type  of  active  control  of  the  flow  described  in  this  facet  of  the  investi¬ 
gations  relies  on  the  forces  between  the  cylinder  and  the  surrounding  flow.  It  is 
this  interaction  that  regulates  the  ejection  of  circulation  into  the  separated  flow 
region  and  actuates  the  entire  flow  field.  This  is  exemplified  by  the  control  of  the 
structure  demonstrated  by  the  pluise  locked  flows  in  Appendix  A  and  Chs.  3  -  5  of 
Appendix  B  as  well  as  the  up-stream  influence  characterized  by  the  lifting  cylinder 
flows  described  in  Appendix  B,  Ch.  5.  In  addition,  this  type  of  control  exploits 
the  natural  tendency  to  have  vortical  structures  in  the  flow.  This  is  in  contrast  to 
the  more  traditional  control  philosophy  that  unsteadiness  is  something  that  should 
be  canceled,  e.g.,  Liepmann,  et  al.  (1982),  Liepmann  &  Nosenchuck  (1982),  and 
Ffowcs-Williams  &  Zhao  (1989). 

The  efficacy  of  oscillatory  cylinder  rotation  as  an  actuation  mechanism  for 
actively  controlling  the  cylinder  wake  is  examined  in  Appendix  A  and  Appendix  B, 
Ch.  3.  It  is  shown  that  considerable  control  can  be  exercised  over  the  structure  in  the 
wake  with  such  forcing.  In  particular,  a  large  increase,  or  decrease,  in  the  resulting 
displacement  thickness,  estimated  cylinder  drag,  and  associated  mixing  with  the  free 
stream  can  be  achieved,  depending  on  the  frequency  and  amplitude  of  oscillation. 
Not  surprisingly,  these  residts  show  that  working  in  a  control  domain  in  which 
the  structures  shed  are  synchronous  with  the  forcing  provides  the  greatest  control 
authority  over  the  wake  structure.  While  these  results  were  obtained  for  moderate 
Reynolds  numbers,  i.c..  Re  ~  10'*.  there  are  indications  that  the  description  of 
the  flow  phenomena  presented  here  is  qualitatively  the  same  over  a  large  range  of 
Reynolds  numbers.  The  streamline  data  in  Appendix  B,  Ch.  4,  show  that  under 
control  (forcing),  the  vortical  stinicturcs  are  formed  very  close  to  the  cylinder.  This 
is  an  indication  that  the  formation  of  the  wake  structure,  under  control,  has  more  to 
do  with  the  controlled  ejection  of  circulation,  and  the  associated  forces  on  the  body, 
than  with  the  stability  characteristics  of  some  generic  unforced  wake.  This  forced 
shedding  is  in  contrast  to  the  natural  (unforced)  shedding  of  vortices,  where  the 
structures  form  only  after  the  vorticity  has  been  shed  into  the  wake,  from  more-or- 
less  steady  separation  points.  Hence,  it  is  possible  that  the  stability  characteristics 
of  the  wake  can  play  a  more  significant  role  in  the  unforced  case  than  in  the  forced. 
In  addition,  it  is  shown  that  under  control,  very  little  of  the  circulation  introduced 
at  the  cylinder  wall  appears  in  the  wake.  Evidently,  most  of  the  positive  and 
negative  vorticity  introduced  into  the  boundary  layer  fluid  combines  before  the 
fluid  is  allowed  to  enter  the  wake. 

A  method  for  estimating  the  mean  lift  of  a  rotating  cylinder  is  presented  in 
Appendix  B,  Ch.  5.  Estimates  based  on  this  method  compare  favorably  with  similar 
data  published  for  steadily  rotating  cylinders.  It  is  also  demonstrated  that  the 
addition  of  forced  rotary  oscillations  to  the  steady  rotation  of  the  cylinder  increases 
the  lift  coefficient  in  the  cases  when;  the  wake  would  normally  be  separated  in  the 
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steadily  rotating  case,  but  decreases  it  otherwise. 

This  research  is  documented  in  more  detail  in  Appendix  A  “Rotary  oscillations 
control  of  a  cylinder  wake”  (Tokumaru  &  Dimotakis  1991),  Appendix  B  “Active 
Control  of  the  Flow  Past  a  Cylinder  Executing  Rotary  Motions”  (Tokumaru  1991). 
Pre-print  drafts  of  the  papers  “The  mean  lift  of  a  rotating  and  oscillating  cylinder” 
(Appendix  C),  and  “Streamline  patterns  in  the  neighborhood  of  a  cylinder  executing 
rotary  oscillations  in  a  uniform  flow”  (Appendix  D),  by  Tokumaru  &  Dimotakis  are 
also  included. 
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5.  Personnel 

In  addition  to  the  principal  Investigators, 

P.  Dimotakis:  Professor  of  Aeronautics  and  Applied  Physics; 

A.  Leonard;  Professor  of  Aeronautics; 

A.  Roshko:  Theodore  von  Kanuan  Professor  of  Aeronautics; 

other  personnel  who  participated  directly  in  the  effort  during  the  course  of  these 
investigations  are  listed  below: 

K.  Chua:  Graduate  Research  Assista’  t  in  Aeronautics;** 

D.  Lisoski:  Graduate  Research  Assistant  in  Aeronautics; 

P.  Tokuinaru:  Graduate  Research  Assistant  in  Aeronautics;^ 

F.  Noca:  Graduate  Research  Assistant  in  Aeronautics; 


**  Presently  at  Continuum  Dynamics,  Inc.,  Princeton. 
1  Presently,  Research  Fellow  in  Aeronautics,  Caltech. 
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Rotary  oscillation  control  of  a  cylinder  wake 

By  P.  T.  TOKUMARU  and  P.  E.  DIMOTAKIS 
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Exploratory  experiments  have  been  performed  on  eireular  eylinders  executing  fo!  eed 
rotary  oscillations  in  a  steady  uniform  flow.  Flow  visualization  and  wake  profile 
measurements  at  moderate  Reynolds  numbers  have  shown  that  a  considerable 
amount  of  control  can  be  exerted  over  the  structure  of  the  wake  by  such  means.  In 
particular,  a  large  increase,  or  decrease,  in  the  resulting  displacement  thickness, 
estimated  cylinder  drag,  and  associated  mixing  with  the  free  stream  can  be  achieved, 
depending  on  the  frequency  and  amplitude  of  oscillation. 


1.  Introduction 

It  has  long  been  appreciated  that  the  dynamics  of  the  flow  over  a  bluff  body  are 
dominated  by  large-scale  vortical  structures  shed  in  its  wake.  Weihs  (1972)  proposed 
a  methodology  for  predicting  the  structure  of  such  a  wake  and  found  that  it  could 
become  wider,  narrower,  or  remain  the  same,  depending  on  the  initial  strength  and 
spacing  of  these  vortices.  His  work  was  primarily  aimed  at  describing  the  wake 
behind  an  oscillating  airfoil  (e.g.  Bratt  1953;  Koochesfahani  1987).  Nevertheless, 
Weihs's  results  do  suggest  the  potential  for  effective  control  of  wakes  in  general. 

The  effects  of  periodic  forcing  on  a  wake  were  inv^'estigated  by  Roberts  (1985),  and 
Roberts  &  Roshko  (1985).  They  demonstrated  that  a  significant  amount  of  control 
over  the  wake  vortex  structure  and  the  associated  mixing  was  possible,  by  pulsing 
one  stream  relative  to  the  other.  Okajima,  Takata  &  x\sanuma  (1975)  examined  the 
forces  acting  on  a  rotationally  oscillating  cylinder  for  Reynolds  numbers,  based  on 
cylinder  diameter,  in  the  range 

40  < <  fi  X  1(F,  (1) 


where  is  the  free-stream  v'elocity  far  ahead  of  the  cylinder,  d  is  the  cylinder 
diameter,  and  c  is  the  kinematic  viscosity.  His  measurements  were  for  a  normalized 
{)eak  rotation  rate. 


d 


e,. 


(-’) 


of  0.2  to  1.0.  and  a  forcing  Strouhal  number, 

s,  =  if /.  (3) 

ranging  from  0.05  to  0.3,  where  is  the  peak  rotational  rate  of  the  cylinder,  is  the 
peak  circumferential  velocity,  and  /  is  the  forcing  frequency.  He  noted  a 
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synchronization  '  similar  to  that  observed,  for  example,  by  Bishop  &  Hassan  (liHU). 
Koopman  (1007).  and  more  recently.  Ongoren  &  Rockwell  (10SSa./>).  and  by 
Williamson  &  Roshko  (1988)  for  a  cylinder  in  transverse  and  in-line  oscillation.  Both 
Okajima  et  <tl.  and  Bishop  &  Ha.ssan  reported  a  hysteresis  in  their  measurements  for 
increasing  ajid  decreasing  Inve.stigations  at  comparable  forcing  frc(|uencies. 
amplitudes  and  lie  were  also  performed  by  Wu,  Mo  &  N’akili  (1989).  Taneda  (1978) 
demonstrated  that  in  the  range  30  ^ /^c  ^  300,  the  ‘dead  water'  region  behind  a 
cylinder  can  be  removed  for  sufficiently  large  (Q^  >  7  to  27.  depending  on  such 
factors  as  the  cylinder  span  and  tank  width).  For  somewhat  larger  Q^,  the 
streamlines  converged  downstream.  Similar  behaviour  was  also  documented  by 
Williams  &  Amato  (1988)  using  a  line  of  unsteady  pulsing  jets  embedded  in  the 
trailing  edge  of  the  cylinder,  at  a  comparable  Reynolds  number  of  370. 

In  the  work  to  be  discussed  here,  we  examined  the  efficacy  of  forced  rotary 
oscillations  of  a  (fixed  axis)  cylinder  at  a  moderate  Reynolds  number  of  1.5  x  10“*,  for 
the  purpose  of  controlling  the  unsteady  separated  flow  in  its  wake.  The  effect  of  the 
cylinder  oscillations  on  the  structure  of  the  wake  was  investigated  using  flow 
visualization,  as  well  as  surveys  of  the  resulting  wakes  using  laser- Doppler  velocity 
measurements. 

2.  Experimental  facility 

The  experiments  documented  here  were  performed  in  the  18  in.  wide,  by  19  in. 
deep,  low  speed  water  channel,  as  well  as  the  20  x  20  in.  free  surface  water  channel 
at  GALCIT.  A  4  in.  diameter  Plexiglas  cylinder  was  supported  10  in.  above  the 
bottom  of  the  channel,  using  0.5  in.  thick  Plexiglas  fairings  placed  flush  to  the 
sidewalls  of  the  channel.  Drive  belts  and  pulleys  were  entirely  enclosed  in  the  fairings 
and  did  not  interfere  with  the  flow.  Rotation  about  the  axis  of  the  cylinder  was 
achieved  using  a  high-performance  servo-controlled  motor,  capable  of  tracking  an 
arbitrary  command  signal  within  the  bandwidth  and  slewing  rates  that  were 
investigated.  This  command  signal  was  generated  by  a  function  generator,  or  a 
computer  digital-to-analog  converter  output  channel.  Flow  visualization  was 
accomplished  by  introducing  dye  upstream  of  the  cylinder,  or  through  small  holes  on 
the  cylinder  surface  at  midspan.  Streamwise  mean  and  r.m.s.  velocity  profiles  were 
measured  using  a  Bragg  cell  frequency-shifted  laser- Doppler  velocimeter. 

3.  Results  and  discussion 

The  cylinder  was  rotated  sinusoidally  in  time  at  a  forcing  frequency  /  and  a 
normalized  peak  rotational  rate  of  12,  (sec  (2)),  i.e. 

=  Q(t)  =  12,  sin  (2Kft).  (4) 

12,  was  chosen  such  that  the  peak  circumferential  velocity  of  the  cylinder  would  be 
comparable  with  the  velocity  just  outside  the  boundary  layer  of  the  cylinder 
(approximately  two  times  the  free-stream  velocity).  Wv  anticipated  that  amplitudes 
of  this  magnitude  or  greater  would  be  necessary  to  send  the  vorticity  stored  in  the 
boundary  la\'er  into  the  wake  in  a  regulated  manner.  For  these  experiments, 
we  examined  the  control  parameters  12,  and  Sf  in  the  range  0  ^  12,  ^  16  and 
0.17  ^  3.3.  The  free-stream  velocity  was  approximately  15  cm/s.  yielding  a 
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Reynolds  number  based  on  cylinder  diameter  o{  Re  =  1.5  x  10^.  In  addition,  wake 
mean  and  r.rn.s.  velocity  profiles  were  measured  at  a  single  streamwise  station, 
located  4.5  diameters  downstream  of  the  cylinder  axis.  See  figure  1. 


3.1.  Displacement  thickness  and  drag  coefficient  estimation 


In  order  to  provide  a  quantitative  measure  of  the  effects  of  forcing  on  the  flow  over 
a  cylinder  and  the  resulting  wake,  the  cylinder  drag  and  wake  displacement  thickness 
were  estimated  from  wake  streamwise  mean  and  r.rn.s.  velocity  profiles.  For  flow 
that  is  two-dimensional,  we  can  define  the  displacement  thickness  in  terms  of  the 
integral 


S* 


■  II 


0  J 


di/. 


(5) 


In  this  expression,  y  ~  y/h  is  the  normalized  vertical  (cross- stream)  position,  u{y)  is 
the  mean  streamwise  ’-olocity,  y  and  h  are  the  vertical  position  and  water  channel 
depth  respectively,  a  ’  fg  is  the  velocity  in  the  (free-stream)  region  outside  the 
cylinder  wake  and  water  channel  boundary  layers.  See  figure  1.  By  analogy,  we  may 
define  the  'displaced  area’  for  a  three-dimensional  flow.  i.e. 


A* 

T 


Uo  J 


dyd^. 


(6) 


where  A  is  the  test-section  cross-sectional  area,  and  ^  =  zjh  is  the  normalized 
spanwisc  coordinate,  with  z  the  spanwise  coordinate,  and  b  the  water  channel  span. 
Continuity  then  yields  for  the  free-stream  velocity  far  ahead  of  the  cylinder, 

^^00  =  1  (7) 


Using  (6)  and  (7)  we  then  find 


A  (/„• 


(8) 


This  expression  is  useful  in  that  it  provides  information  about  the  flow  over  an  entire 
cross-section  of  the  water  channel,  while  requiring  onlv  the  measurement  of  and 

In  order  to  remove  the  contribution  of  the  initial  boundary  layer  (in  the  absence 
of  the  cylinder)  to  leading  behaviour,  it  is  useful  to  compare  (with  the  cylinder  in 
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Figi  rk  2.  Variation  of  A*/A  with  Sf  and 


place)  with  the  undisturbed  flow  velocity  in  the  absence  of  the  cylinder  Uq,  at  the 
same  strearnwise  location.  This  yields  an  expression  for  the  approximate  displaced 
area  of  the  wake : 


A  ~  f/„- 


(9) 


Figure  2  depicts  the  variation  of  A* /A  with  and  ^2■^.  The  dashed  curve  in  this  figure 
denotes  the  A* /A  level  for  the  unforced  case.  The  parameter  space  above  this  curve 
represents  an  increa.se  in  A* /A  over  the  unforced  case,  and  below  represetits  a 
decrease.  Note  that  the  displacement  area  can  be  made  substantially  larger  or 
smaller  by  varying  the  frequency  and  amplitude  of  oscillation.  The  minimum  for 
A* /A  over  the  range  shown  occurs  in  the  neighbourhood  of  1.  and  %  3.  and 
is  roughly  half  of  the  unforced  case.  Measurements  for  A* /A  were  taken  in  the  range 
0.17  <  ^  3.3  in  steps  of  about  0.17.  and  0  ^  ^  16  in  steps  of  I .  Each  of  the  data 

points  presented  in  figures  2  and  3  were  averaged  over  2  minutes  (approximately  40 
natural  unforced  shedding  cycles).  Wake  midspan  velocity  profiles  were  measured  for 
Sj  in  the  range  of  O.lo  to  1.4.  holding  fixed  at  2.  Figure  3  compares  the  wake 
profiles  for  several  S^.  and  =  2. 

For  finite  test  sections  and  flow  that  is  two-dimensional  in  the  mean,  one  can  show 
that  the  sectional  drag  coefficient  can  be  estimated  by  the  expression  (e.g.  Dimotakis 
1978): 


2h/d 

{l-S*/h)^ 


(10) 


In  this  expression,  u'  and  v'  are  the  strearnwise  and  cross-stream  r.m.s..  velocities  and 
8*  is  the  displacement  thickne.ss  defined  in  (5).  The  displacement  thickness  and  drag 
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cocfficiont  ostimatos  suinmarizod  in  fi^un*  4  wero  calculated  from  mean  mids{)an 
velocity  f)rofilos.  using  (o)  and  (10).  negle<  ting  the  fluctuating  terms.  As  can  be  .seen, 
there  is  a  broad  minimum  in  and  S*/f(  aroutul  =  1.  The  estimated  for  the 
unforced  case  is  a  factor  of  si.x  greater  than  for  this  forced  ease.  Similarly,  there  i.s  a 
factor  of  five  reduction  in  S*/h.  Xote  that  the  value  of  S*/k  for  a  particular  Sj-  and 
is  generally  less  than  the  corresponding  value  of  A* /A,  indicating  that  there  is  an 
increase  in  8*/h  outside  the  midspan  region.  Xote  also  that  since  C’d  increases  with 
8*/h  it  must  also  increase  outside  this  midspan  region.  This  spanwise  variation  in 
8*/h  is  due  to  the  interaction  of  the  .sidewalls  with  the  .shed  vortices.  See  Roberts 
(1985),  Roberts  &  Roshko  (1985).  Kooche.sfahani  (1987).  and  Kurosaka  cf  al.  (1988). 
Another  local  minimum  e.xists  near  ^  9.2.  The  How  visualization  in  Hgure  5 
suggests  that  this  corresponds  to  a  mode  in  which  the  cylinder  releases  tiro  vortict's 
of  opposite  sign  during  each  half-cycle.  When  one  of  the.se  two  vortices  .shed  from  the 
bottom  ap{)ears  on  tof).  a  di[)  is  ob.scrved  in  the  corresponding  velocity  proHle.  Sei' 
the  circles  on  figure  3. 

We  should  mention  that  Okajima  cf  al.  (1975)  observed  a  hysteresis  with  itn-reasing 
and  decreasing  Sj  and  =  0.2.  We  al.so  observed  some  hysteresis  for  small  values  of 
but  none  was  discernible  in  e.xperiments  with  ^  2. 

Finally,  we  suspect  that  had  it  not  been  for  the  finite  height  of  the  test  section,  the 
wake  could  have  been  made  even  larger  by  forcing  at  Sj  <  0.2. 

3.2.  Flow  visualization 

A  clear  picture  of  the  control  that  can  be  exerci.sed  over  the  .separated  flow  in  the 
wake  of  the  cylinder  emerges  from  our  flow  visualization  data.  These  demonstrate 
the  overall  narrowing  or  widening  of  the  wake,  depending  on  the  values  of  and  f?,. 
the  control  parameters,  as  well  as  the  formation  of  distinct  vortical  structures  in  the 
wake. 

In  figure  6  (plate  1).  dye  is  introduced  ahead  of  the  cylinder  by  a  comb  of  dye 
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Figure  6.  Close  up  and  full  flow  visualization  views  with  and  without  active  control, 
(a)  Unforced  (non-rotating)  cylinder  in  uniform  flow;  (fc)  5/  ~  1,  ~  3. 
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Fi<!1  KK  7.  Forced  sheddinfj  at  fixed  =  1  (lie  =  1.5  x  HH).  (n)  Q^  =  H.  (h)  i?,  =  l(i.  (r)  =  :{2. 

(d)  =  (»4. 


injt'ctors.  In  the  unforced  ease,  figure  ()(«).  the  dye  can  f)e  seen  to  be  dispersed  and 
mixed  across  the  full  height  of  the  test  section,  within  a  few  diameters  downstream 
of  the  cylinder.  In  contrast,  the  dye  marker  in  the  forced  ca.se  corresponding  to  the 
minimum  wake  width  found  in  §3.1.  figure  (}(/*)■  occupies  ap|)ro.\imately  the  same 
fraction  of  the  test  .section  height  oji  leaving  as  it  did  upon  entering. 

The  data  in  figures  7  and  S  (recorded  witli  dye  issuing  from  the  surface  of  the 
evlinder)  illustrate  how  the  wake  mav  be  made  witler  or  narrower  bv  varvitio  the 
control  parameters  and  12,.  Wcihs  (1972)  examined  this  phenomenon  by 

considering  the  mutual  influence  of  the  vortices  in  the  wake  and  those  being  shed. 
Although  his  analysis  is  highly  idealized,  considering  only  point  vortices  in  potential 
flow,  it  does  provide  an  intuitive  physical  argument  for  the  various  trajectories  taken 
f)y  the  vortices.  Additional  flow  visualization  data  (al.so  recorded  with  dye  issuing 
from  the  surface  of  the  evlinder).  at  fixed  =  8  and  increasing  Sj.  are  def)ieted  in 
figuri's  9  and  19. 
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Fujurk  8.  .Mode  II  .  wake  structure  is  synchronized  with  the  forcing  frequency 
(i2,  =  8.  /?e  =  1.5  X  li)*).  (a)  =  0.3. '(6)  =  0.5.  (f)  =  0.7,  (d)  N,  =  0.9.  ' 


Several  qualitatively  different  vortex  shedding  modes  were  observed.  They  are 
presented  here  as  modes  I  to  IV.  in  order  of  increasing  frequency.  In  the  forced  mode 
I  (e.g.  figure  5),  the  cylinder  releases  two  vortices  of  the  same  sign  per  half-cycle.  This 
second  vortex  may  be  attributed  to  an  additional  separation  caused  by  the  most 
recently  generated  vortex.  In  mode  II  (e.g.  figure  8).  the  wake  .structure  is 
synchronized  with  the  forced  cylinder  oscillation  and  fx?rsists  beyond  the  end  of  the 
test  section.  In  mode  III  (e.g.  figure  9).  the  near-wake  structure  is  also  synchronized, 
but  becomes  unstable  and  evolves  into  a  structure  with  a  lower  spatial  frequency 
some  distance  downstream  of  the  cylinder.  In  mode  IV'  (e.g.  figure  10).  the  effect  of 
forcing  is  primarily  observed  in  the  shear  layers  separating  from  the  cylinder.  Note 
that  while  the  flow  visualization  photographs  of  the  forced  cylinder  wake,  at  the 
largest  values  of  (cf.  figure  10.  bottom),  resemble  those  for  the  unforced  ca.se  (cf. 
figure  11).  we  find  that  A*/A  in  these  forced  cases  is  still  noticeably  less  than  in  the 
unforced  case. 
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Tlu'  transition  between  modes  II  and  IV  occurs  gradually  through  mode  III. 
Figures  8-10  demonstrate  how  the  ‘breakup’  of  the  synchronized  structure  in  this 
mode  moves  toward  the  cylinder  with  increasing  until  finally  only  the  shear  layers 
are  forced.  This  breakup  in  the  flow  structure  in  modes  III  and  IV  may  be 
attributable  to  the  stability  characteristics  of  the  evolving  mean  velocity  profiles  (as 
suggested  by  Cimbala,  Nagib  &  Roshko  1988  for  lower  Reynolds  number  flows  far 
downstream),  vortex  coalescence,  three-dimensional  effects,  or  some  combination 
thereof. 

Behaviour  similar  to  modes  I,  II,  and  III  was  previously  observed  by  Roberts 
(198.5).  and  Hol)erts  &  Roshko  (1985)  for  the  case  of  forcing  the  wake  behind  a 
splitter  f)late.  .More  generally,  the  bluff-body  wake  can  be  compared  to  a  nonlinear 
oscillator  (e.g.  Provansal,  Mathis  &  Boyer  1987)  which,  when  forced,  can  reveal 
synchronized  or  phase -locked  regions  in  their  parameter  spaces,  as  well  as  regions 
which  may  be  periodic,  quasi-periodic,  chaotic,  unstable,  or  even  stable.  See.  foi- 
example,  (iuckenheimer  &  Holmes  (1983,  Ch.  2).  Note  that  both  Taneda  (1978).  and 
Williams  &  Amato  (1988)  show  stable  (non-vortex-shedding)  regions  in  their 
frcMiiicncy  amplitude  parameter  spaces. 

Figure  12  shows  the  variation  of  the  normalized  phase  speed. 


with  A  is  the  initial  streamwise  vortex  spacing,  measured  from  photographs  like 
those  in  figures  8-10.  Note  the  rapid  change  in  5^  around  Sf=i.  This  change  is 
reminiscent  of  Tritton’s  (1959)  high-speed  and  low'-speed  vortex  shedding  modes. 
Our  mode  II  resembles  Tritton’s  high-speed  mode  and  our  mode  III.  his  low-speed 
mode. 

Observations  of  unforced  vortex  shedding  show  that  vorticity  is  introduced  into 
the  How  from  a  more-or-less  steady  separation  point,  with  the  resulting  vortical 
structures  forming  in  the  wake  of  the  cvlind(“r.  This  is  in  contrast  to  the  forced 
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slu'tidiiiji  process.  r('prcsciitativc  of  nuxlcs  1  111.  which  is  charactcrizt'd  by  an 
unsteady  separation  point  .  The  forced  sheddinrt  j)roeess  for  I  and  O,  ^  2  is 
illustrated  in  figure  Hi.  The  shaded  areas  in  this  figure  denote  fluid  inarketl  with  dye 
introduet'd  just  ahead  of  th(‘  ten  o'trloek  jiosition  on  the  cylinder.  This  forcing 
replaces  the  natural  shedding  ineehanisin  with  one  that  gathers  vortieity  close  to  the 
cylinder  surface,  and  t/tpn  relea.ses  it  into  the  How  as  a  vorte.x  ’.  It  is  jdausible  that 
this  forcing  mechanism  will  continue  to  oj)erate  for  much  high(‘r  Reynolds  numbers, 
i.e.  at  least  as  high  as  the  drag-eritieal  Reynolds  number.  Re  %  >  10^.  and 

possibly  higher. 

The  results  ref)orted  here  were  retrorded  at  moderate  Reynolds  numbers,  i.e.  Re  — 
f.ox  10^.  In  an  effort  to  a.s.sess  Reynold.s-number  effects,  .some  flow  visualization  and 
wake  measurements  were  also  {)erformed  on  a  cylinder  whose  surface  was  roughened 
by  means  of  0.017  in.  diameter  glass  l)eads.  This  roughness  had  no  diseernible  effect 
on  the  wake  flow  visualization  or  velocity  measurements,  in  addition,  some  flow 
visualization  was  also  perh)rmed  on  a  smooth-surface  cylinder,  at  a  lower  Reynolds 
number  of  3.3  x  lO'*.  The  resulting  flows  «(>!•('  found  to  Ix'  viuy  similai-.  It  may  be 
interesting  to  note,  however,  that  the  effects  at  the  higher  forcing  Strouhal  numbers 
were  somewhat  less  pronounced  at  this  lower  Reynolds  nuinlxM. 

3.3.  ..-I  note  about  the  interpretation  of  linear  Mahilitif  analp.se'i 

There  has  been  much  interest  in  the  u.se  of  linear  stability  analyses  to  (h'serilx' 
unforcetl  bluff  Ixxly  flows.  See  K(X‘li  (IffSo).  Triantafyllou  &  ('hrvs.sostomidis  (IDSti). 
•Monkewitz  &  Nguyen  (1US7).  Rrovansal  ef  al,  (11187).  ('hf)maz.  Huerre  iV’  lfedekof)p 
(1988).  Monkewitz  (1988).  and  Karniadakis  &  Triantafyllou  (1989).  L'sing  linear 
theory,  it  has  been  noted  that  the  response  to  pericxlic  forcing  of  the  absolutely 
unstable  near-wake  region  of  a  bluff  Ixxly  will  be  overwhelmed  by  the  intrinsic' 
vorte.x  .shedding  nuxle  observed  in  the  abstuiee  of  external  forcing.  See  Huerre  & 
Monkewitz  (198o).  Monkewitz  &  Nguyen  (1987).  Choinaz  et  al.  (1988).  and 
Karniadakis  &  Triantafyllou  (1989).  This  argument  is  based  on  a  linearized  system 
and  suggests  that  the  pericxlic  resfx)nse  will  be  overwhelmed  because  the  transient 
n'sponse  grows  exponentially  in  tim('  while  the  forcing  is  only  fx'riodic  in  time.  Sinc<' 
th('  svst('m  cannot  lx‘  consid<'n'd  linear  for  larire  (‘xcursions  from  the  ‘  l)as('  state-. 
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ht)\\ evfi'.  tills  arj^uiiUMit  can  only  he  relied  ujioii  for  iiiHnitesiinal  foicinj:  amplitudes 
and  short  times.  It  is  not  surprising  perhaps,  that  this  interpri'tation  of  tlu‘  r(‘sult.s 
of  linear  theory  is  not  iii  accord  with  (‘xperiments  aiul  simulations  which  show  that 
finiti'-amplitude  forcing  of  tlu'  wake  of  a  hluff  body  can  result  in  a  wake  structure 
that  is  symhronous  with  the  forcing,  over  a  range  of  foicing  fiiMiueneies  and 
amplitudi's.  leg.  Koojiman  (llHiT).  Provan.'-al  <‘t  (tl.  (IhoT).  and  Karniadakis 
Triantafvllou  ( 1081)). 

i’rovansal  H  al.  (1087).  examined  th('  (nonlinear)  Landau  stahilit\  model  (Landau 
1044.  Stuart  1058.  lOhO)  and  showed  how  it  could  be  usial  to  (lualitatively  deserilx' 
the  conditions  nec-essary  for  synchronization  at  small  but  finite  forcing  am|)litudes. 
and  Reynolds  numbers  in  the  vicinity  of  th('  on.set  of  vortex  shedding.  See  also 
Landau  &  Lifshitz  (1087.  §§2b‘.  30).  Kvidently.  while  the  liiu'ar  stability  analysis,  and 
hence  the  notion  of  th('  bluff-body  wake  as  a  globally  unstable  How  with  a  region  of 
ab.solute  instability,  is  useful  in  <leseribing  the  system  dyi>amics.  given  an  a  priori 
knowledge  of  the  final  averaged  flow  (Karniadakis  &  Triantafyllou  1080).  prediction 
of  the  final  dynamic  state,  especially  in  the  presence  of  finite  external  forcing  recpiires 
a  ditten'iit  ajiproaeh. 


4.  Conclusions 

In  these  experiments,  we  examined  the  efficacy  of  oscillatory  cylinder  rotation  as 
an  actuation  mechanism  for  actively  controlling  the  cylinder  wake.  We  conclude  that 
it  is  indeed  possible  to  exercise  considerable  control  over  the  cylinder  wake  via 
oscillatory  rotary  forcing.  Not  surprisingly,  our  experiments  show  that  working  in  a 
control  domain  in  which  the  structures  shed  are  synchronous  with  the  forcing 
provides  the  greatest  control  authority  over  the  wake  structure.  In  addition,  while 
these  results  were  obtained  for  moderate  Reynolds  numbers,  i.e.  Re  =  1.5x10^, 
preliminary  evidence  suggests  that  the  description  of  the  flow  phenomena  presented 
here  is  qualitatively  the  same  over  a  large  range  of  Reynolds  numbers.  The  similarity 
of  our  results  with  tho.se  of  others  (e.g.  Koopman  1967;  Roberts  1985;  Roberts  & 
Roshko  1985;  Ongoren  &  Rockwell  1988a.  6;  Williams  &  Amato  1988.  and 
Karniadakis  &  Triantafyllou  1989).  using  various  bluff  bodies  and  forcing  techniques. 
su[)[)orts  the  projiosition  that  the  mechanisms  by  which  the  dynamics  in  the  wake 
can  be  controlled  are  largely  generie  and  have  more  to  do  with  the  ejection  of 
circulation  into  the  flow,  rather  than  with  the  behaviour  of  the  flow  ob.scrved  in  the 
absence  of  forcing. 

Finally,  we  note  that  the  type  of  active  progronnne  control  described  in  this  paper 
relies  on  the  regulated  injection  of  circulation  into  the  s(>parated  flow  region,  and  the 
subsequent  large-scale  vortical  interactions.  This  type  of  control  which  exploits  the 
natural  unsteatlinesi’  in  the  flow  instead  of  reacting  to  it  is  in  contrast  to  the  feed¬ 
forward  control  employed  by  Liepmann.  Brown  &  Xo.senehuek  (1982).  and  Liepmann 
&  Xosenchuck  ( 1982).  and  thv feedback  control  employetl  by  Ffowes- Williams  &  Zhao 
(1989).  for  example.  We  believe  that  to  control  high-Reynolds-number.  fully 
developed,  separated  turbulent  flows,  it  may  prove  the  more  efficient  and  effective 
to  programme  the  inevitable  large-.s(;ale  structures  that  dominate  the.se  flows,  rather 
than  to  attempt  to  suppress  them.  The  large  reductions  in  the  drag  coefficient  and 
the  normalized  wake  displacement  thickness  that  were  achieved  by  the  simple 
control  scheme  emjiloved  here  supports  this  ajiproach. 

We  would  like  to  thank  Hans  Hornunir.  Tonv  Li'onard.  Richard  .Minkc-Lvc*.  Herb 


Rotary  o-'^ciUation  control  of  a  cylinder  icak'r  SI) 

Gaebler.  Harry  Hamaguchi.  and  Pavei  Svitek  of  GALGIT.  as  well  as  John  Doyle  of 
the  Dept,  of  Electrical  Engineering  at  C'altech,  for  their  help  and  disi  ussions.  This 
research  was  sponsored  by  the  Air  Force  Office  of  Scientific  Hesc'arch.  TRI  AFOSK 
Grant  No.  F49620-8fi-r  0134. 


rkfkrkx(t:s 

Risiioi’.  R.  H.  I).  &  H.ass.^n.  a.  V.  1{)(34  The  lift  and  drag  on  a  cireular  cylindei  oseillating  in  a 
How  field.  Proc.  R.  Sac.  Loud.  A  277.  51-75. 

Bratt.  .1.  B.  1953  Flow  patterns  in  the  wake  of  a  oscillating  aerofoil.  Aero.  Rex.  t’ounr.  R.  &  M. 
2773. 

CnoMAZ.  .1.  M..  Hi  erre.  P.  &  Redekopp.  L.  O.  1988  Bifurcations  to  local  and  global  modes  in 
spaciallv  developing  flows.  Phys.  Rer.  Lett.  60,  25-28. 

PiMBAEA,  J.  M..  Xac.ib,  H.  M.  &  RosHKO.  A.  1988  Large  .structure  in  the  far  wakes  of  two- 
dimensional  bluff-bodies.  J.  Fluid  Mech.  190.  265-298. 

Dimot.akis.  P.  E.  1978  Laser-Doiipler  velocimetry  momentum  defect  measurements  of  cable  drag 
at  low  to  moderate  Reynolds  numbers.  XCliC  Rep.  Contract  X62583/77-M-R541 . 

Fkowcs-Wim.iams.  J.  E.  &  Zhao.  R.  C’.  1989  The  active  control  of  vorte.v  shetlding.  J .  Fluidx 
Struct.  3.  115-122. 

(IcCKENUEi.MER.  .1.  &  HoE.MES.  P.  1983  Sonlinear  OscUlations,  Dynamicut  Sydetux.  and  Bijurcationx 
of  Vector  Fieldx.  Springer. 

Hcerpe,  P.  &  Monkewitz,  P.  A.  1985  Absolute  and  ci  nveotivc  instabilities  in  free  shear  layers. 
J.  Fluid  Meek.  159.  151-168. 

Karmadakis.  (J.  E.  &  Triantapylloc.  Cl.  S.  1989  Frequency  selection  anil  asymptotic  states  in 
laminar  wakes.  ./.  Fluid  Mech.  199.  441-469. 

Kocii.  W.  1985  Local  instability  characteristics  and  frequency  determination  of  self-excited  wake 
flows.  J.  Sound  Vib.  99,  53-83. 

Koocheseahani.  M.  M.  1987  Vortical  pattern  in  the  wake  of  an  oscillating  airfoil.  2o<A  A/AA 
Aeroxpace  Sciences  Meeting,  12-1-5  January  1987,  Reno,  Nevada,  AIAA  Paper  87-0111. 

Koocheseahani.  M.  M.  &  Dimotakis.  P.  E.  1988  A  cancellation  experiment  in  a  forced  turbulent 
shear  layer.  In  Proc.  First  Natl  Fluid  Dynamics  Congr.  25-28  July  1088.  C’incinatti.  Ohio.  vol. 
II.  pp.  1204-1208.  J/.4.‘l-88-3713-CP. 

Koop.man.  (1.  H.  1967  The  vortex  wakes  of  vibrating  cylinders  at  low  Reynolds  numbers.  J.  Fluid 
Mech.  28.  501-512. 

Kcro-saka,  .m..  Christiansen,  W.  H..  (Joodman.  J.  R..  Tirres.  L.  &  Wohi.man.  R.  A.  1988 
(^ro.ssflow  trans})ort  induced  by  vortices.  J/d.-l  J.  26.  1403-1405. 

Landau.  L.  I).  1944  On  the  problem  of  turbulence.  Dokl.  Acad.  Nauk.  SSSR  44.  311-314. 

Landau,  L.  D.  &  Lifshitz.  E.  .\I.  1987  Fluid  Mechanics.  2nd  eiln.  Pergamon. 

Liepmann.  H.  W.,  Brown.  G.  L.  &  Xosenchuuk.  I).  M.  1982  Control  of  laminar-instability 
waves  using  a  new  technique.  J.  Fluid  Mech.  118.  187-200. 

Liepmann.  H.  W.  &  Xosenchuck.  I).  M.  1982  Active  control  of  laminar-turhulent  transition.  J. 
Fluid  Mech.  118.  201-204. 

Monkewitz.  P.  A.  1988  The  absolute  and  convective  nature  of  instability  in  two-dimensional 
wakes  at  low  Reynolds  numbers.  Phys.  Fluids  31.  999-l(K)6. 

Monkewitz.  P.  A.  &  Xouyen,  L.  X.  1987  Absolute  instability  in  the  near-wake  of  two- 
dimensional  bluff  bodies.  J.  Fluids  Struct.  1.  16,5-184. 

Oka.jima,  a..  Takata.  H.  &  Asanuma,  T.  1975  Viscous  flow  around  a  rotationally  oscillating 
circular  cylinder.  Inst.  Space  d-  Aero.  Sci.  (Cniversity  of  Tokyo),  Rep.  532. 

Ongoren,  a.  &  Rockweel,  D.  1988a  Flow  structure  from  an  oscillating  cylinder.  Part  I. 
Mechanisms  of  pha.,e  shift  and  recovery  in  the  near  wake.  J.  Fluid  Mech.  191.  197  223. 

Ongoren.  A.  &  Ror^wr..  D.  1988ft  Flow  structure  from  an  oscillating  cylinder.  Part  2.  Mode 
competition  in  the  ne,*r  wake.  J.  Fluid  Meek.  191,  22,5-245. 

Provansai,,  .M..  Mathis.  C.  &  Boyer.  L.  1987  Benard-von  Karman  instability:  transient  and 
forced  regimes.  ,/.  Fluid  Mech.  182,  1-22. 


!M>  7'.  J'ok'niHuni  nml  !*.  E.  I )i inoidkis 

Kohkkts.  I'.  A.  1!KS;)  Kfffctsot  a  <listiu  l)aiH-t‘ on  struct  me  and  tnixiicj:  in  tiirl>ulcnt  slicar 

label's  and  wakes.  IMi.l).  thesis.  Califoinia  Institute  of  Teehnoloey. 

lIouKKTs.  F.  A.  &  ItosiiKo.  A.  IttS.')  Fffe<-ts  of  [ifi  iodii- foieine  on  mi.xin;::  in  turiiuleiit  shear  la\  et  > 
and  wakes.  .-I/.I.-I  Slmir  Flmr  t'ltulrol  f'lmf.  li?  14  March  ntsr).  lioiihhi'.  ('().  .4/.I.I  s.”)- 

doTO. 

S'l'iWKT.  .1.  d’.  I'.ldS  ( )n  the  non-linear  meehanies  of  h\ (irod\ iiatnie  stal)ilit\  .  ./.  Fhiitl  Mn  h.  4. 
I  21. 

S'i’i  AKT.  d.d'.  liltiO  On  the  non-linear  ineel)anie>  of  wave  (list url)anees  in  stable  and  unstable 
parallel  Hows.  ./.  Fhiitl  MitIi.  9.  ‘loH 

'r.vNKD.A.  S.  HITS  \’isual  observations  of  the  tlow  past  a  circular  e\  linder  performiiii;  a  rotatoi  v 
oscillation.  •/.  Flii/'<.  Snr.  Japmi  45.  KKIS  |()4;l. 

Tki.vnt.vkvi.i.oc.  (;.  N..  'rKi.xNT.VKVi.i.ot  .  M .  S.  A-  (  iikvsso.stomidis.  ('.  litSC)  Oil  the  formation  of 
vortt'x  streets  behind  stationary  cylinders.  ./.  Fluid  Mcrh.  170.  4(51-477. 

Tkitton.  1).  .1.  Ilir)!!  Kxperiments  on  the  How  jiast  a  eireukir  cylinder  at  low  Reynolds  numbers 
•/.  Fluid  Merh.  6.  “>47  o()7. 

Whims.  I).  1072  Sfuni-infinite  vortex  trails,  and  their  relation  to  oseillatine  airfoils.  ./.  Fluid  Mnh. 
54.  (571)-()‘MI. 

W’li.l.i.x.Ms.  1).  R.  &  Am.x'I’o.  ('.  W.  IPSS  I'n.steady  pulsine  of  cylinder  wakes.  In  Finr.  First  Xull 
Fluid  l>!/u(nuirs  ('nuijr..  2->  2S  dull/  IDS.s.  Cineinatti.  Ohio.  pp.  7)51  7)57  { AI A A-SS-)5.“))52-( 

W’ll.i.i.vMsoN .  ('.  H.  K.  A  Rosmko.  .\.  HISS  X’ortex  formation  in  the  wake  of  an  oseillat iiiL: 
cylinder.  ./.  Fluids  Shuii.  2.  )5.>.“>  )5S1. 

Wf.  .1..  Mo.  .1.  &  X'akii.i.  IHSH  On  the  wake  of  a  e\  linder  with  rotational  oscillations.  .4/.4.4 
‘Jud  Shear  Flow  Couf..  Fi  IG  Mnrrh  lU^iU.  Teiuiie.  AZ.  .4/.-I.-I  Paper  80-1024. 


APPENDIX  B 


Active  Control  of  the  Flow  Past  a 
Cylinder  Executing  Rotary  Motions 


ACTIVE  CONTROL  OF  THE  FLOW  PAST 


A  CYLINDER  EXECUTING  ROTARY  MOTIONS 


Thesis  by 

Phillip  Takeo  Tokumaru 


In  Partial  Fulfillment 
of  the  Requirements  for  the  degree  of 
Doctor  of  Philosophy 


California  Institute  of  Technology 
Pasadena.  California 

1991 


(Submitted  May  31,1991) 


-  11  - 


©1991 

Phillip  Takeo  Tokumaru 
All  Rights  Reserved 


-  ni  - 


Dedicated  to  my  parents, 
Tadasu  and  Anne-Marie  Tokumaru 


-  IV  - 


Acknowledgments 

First,  I  would  like  to  thank  Mom  and  Dad,  my  sisters,  Isabelle  and  Maria, 
and  my  advisor,  Paul  Dimotaikis.  I  could  not  have  performed  this  research  without 
their  support,  patience,  and  understanding. 

I’d  also  like  to  thank  the  rest  of  my  committee,  Don  Coles,  Hans  Homung, 
Jim  Knowles,  and  Anatol  Roshko;  as  well  as  Tony  Leonard,  Dave  Dowling,  Manooch 
Koochesfahani,  Dan  Lang,  Rick  Miake-Lye,  Greg  Smedley,  and  Derek  Lisoski,  for 
lending  their  hands,  ears,  talents,  and  insights;  Herb  Gaebler,  Harry  Hamaguchi, 
and  Pavel  Svitek,  in  the  Hydro  Lab,  who  helped  me  setup  and  run  my  experiments; 
and  George  Lundgren,  George  Wilson,  Phil  Wood,  of  the  Aero  Shop,  who  helped 
me  design  and  build  my  apparatus.  Finally,  this  research  was  supported  by  the  Air 
Force  Office  of  Scientific  Research,  URI  AFOSR  Grant  No.  F49620-86-C-0134. 


Abstract 


Exploratory  experiments  have  been  performed  on  circuleir  cylinders  execut¬ 
ing  forced  rotzu-y  motions  in  a  steady  uniform  flow.  These  motions  include  harmonic 
oscillations,  steady  rotation,  2Lnd  combinations  of  the  two.  Flow  visualization  and 
laser-Doppler  velocimetry  measurements  were  used  to  characterize  the  waJce  struc¬ 
ture,  and  to  estimate  the  convection  speed,  spacing,  and  strength  of  the  vortical 
structures.  Laser-Doppler  velocimetry  measurements  were  also  made  to  estimate 
the  cylinder  drag  coefficient  and  wake  displacement  thickness.  In  addition,  the  peri¬ 
odic  flow  close  the  cylinder  and  in  the  near  waJce  region  was  mapped  for  a  particular 
forced  case.  The  data  show  that  a  considerable  amount  of  control  can  be  exerted 
over  the  flow  by  such  means.  In  particular,  a  large  increase,  or  decrease,  in  the 
resulting  displacement  thickness,  estimated  cylinder  drag,  and  associated  mixing 
with  the  free  stream  can  be  achieved,  depending  on  the  frequency  and  amplitude 
of  oscillation.  In  order  to  assess  the  effects  of  oscillatory  forcing  on  a  cylinder  with 
a  net  (mean)  rotation  rate,  a  novel  method  for  estimating  the  steady  lift  forces  was 
employed.  Using  this  method,  it  was  also  found  that  the  addition  of  forced  rotary 
oscillations  to  the  steady  rotation  of  the  cylinder  helped  to  increase  Ci  in  the  cases 
where  the  wake  would  normally  be  separated  in  the  steadily  rotating  case,  and 
decrezise  it  otherwise.  Results  obtained  for  a  steadily  rotating  cylinder  (no  forced 
oscillations)  compare  favorably  with  similar  data  published  in  the  literature. 
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=  16,  (c)  fli  =  32,  (d)  Qi  =  64 . p.  3.6 

Fig.  3.3  Location  of  figures  in  the  (5/,  )-paxameter  space . p.  3.7 

Fig.  3.4  Dual  shedding:  Forced  shedding  of  two  saune  sign  vortices  during  a  half 
forcing  cycle,  S/  =  0.2,  Dj  =2,  Re  =  1.5  x  10^ . p.  3.8 

Fig.  3.5  Global  locking:  Global  wake  structure  is  synchronized  with  the  forcing 
frequency,  (a)  S/  =  0.3,  (b)  5/  =  0.5,  (c)  S/  =  0.7,  (d)  5/  =  0.9. 
(D,  =8,  i?c  =  1.5  x  10'‘) . p.  3.9 

Fig.  3.6  Local  locking:  Local  waJce  structure  is  synchronized  with  the  forcing 
frequency,  (a)  5/  =  1.1,  (b)  5/  =  1.5,  (flj  =8,  Re  =  1.5  x  10^). 
p.  3.10 

Fig.  3.7  Shear  layer  forcing;  (a)  5/  =  2.0,  (b)  5/  =  3.3.  {Q.\  =  S,  Re  = 
1.5  X  10^) . p.  3.10 

Fig.  3.8  Natural  (unforced)  shedding.  Flow  visualized  with  dye  issuing  as  in 
Fig.  3.2  -  Fig.  3.8 . p.  3.11 
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Fig.  3.9  Diagr2im  of  (fr)o,  (^z)i ,  (^y)o,  and  {£y)i . p.  3.11 

Fig.  3.10  Vau-iation  of  initial  with  5/  (fli  w  8).  Squeires:  A^  for  the  first 
psiir  of  structures,  circles,  for  second  pair . p.  3.13 

Fig.  3.11  Variation  of  initial  Xy  with  S/  (fli  w  8).  Squares:  A^  for  the  first 
pair  of  structures,  circles,  for  second  pair . p.  3.14 

Fig.  3.12  Variation  of  initial  ratio,  XyfXi  with  5/  (Qi  ss  8).  Squares:  Ay/A. 

for  first  pair  of  structures,  and,  circles,  for  the  second  pair.  . . .  p.  3.15 

Fig.  3.13  Variation  of  A.,  Xy,  and  A^/A.  with  5/  at  x/a  =  9  with  =  8. 

Triangles:  A.,  squares:  Ay,  and  diamonds  Ay/A. . p.  3.16 

Fig.  3.14  Variation  of  A.,  Ay,  and  Ay/A.  with  5/  at  x/a  =  9  with  Q]  =  4. 

Triangles:  A.,  squares:  Ay,  and  diamonds  Ay/A. . p.  3.17 

Fig.  3.15  T^/Uooix  w.  S/.  (fii  =  8) . p.  3.19 

Fig.  3.16  Ty/UooO’  vs.  Sf.  (flj  =8) .  . p.  3.20 

Fig.  3.17  Variation  of  A* /A  with  5/  and  fii .  Dashed  lines  mark  the  region 
where  the  wake  displacement  thickness  is  less  than  in  the  unforced  case 
(note  the  minimum  around  S/  =  1  and  flj  =3).  Solid  lines  mark  the 
region  where  it  is  greater  than  or  equal  to  the  unforced  case.  .  p.  3.25 

Fig.  3.18  Cylinder  wake  mean  velocity  profiles,  measured  at  midspan  (flj  =2). 

Squares:  unforced,  circles:  5/  ~  0.2,  and  triangles:  S/  ~  1.  . .  p.  3.26 

Fig.  3.19  Cylinder  wake  displacement  thickness  variation  with  5/  (flj  =  2). 
p.  3.27 

Fig.  3.20  Drag  coefficient  \'ariation  with  5/  (fij  =  2) . p.  3.2S 

Fig.  3.21  Cd  vs.  6*/h.  Circles:  forced  (fii  =  2),  square:  unforced  (Qi  =  0). 
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p.  3.29 

Fig.  4.1  Measurement  locations . p.  4.3 

Fig.  4.2  Streamlines  averaged  over  all  phases.  Note  the  delayed  average  sepa¬ 
ration  point  and  the  small  averaged  “recirculation”  region  realized  in 
this  forcing  regime.  This  looks  very  much  like  the  (unforced)  time  aver¬ 
aged  flow  around  the  cylinder  at  higher  Reynolds  numbers  (  >  5  x  10® ), 
which  also  has  a  similar  drag  coefficient  of  about  0.2  to  0.3.  . . .  p.  4.5 

Fig.  4.3  ft  =  0.0 .  At  this  point,  the  cylinder  has  just  finished  shedding  a  vortex 
(at  x/a  =  1.5),  and  the  kernel  (bump)  of  the  next  vortex  to  be  shed 
can  be  seen  at  about  45  degrees  on  the  cylinder  surface . p.  4.8 

Fig.  4.4  ft  =  0.1.  The  bump  has  grown  larger,  note  also  the  beginnings  of  a 
second  bump  on  the  cylinder  surface  at  about  90  degrees . p.  4.9 

Fig.  4.5  ft  =  0.2.  Figure  4.5.  Both  the  main  bump  zmd  secondary  bump  have 
grown  larger,  and  moved  around  the  cylinder  in  the  clockwise  direction. 
Note  that  there  is  also  a  large  vortical  structure  being  shed  off  the  lower 
surface  at  this  time . p.  4.10 

Fig.  4.6  ft  =  0.3.  The  bumps  continue  to  grow.  Note  that  the  main  bump 
seems  to  be  held  in  place  relative  to  the  cylinder  surfaw:e,  possibly  by 
the  flow  induced  by  the  previously  shed  vortices.  Note  that  the  cylinder 
motion  is  in  the  clockwise  direction  which,  intuitively,  would  tend  to 
have  the  opposite  effect . p.  4.11 

Fig.  4.7  ft  =  0.4.  Again,  the  bumps  continue  to  grow,  but  now  are  becoming 
more  rounded.  The  large  bump  almost  looks  as  if  it  is  being  “pinched” 
into  the  flow . p.  4.12 

Fig.  4.8  ft  =  0.5.  The  presence  of  the  secondeiry  bump  is  less  noticeable;  most 
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Fig.  4.9 

Fig.  4.10 
Fig.  4.11 

Fig.  4.12 
Fig.  4.13 

Fig.  4.14 

Fig.  4.15 

Fig.  5.1 


likely  it  has  been  merged  with  the  primary  structure . p.  4.13 

ft  =  0.6 .  The  cylinder  begins  to  move  in  the  opposite  direction  to  the 
flow  on  that  side.  This  appears  to  be  when  the  vortex  begins  to  be 
released  into  the  flow . p.  4.14 

ft  =  0.7.  Note  the  beginnings  of  the  next  vortex  to  be  shed.  .  p.  4.15 

ft  =  0.8.  The  main  vortex  begins  to  move  off  to  the  right.  The  next 
vortex  to  be  shed  continues  to  be  held  in  place  both  by  the  vortices 
which  have  been  shed,  and  the  rotation  of  the  cylinder . p.  4.16 

ft  =  0.9.  The  main  vortex  continues  to  move  off  to  the  right,  p.  4.17 

Variation  of  the  normalized  circulation,  Fc/f/oca,  taken  round  contours 
of  varying  radius,  r/a  {ft  =  0.54).  The  rise  in  the  curve  between 
r/a  =  1.2  and  1.7  marks  a  patch  of  positive  (clockwise)  circulation 
and  the  drop  between  r/a  =  1.7  and  2.3  a  negative  (counter-clockwise) 
patch,  cf.  Fig.  4.14 . p.  4.19 

Streamline  pattern  associated  with  the  data  in  Fig.  4.13  {ft  =  0.54). 
p.  4.20 

Fc/f^ooti  vj.  r/a  zmd  ft.  The  daished  lines  denote  net  counter-clockwise 
rotation;  the  solid  lines,  clockwise.  The  value  Tc/U^a  at  contour  level 
‘1’  =  0.6,  ‘2’  =  0.5,  ‘3’  =  0.4  . . . ,  down  to,  T3’  =  —0.6.  Contour  level 
‘7’,  solid  line,  marks  the  zero  contour . p.  4.22 

Diagram  of  periodic  image  vortices.  Only  the  image  vortices  imme¬ 
diately  above  and  below  are  pictured.  The  shaded  circle  denotes  the 
cylinder  (not  to  scale),  and  the  smaller  circles  the  virtual  vortex  and 
its  images . p.  5.5 
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Fig.  5.2  Normalized  transverse  velocity  ahead  of  a  cylinder  executing  steady 
rotation  (fii  =  0):  2/i  v/F  vs.  (i  —  io)/a,  Rt  =  3.8  x  10^  .  flo  =  0.5  ( 
□  ),  1.5  (  O  ).  2  (  A  ),  2.5  (  +  ),  3  (  ^  ),  4  (  o  ),  5  (  *  ),  6  (  A  ), 
8  (  I  ),  10  (  —  ),  and  fli  =  0.  The  dashed  line  is  csch(7r(x  —  xo)//i)- 
p.  5.6 

Fig.  5.3  Cl  based  on  data  fit  in  Fig.  5.2  Re  =  3.8  x  10^ ,  =0,  AR  =  18.7 

(□).  Data  of  Reid  (1924)  AR  =  13.3,  Re  =  4.0  x  10*  (  Q  ),  R^  = 
5.6  X  10*  (  A  ),  Re  =  7.9  x  10*  (  +  ),  Re  =  1.2  x  10^  )•  Data 

of  Prandtl  AR  =  4.7,  Re  =  5.2  x  10*  ( - ) . p.  5.7 

Fig.  5.4  Centroid  of  the  virtual  vortex,  xq  ,  based  on  the  data  fit  in  Fig.  5.2. 
p.  5.8 

Fig.  5.5  Comparison  with  published  data  of  Cl  vs.  Qq  for  steadily  rotating 
cylinders  whose  ends  are  flush  with  the  test  section  walls.  Present  data 
(circles)  Re  =  4.2x10^  ,  (squares)  Re  =  6.8x10^,  i4R  =  18.7,  in  water. 
Data  of  Reid,  (dots)  Re  =  4.0  x  10*  <  Re  <  1.2  x  10®,  AR  =  13.3, 
in  air.  Data  of  Prandtl  (dashes).  Re  =  5.2  x  10* ,  AR  =  4.7 ,  in  air. 
p.  5.10 

Fig.  5.6  Flow  visualization  for  various  mean  rotation  rates  ( flo  )•  With  oscilla¬ 
tions  (left,  0,1  =  2.1,  Sf  =  0.7),  and  without  (right,  flj  =  0),  cf.  stars 
and  circles  in  Fig.  5.8.  (a)  Do  =  0,  (b)  Do  =  1,  (c)  Do  =  2,  (d)  Do  = 
3,  (e)  Do  =  4,  (f)  Do  =  5.  Re  =  4.1  x  10® . p.  5.12 

Fig.  5.7  Comparison  of  Cl  vs.  Do  data  for  severed  amplitudes  of  oscillations, 
with  constant  forcing  frequency  and  Reynolds  number  {Sf  =  0.7 ,  Re  = 
6.8  X  10®).  Squares:  Di  =  0,  crosses;  1.2,  and  diamonds:  2.3.  p.  5.11 

Fig.  5.8  Comparison  of  Cl  vs.  Dq  data  for  the  steady  case  and  an  oscillating 
case,  at  two  Reynolds  numbers.  Steady  rotation:  (circles)  Re  =  4.2  x 
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10^,  ill  =  0,  (squares)  Re  =  6.8  x  10^,  Sli  =  0.  Net  rotation  with 


oscillations,  (stars)  Re  =  4.2  x  10^,  ili  =  2.1,  S/  =  0.7,  (diamonds) 
Re  =  6.8  X  10=*,  fii  =  2.3,  5/  =  0.7 . p.  5.13 

Fig.  5.9  The  effect  on  Cl  of  var3ring  S/.  Circles:  ilo  =  0.3,  squares:  ilo  =  1-4 
{Re  =  6.8  X  10^  ill  =2.3) . p.  5.15 

Fig.  A.l  3-D  Geometry.  . p.  A. 2 

Fig.  A.2  2-D  Geometry.  . p.  A.4 

Fig.  B.l  Time-averaged  streamlines . p.  B.l 


Fig.  B.2  Streamlines  for  flow  averaged  at  constant  phase.  Sequence  on  left: 

0.00  <  ft  <  0.08,  right:  0.10  <  ft  <  0.18 . p.  B.2 

Fig.  B.3  Streamlines  for  flow  averaged  at  constant  phase.  Sequence  on  left: 

0.20  <  ft  <  0.28,  right:  0.30  <  ft  <  0.38 . p.  B.3 

Fig.  B.4  Streamlines  for  flow  averaged  at  constant  phase.  Sequence  on  left: 

0.40  <  ft  <  0.48,  right:  0.50  <  ft  <  0.58 . p.  B.4 

Fig.  B.5  Streeunlines  for  flow  averaged  at  constant  phase.  Sequence  on  left: 

0.60  <  ft  <  0.68,  right:  0.70  <  ft  <  0.78 . p.  B.5 

Fig.  B.6  Streamlines  for  flow  averaged  at  constant  phase.  Sequence  on  left: 

0.80  <  ft  <  0.88,  right:  0.90  <  ft  <  0.98 . p.  B.6 


CHAPTER  1 


Introduction 


1.1  Preface 

The  title  of  the  thesis,  “Active  Control  of  the  Flow  Past  a  Cylinder  Executing 
Rotary  Motions,”  can  mean  different  things  to  different  people.  In  this  case,  it  refers 
to  the  flow  past  a  circular  cylinder  in  a  water  tunnel.  This  cylinder  was  bounded 
at  both  ends  and  below  by  test  section  walls,  and  from  above  by  the  water  tunnel 
free  surface.  See  Fig.  1.1.  To  control  the  structure  of  the  flow,  the  cylinder  was 
programmed  to  execute  rotary  motions  about  its  axis.  This  control  strategy  is  called 
“active,”  because  energy  must  be  (actively)  supplied  to  execute  these  motions,  cf. 
“passive”  controls,  such  as  spoilers  or  riblets,  which  require  no  special  input  of 
energy.  Since  no  feedback  was  necessary  (or  used)  the  forcing  was  under  “open- 
loop,”  or  “program”  control. 

In  the  past  this  has  been  called  “forcing  the  flow,”  and  if  the  flow  responded 
in  a  periodic  way,  terms  like  “resonance,”  “locking,”  and  “synchronization”  were 
used.  The  notion  of  controlling  unsteadiness  was  reserved  for  the  elimination  of 
instabilities.  The  aim  of  the  present  investigations,  however,  is  to  program  the  in¬ 
evitable  vortical  structures  that  dominate  unsteady  flow  behind  a  circular  cylinder. 


Fig.  1.1  Configuration 


1.2  Background 

Traditionally,  experiments  involving  a  cylinder  executing  rotary  motions  in 
a  uniform  flow  have  been  performed  with  cylinder  motions  that  are  either  harmonic 
oscillations  or  steady  rotation.  It  is  believed  that  the  present  investigations  are  the 
first  to  examine  cylinder  motions  combining  the  two.  Because  of  this,  background 
material  for  the  oscillating  and  steadily  rotating  cases  falls  naturally  into  two  parts 
presented  here  as  Secs.  1.2.2  “Rotary  oscillations”  and  1.2.3  “Steady  rotation.”  Sec¬ 
tion  1.2.1  introduces  a  class  of  forcing  functions  which  includes  (rotary)  oscillations 
plus  net  rotation  of  the  cylinder. 


-  1.3- 


1.2.1  Forcing  parameters 

The  motion  of  the  cylinder  in  the  present  experiments  can  be  described  by 
a  single  equation;  i.e., 

=  Qo  +  fii  sin(27r/<)  ,  (1.2.1) 

where 

a  =  p-.  (1.2.2) 

'-'oo 

is  the  normalized  rotation  rate  of  the  cylinder,  a  is  the  cylinder  radius,  Uoo  is  the 
free-stream  velocity,  t  is  the  time,  6  is  the  angular  velocity  of  the  cylinder,  /  is 
the  forcing  frequency,  and,  flo  and  fli  are  amplitudes  of  the  steady  and  harmonic 
components  of  the  cylinder  motion.  The  normalized  forcing  frequency  is 

S,  s  ^  ,  (1.2.3) 

^OO 

t.e.,  the  forcing  Strouhal  number. 


Three  forcing  parameters  are  considered  in  the  present  experiments,  flo  > 
zmd  S/.  In  addition,  another  important  parameter  is  the  Reynolds  number, 


Re  = 


2aUo 


(1.2.4) 


where  v  is  the  kinematic  viscosity. 


Another  dimensionless  quantity  which  might  have  been  used  is  the  angular 
position  of  the  cylinder,  0{t).  In  the  present  experiments,  however,  fl(<)  was  used 
in  favor  of  ${t)  to  provide  consistent  notation  for  the  normalized  mean  ( flo )  and 
harmonic  ( fli )  components  of  the  rotation  rate.  Using  angular  position  to  describe 
the  motion  of  a  cylinder  with  a  net  rotation  rate  is  awkward.  Note  that  with  fio  =  0 
(no  net  rotation  rate),  5/  =  1,  and  fli  =  2,  the  amplitude  of  the  oscillations  in  the 
angular  position  of  the  cylinder  is  6i  s;  36.5®  . 
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1.2.2  Rotary  oscillations 

The  dynamic  behavior  of  the  flow  over  a  bluff  body  is  dominated  by  large 
scale  vortical  structures  shed  in  its  wake.  Weihs  (1972)  proposed  a  methodology 
for  predicting  the  structure  of  such  a  wake  and  found  that  it  could  become  wider, 
narrower,  or  remain  the  same,  depending  on  the  initial  strength  and  spacing  of  these 
vortices.  His  work  was  primarily  aimed  at  describing  the  wake  behind  an  oscillating 
airfoil  (e.^^.,  Bratt  1953,  Koochesfahani  1987).  Nevertheless,  Weihs’s  results  do 
suggest  the  potential  for  effective  control  of  wakes  in  general. 

Okajima  et  al.  (1975)  examined  the  forces  acting  on  a  rotationally  oscillating 
cylinder,  for  Reynolds  numbers  based  on  cylinder  diameter  in  the  range  of  40  < 
Re  <  6  X  10* .  Their  measurements  were  for  a  normalized  peak  rotation  rate  in 
the  range  0.2  <  flj  <  1.0,  and  a  forcing  Strouhal  number  in  the  range  0.05  < 
Sf  <  0.3.  They  noted  a  “synchronization”  similar  to  that  observed,  for  example, 
by  Bishop  &  Hassan  (1964),  Koopman  (1967),  and,  more  recently,  Ongoren  & 
Rockwell  (1988a,b),  and  Williamson  &  Roshko  (1988)  for  a  cylinder  in  transverse 
and  in-line  oscillation.  Both  Okajima  and  Bishop  &  Hassan  reported  a  hysteresis 
in  their  measurements  for  increasing  and  decreasing  S / .  Investigations  at  similarly 
low  normalized  forcing  frequencies,  amplitudes  and  Reynolds  numbers  were  also 
performed  by  Wu  et  al,  (1989).  Taneda  (1978)  demonstrated  that  in  the  range 
30  <  Re  <  300,  the  “dead  water”  region  behind  a  cylinder  can  be  removed  for 
sufficiently  large  (fli  >  7  to  27,  depending  on  such  factors  as  the  cylinder  span 
and  tank  width).  For  somewhat  larger  fli  vortex  shedding  was  inhibited.  Similar 
behavior  was  also  documented  by  Williams  &  Amato  (1988)  using  a  line  of  unsteady 
pulsing  jets  embedded  in  the  trailing  edge  of  the  cylinder,  at  a  comparable  Reynolds 
number  of  370. 
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The  present  experiments  covered  in  Chs.  3  and  4  were  performed  on  circular 
cylinders  executing  forced  rotary  oscillations  in  a  uniform  flow  (no  mean  rotation 
rate).  The  Reynolds  numbers  are  higher  than  in  previous  investigations  at  Re  = 
1.5  X  10^  in  Ch.  3  and  Re  =  2.1  x  10^  in  Ch.  4.  Okajima  et  al.  approached  these 
Reynolds  numbers  ( <  6  x  10^ ),  but  their  experiments  examined  perturbations  of 
the  flow  employing  peak  forcing  amplitudes,  ,  sm  order  of  magnitude  lower  than 
in  the  present  investigations. 

In  Ch.  3,  flow  visualization  data  were  used  to  characterize  the  wake  structure, 
and  to  estimate  the  convection  speed,  spacing,  and  strength  of  the  vortical  struc¬ 
tures.  In  addition,  LDV  measurements  were  made  to  estimate  the  cylinder  drag 
coefficient,  and  wake  displacement  thickness.  The  data  show  that  a  considerable 
amount  of  control  can  be  exerted  over  the  wake  flow  by  such  means.  In  addition, 
a  large  increase,  or  decrease,  in  the  associated  displacement  thickness,  estimated 
cylinder  drag,  and  mixing  with  the  free  stream  was  achieved,  depending  on  the 
frequency  and  amplitude  of  oscillation.  See  also  Tokumaru  &  Dimotakis  (1991). 

Chapter  4  examines  a  particular  case  of  a  cylinder  executing  forced  rotary 
oscillations.  An  estimate  of  the  unsteady,  phase-averaged  streamfunction  was  made. 
It  was  found  that  the  forced  shedding  mechanism  gathers  vorticity  close  to  the 
cylinder  surface  and  then  releases  it  into  the  flow.  This  is  in  contrast  to  natural 
(unforced)  shedding,  where  the  vorticity  is  released  into  the  flow  from  more-or-less 
steady  separation  points,  allowing  the  vortices  to  form  in  the  wake.  The  speed 
and  strength  of  the  forced  vortical  structures  was  estimated  and  it  was  found  that 
very  little  of  the  asymmetry  (quamtified  by  the  circulation  generated  at  the  cylinder 
surface)  is  found  in  the  wake. 
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1.2.3  Steady  rotation 

A  rotating  cylinder  moving  in  a  uniform  stream  experience  a  force  normal 
to  the  direction  of  motion.  Goldstein  (1938)  makes  several  historical  references  for 
both  rotating  spheres  and  cylinders  and  credits  Magnus  (1853)  with  the  first  labo¬ 
ratory  experiments  examining  the  lift  on  a  rotating  cylinder.  Early  in  this  century, 
experiments  on  a  circular  cylinder  rotating  about  its  axis  in  a  uniform  flow  were 
performed  by  Reid  (1924),  Prandtl  (1925),  Thom  (1926,  1931),  for  example.  Now, 
on  the  verge  of  the  21**  century  researchers  continue  to  find  this  flow  interesting, 
t.g.,  the  experiments  and  simulations  of  Badar,  tt  al.  (1990)  for  the  flow  past  a 
cylinder  impulsively  started  in  both  rotation  and  translation. 

It  can  be  shown  that  the  mean  lift  coefficient  of  a  2-D  body  can  be  written 
as 


where  p  is  the  fluid  density,  U^o  is  the  free-stream  velocity,  a  is  the  cylinder  radius, 
L  is  the  lift  per  unit  span,  and  F  is  the  circulation  taken  round  a  contour  enclosing 
the  lifting  body.  See  for  example  Taylor  (1925),  Thwaites  (1960,  §  V.9).  That 
this  circulation  could  be  measured  around  contours  close  to  the  cylinder  was  shown 
experimentally  by  Thom  (1931). 

In  the  potentizd  flow  prototype  of  a  rotating  cylinder  in  a  uniform  free  stream, 
the  rotation  of  the  cylinder  is  modeled  by  placing  a  point  vortex  of  strength  F  at 
its  center.  When  modeling  a  flow,  F  should  be  considered  as  a  parameter  to  be  ad¬ 
justed  so  that  the  flow  is  properly  represented;  e.g.,  to  satisfy  the  Kutta  Condition. 
Proposing  that  Cl  can  be  made  arbitrarily  large  by  arbitrarily  increeising  F  is  not 
very  helpful;  F  is  the  mean  lift;  cf.  Eq.  1.2.5. 
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Pr2uidtl  (1925)  argued  that  the  maximum  circulation,  Tmix ,  which  could 
be  reedized  about  the  rotating  cylinder  was  equal  to  the  circulation  at  which  the 
upstream  and  downstream  stagnation  points  join  on  the  bottom  of  the  cylinder. 
From  the  potential  flow  analog  of  that  case,  a  value  of 

Fniix  =  ^irUooO  (1.2.6) 

is  used.  See  also  the  discussion  in  Goldstein  (1938  I,  pp.  81-84).  Beyond  this  point, 
Prandtl  argued  that  no  circulation  could  be  shed  by  the  cylinder  to  infinity  as  it 
is  rsimped  up  from  rest  to  some  constant  rotation  rate  of  flo  >  4 ,  fixing  the  total 
circulation  at  infinity  and  the  lift  of  the  rotating  cylinder  thereby.  If  this  argument 
were  correct,  the  maximum  (steady  state)  lift  coefficient  that  can  be  realized  would 
be  given  by, 

=  4,r  «  12.6  .  (1,2.7) 

It  is  not  at  all  clear  from  the  discussions  in  the  literature,  but  Prandtl  relates  the 
real  and  potential  flow  case  of  the  rotating  cylinder,  not  by  matching  the  circulation 
calculated  by  taking  a  contour  round  the  cylinder  at  the  surface,  as  one  would 
expect,  but  by  matching  the  peak  circumferential  velocities  in  the  two  cases.  Not 
coincidently,  this  is  also  the  choice  that  “looks  right.” 

Chapter  5  examines  the  effect  of  cylinder  rotation  on  the  flow  ahead  of  the 
cylinder.  In  order  to  assess  the  effects  of  forcing,  a  novel  method  for  describing 
the  flow  ahead  of  the  cylinder,  in  terms  of  the  strength  and  position  of  a  virtual 
vortex,  was  employed.  The  strength  of  the  virtual  vortex  was  then  related  to  the 
lift  coefficient  of  the  cylinder.  Results  obtained  for  a  steadily  rotating  cylinder  (no 
forced  oscillations)  compare  favorably  with  simileir  data  published  in  the  literature. 
It  was  also  found  that  the  addition  of  forced  rotary  oscillations  to  the  steady  rotation 
of  the  cylinder  helped  to  increase  the  lift  coefficient  in  the  cases  where  the  wake 
would  normally  be  separated  in  the  steadily  rotating  caise,  and  decrease  it  otherwise. 
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It  is  believed  that  this  is  the  first  investigation  of  the  flow  past  a  cylinder  with  both 
net  rotation  and  oscillations. 


1.3  A  note  on  the  transport  of  vorticity  about  the  cylinder 

The  concept  of  separation  is  useful  because  it  gives  an  intuitive  feel  for  the 
transport  of  vorticity  into  a  flow.  For  the  generic  case  of  flow  about  a  bluff  body 
in  steady  uniform  flow,  the  vorticity  is  shed  into  the  wake  from  “slowly”  moving 
separation  points.  It  might  be  appropriate  to  cadi  this  kind  of  separation  quasi- 
sitady.  In  contrast,  when  a  body  is  accelerated  in  rotation,  or  translation,  as  in  the 
present  experiments,  separations  can  appear,  move,  and  disappear  over  length  and 
time  scales  comparable  those  of  the  body  motion  and  the  surrounding  flow.  In  the 
latter  case  it  is  even  difficult  to  define  “separation” . 

r 


Uoo 


Fig.  1.2  Flow  past  a  cylinder. 


An  ad  hoc  approach  is  to  examine  the  trzuisport  of  vorticity  close  to  the  body. 
In  particuleur,  for  the  ceise  of  a  circular  cylinder,  the  convective  vorticity  flux  across 
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a,  6  =  constant  ray  (traversing  the  boundary  layer)  can  be  written  as 

J  =  I  Ljug  dr  , 

where  6  is  the  boundary  layer  thickness.  In  dimensionless  form 

rs/a 


iVl  I  Ul  “'(a)  • 


If  the  vorticity  in  the  boimdary  layer  can  be  approximated  by 

dug 


(jJ 


dr 


then  Eq.  1.3.2  can  be  approximated  by 


M'^)\ 
Uoo  J 


(1.3.1) 


(1.3.2) 


(1.3.3) 


(1.3.4) 


Looking  at  the  vorticity  convected  in  the  boundary  layer  at  9  =  90“  gives 
an  indication  of  the  vorticity  convected  into  the  (dead  water)  region  behind  the 
cylinder.  See  Fig.  1.2. 


Consider  the  czise  of  a  cylinder  executing  rotary  oscillations  with  no  mean 
rotation  rate.  The  potential  flow  model  for  the  flow  past  a  circular  cylinder  has 


ug  =  -2Uoo  at  90°  , 


(1.3.5) 


which  is  close  to  the  velocity  outside  the  boundarj'^  layer  in  the  resJ  case.  Using 
Eqs.  1.3.5,  1.2.2,  and  1.2.1  with  fto  =  0  (no  meem  rotation  rate)  then  yields  for  the 
vorticity  convected  past  the  90“  point. 


.  4  -  Q? 


11/2 

2^00 


Ilf  sin^{2n ft)  . 


(1.3.6) 
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Hence  any  rotaxy  oscillation  of  the  cylinder  will  tend  to  decrease  the  convection  of 
vorticity  past  the  90“  point.  Note  that  when 

n?  sin^(27r/<)  =  4  ,  (1.3.7) 

(instantaneously)  no  vorticity  is  convected  past  the  90“  point  as  would  be  the  case 
in  potential  flow.  Equation  1.3.7  suggests  that  a  value  of 

=  2  (1.3.8) 

may  be  needed  to  effectively  control  this  flow. 

If,  on  average,  no  vorticity  is  to  be  convected  past  the  90“  point,  then 
will  satisfy 

fti  w  2v/2  «  2.8  .  (1.3.9) 

It  is  shown  (experimentally)  in  Ch.  3  that  (5/,fii)  «  (1,3)  marks  the  minimum 
wake  displacement  thickness  in  the  (5/,fIj  )-parzuneter  space. 


CHAPTER  2 


Experimental  facilities 


2.1  Flow  facility 

The  initial  phase  of  the  current  experiments,  ».c.,  W2ike  mean  velocity  mea¬ 
surements  and  preliminary  flow  visualization,  was  carried  out  in  the  18"  wide  by  19" 
deep  Low  Speed  Water  Channel  (LSWC)  located  in  the  basement  of  Guggenheim.  * 
Subsequent  flow  visuzJization  and  lift  measurements  were  preformed  in  the  20"  x  20" 
Free  Surface  Water  Tunnel  (FSWT)  in  the  Hydrodynamics  Lab.  See  Ward  (1976). 


2.2  The  cylinder  apparatus 

The  cylinder  support  structure,  Fig.  2.1,  holds  a  cylinder  10"  above  the 
bottom  of  the  water  channel,  between  0.5"  thick  Plexiglas  fairings  placed  flush  to 
the  side  walls  of  the  channel.  Power  is  transferred  from  the  motor  to  the  cylinder 
using  steel  reinforced  neoprene  timing  belts.  They  are  entirely  enclosed  in  the 
fairings  zmd  do  not  interfere  with  the  flow,  cf.  Figs.  2.1  and  2.2. 

*  Also  known  as  “The  Student  Channel.” 
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When  choosing  a  motor  and  cylinder  diameter,  the  minimum  design  goal  was 
to  be  able  to  oscillate  the  cylinder  with  peak  rotational  rates  greater  than  Qi  =  2 
(as  discussed  in  Sec.  1.3)  and  at  Reynolds  numbers  around  Re  =  10^.  It  was  felt 
that  frequencies  an  order  of  magnitude  faster  than  the  natural  (unforced)  shedding 
frequency  would  be  fast  enough  ( 5/  ss  2 ).  In  addition,  a  five  times  “engineering” 
factor  was  used. 

The  cylinders  were  driven  by  a  high  performance  JR16M4CH  ServoDisc‘"* 
DC  motor,  from  PMI  Motion  Technologies,  capable  of  tracking  an  arbitrary  com¬ 
mand  sign2Ll,  within  the  bandwidth  and  slewing  rates  that  were  investigated.  This 
command  signal  was  generated  by  a  function  generator,  or  a  computer  digital-to- 
analog  converter  output  channel.  In  addition,  the  angulaur  position  of  the  cylinder 
could  be  monitored  through  a  13-bit  M25  Series  absolute  position  encoder  from  BEI 
Motion  Systems  Co. 

Cylinders  with  a  radius  rzmging  from  0.25”  to  2”  and  span  from  17”  to 
19”  were  used.  Cylinders  were  made  out  of  both  anodized  aluminum  and  Plexiglas. 
Two  2”  radius  cylinders  were  used  in  the  present  experiments.  One  was  made  of 
smooth  anodized  2Juminum  tubing.  The  other  wzis  made  of  machined  Plexiglas 
tubing  equipped  with  ports  so  that  dye  could  be  introduced  into  the  flow  from  the 
cylinder  surface.  Additional  flow  visualization  was  accomplished  by  introducing  dye 
upstream  of  the  cylinder,  through  hypodermic  tubing. 


Fig.  2.1  The  cylinder  apparatus. 


Fig.  2.2  The  cylinder  drive  mechanism. 


-2.5- 


2.3  Measurement  apparatus 

Measurements  were  made  using  a  number  of  laser  Doppler  velocimetry  sys¬ 
tems.  More  details  are  given  in  the  “experimental  setup”  sections  for  each  chapter. 
In  the  LSWC  a  u -channel  Bragg  cell  frequency-shifted  He  Ne  based  system  was 
used.  Translation  of  the  LDV  (and  measurement  volume)  was  effected  by  a  manusd 
x-y  traverse  system. 

In  the  FSWT  u,  and  vg,  Argon-ion  based  systems,  derived  from  the 

LDV  system  of  Lang  (1985)  were  used.  All  but  the  u  in  the  {u,v)  arrangement, 
were  Bragg  cell  frequency-shifted.  Translation  of  the  LDV  was  performed  by  a 
digital  x-y-z  traverse,  and  rotation  of  the  focal  volume  for  the  vg  measurements 
was  effected  by  a  rotatable  dove  prism. 

The  data  were  gathered  using  LSl-11  data  acquisition  computers,  and  then 
transferred  to  DEC  VAX,  Apple  Macintosh,  and  NeXT  NeXTstation  computers  for 
further  reduction  emd  presentation. 


CHAPTER  3 


Non-lifting  cylinder 


3.1  Preface 

This  chapter  describes  exploratory  experiments,  performed  on  circular  cylin¬ 
ders  executing  rotary  oscillations,  in  a  steady  uniform  flow.  Flow  visualization  data 
were  used  to  characterize  the  wake  structure  and  to  estimate  the  convection  speed, 
spacing,  and  strength  of  the  vortical  structures.  In  addition,  LDV  measurements 
were  made  to  estimate  the  cylinder  drag  coefficient  and  wake  displacement  thick¬ 
ness.  The  data  show  that  a  considerable  amount  of  control  can  be  exerted  over  the 
flow  by  such  means.  In  particulzu:,  a  large  increase,  or  decrease,  in  the  resulting  dis¬ 
placement  thickness,  estimated  cylinder  drag,  and  associated  mixing  with  the  free 
stream  can  be  achieved,  depending  on  the  frequency  and  amplitude  of  oscillation. 
See  also  Tokumairu  &  Dimotakis  (1991). 
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3.2  Experimental  setup 

For  the  experiments  detailed  in  this  chapter,  fto  was  set  to  zero  and  was 
chosen  such  that  within  the  range  of  parameters  examined,  the  peak  circumferen¬ 
tial  velocity  of  the  cylinder  would  be  comparable  to  the  velocity  just  outside  the 
boimdary  layer  of  the  cylinder  (approximately  twice  the  free-stream  velocity).  It 
was  anticipated  that  amplitudes  of  this  magnitude,  or  greater,  would  be  necessary 
to  effect  a  significant  change  in  the  flow,  as  was  discussed  in  Sec.  1.3. 

The  experiments  presented  in  Secs.  3.3  and  3.4  examine  the  control  parame¬ 
ters  fli  and  5/  over  a  range  of  values.  The  free-stream  velocity  was  approximately 
15  cm/s,  yielding  a  Reynolds  number  based  on  cylinder  diameter  of  Re  =  1.5  x  10^ . 
In  addition,  wake  mean  and  rms  velocity  profiles  were  measured  at  a  single  stream- 
wise  station,  located  4.5  diameters  downstream  of  the  cylinder  axis.  The  remaining 
sections  in  this  chapter  touch  on  some  other  issues  pertaining  to  the  present  exper¬ 
iments. 

The  wake  streamwise  mean  velocity  measurements  were  performed  in  the 
GALCIT  18"  wide  by  19"  deep  Low  Speed  Water  Channel,  using  a  He-Ne  based 
laser  Doppler  velocimeter  (LDV).  The  cylinder  used  for  the  wake  velocity  meaisure- 
ments  is  an  anodized  and  machined  aluminum  tube,  4"  in  diameter,  and  mounted 
and  driven  by  the  apparatus  described  in  Ch.  2.  The  command  signal  for  the 
cylinder  velocity  was  generated  by  a  computer  controlled  function  generator. 

The  (15  mW  He-Ne)  laser  beaun  for  the  LDV  velocity  measurements  was  split 
into  two  beams  of  nearly  equal  intensity  using  a  cube  beam  splitter.  The  beams 
were  then  Bragg  cell  frequency-shifted  to  accommodate  flow  reversals.  The  Bragg 
cell  offset  between  the  two  beams  was  100  to  200  kHz  for  the  data  presented  in 
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this  chapter.  An  achromat  was  used  to  focus  the  beams  in  a  focal  (measurement) 
volume  located  in  the  midspan  plane  of  the  LSWC.  The  light  scattered  from  the 
focal  volume  was  gathered  by  an  achromat  and  focused  on  a  photodiode  with  an 
integral  low  noise  pre-amplifier  designed  by  Dan  Lang.  The  signal  from  the  pre¬ 
amplifier  was  then  band-pass  filtered  above  and  below  the  frequency  corresponding 
to  the  mean  flow  velocity.  The  band-pass  range  was  chosen  to  pass  instantsmeous 
velocity  fluctuations.  A  tracking  phase- locked  loop  designed  by  Dan  Lang  and  Paul 
Dimotakis  was  then  used  to  lock  a  TTL  square  wave  to  the  dominant  frequency.  The 
TTL  signal  was  then  read  by  a  counter-timer  board  on  a  data  acquisition  computer 
zmd  stored  on  disk  for  later  processing.  The  transmitting  and  receiving  hardware 
was  mounted  on  a  manual  x~z  traverse  (positioner). 

Photographs  were  taken  in  the  G  ALCIT  20"  x  20"  Free  Surface  Water  Tunnel 
(FSWT)  using  a  35mm  camera.  Red  dye  was  used  in  all  figures  except  for  Fig.  3.1. 
A  blue  filter  was  used  (except  for  Fig.  3.1)  to  darken  the  red  dye  marker  relative  to 
the  white  background.  Illumination  was  provided  from  behind  by  a  Izurge  (backlit) 
white  sheet  of  paper,  suspended  outside  of  the  tunnel,  behind  the  model.  The 
cylinder  used  for  the  flow  visualization  data  is  a  machined  Plexiglas  tube,  4"  in 
diameter.  Holes  for  introducing  dye  into  the  flow  were  drilled  through  the  surface 
at  midspan.  Additional  flow  visualization  was  accomplished  by  introducing  dye 
through  hypodermic  tubing  into  the  flow  aihead  of  the  cylinder.  The  cylinder  wzis 
mounted  and  driven  by  the  apparatus  described  in  Ch.  2.  The  command  signal  for 
the  cylinder  velocity  was  generated  by  a  computer- controlled  function  generator. 
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3.3  Flow  visualization 

A  useful  picture  of  the  control  that  can  be  exercised  over  the  separated 
flow  in  the  wake  of  the  cylinder  emerges  from  the  flow  visualization  data.  These 
demonstrate  the  overall  narrowing,  or  widening,  as  well  as  the  formation  of  distinct 
vortical  structures  in  the  wake,  depending  on  the  values  of  the  control  parameters 
5/  and  fli . 


In  Fig.  3.1,  dye  is  introduced  ahead  of  the  cylinder  by  a  comb  of  dye  injectors. 
In  the  unforced  case  (Fig.  3.1,  top),  the  dye  can  be  seen  to  be  dispersed  and  mixed 
across  the  full  height  of  the  test  section,  within  a  few  dizuneters  downstreaun  of  the 
cylinder.  In  contrast,  the  dye  marker  in  the  forced  case,  corresponding  to  a  reduced 
wake  width  (Fig.  3.1,  bottom),  occupies  approximately  the  same  fraction  of  the  test 
section  height  on  exit  as  it  did  entering. 

The  data  in  Fig.  3.2  and  Fig.  3.5  (recorded  with  dye  issuing  from  the  surface 
of  the  cylinder)  illustrate  how  the  wake  may  be  made  wider,  or  narrower,  depending 
on  the  control  parameters  5/  and  f2i .  Weihs  (1972)  examined  this  phenomenon 
by  considering  the  mutual  influence  of  the  vortices  in  the  wake  amd  those  being 
shed.  Although  his  analysis  is  highly  idealized,  considering  only  point  vortices  in 
potential  flow,  it  does  provide  a  physical  argument  for  the  various  trajectories  taken 
by  the  vortices.  Additional  flow  visualization  data  (also  recorded  with  dye  issuing 
from  the  surface  of  the  cylinder),  at  fixed  fli  =  8  and  increaising  5/ ,  are  depicted 
in  Fig.  3.6  and  Fig.  3.7. 
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Fig.  3.2  Forced  shedding  at  fixed  5/  =  1  {Re  =  1.5  x  10^).  (a)  fii  =  8,  (b) 
ft,  =  16,  (c)  ft,  =  32,  (d)  ft,  =  64. 

3.3.1  Wake  structure 

Several,  qualitatively  different,  vortex  shedding  modes  were  observed.  They 
are  presented  here  in  order  of  increasing  frequency.  Figure  3.3  locates  the  various 
modes  in  the  ( 5/, ft,  )-paranneter  space.  In  the  dual  mode  (e.^..  Fig.  3.4),  the 
cylinder  releases  two  vortices  of  the  same  sign  per  half-cycle.  This  mode  was  not 
easy  to  set  up  and  should  not  be  considered  robust  enough  for  use  in  active  program 


control. 
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S  f  =  I  flow  visualization  data 


Fig.  3.3  Location  of  figures  in  the  )-paraineter  space. 


In  the  global  mode  {e.g.,  Fig.  3.5),  the  wake  structure  is  synchronized  with 
the  forced  cylinder  oscillation,  auid  persists  beyond  the  end  of  the  test  section. 
This  mode  was  relatively  easy  to  set  up  because  it  occurred  over  a  wide  range  of 
forcing  parameters.  In  the  local  mode  (e.^..  Fig.  3.6),  the  near  wake  structure  is 
synchronized,  but  becomes  unstable  and  evolves  into  a  structure  with  lower  spatial 
frequency  some  distance  downstream  of  the  cylinder.  In  the  shear  layer  mode  (e.^.. 
Fig.  3.7),  the  wake  structure  is  not  synchronized  with  the  forcing.  The  effect  of 
forcing  is  primarily  observed  in  the  shear  layers  sepeu-ating  from  the  cylinder.  Note 
that  while  the  flow  visualization  photographs  of  the  forced  cylinder  wake  at  the 
largest  values  of  Sf  (cf.  Fig.  3.7,  bottom)  resemble  those  for  the  unforced  case  {cf. 
Fig.  3.8),  it  was  found  that  the  wake  displacement  thickness  (Eq.  3.4.2)  in  these 
forced  cases  is  still  noticeably  less  than  in  the  unforced  case. 
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Fig.  3.4  Dual  shedding:  Forced  shedding  of  two  same  sign  vortices  during  a  half 
forcing  cycle,  5/  =  0.2,  Di  =2,  ile  =  1.5  x  10^ . 

The  transition  between  the  global  and  shear  layer  modes  occurs  gradually 
through  the  local  mode.  Figures  3.5  to  3.7  demonstrate  how  the  “break-up”  of  the 
synchronized  structure  in  the  local  mode  moves  toward  the  cylinder  with  increasing 
5/.  The  shear  layer  mode  is  then  found,  at  even  higher  5/,  when  none  of  the 
wake  structure  (save  small  perturbations  of  the  separating  shear  layers)  is  shed  syn¬ 
chronously  with  the  forcing.  This  “break-up”  in  the  flow  structure  in  the  local  and 
sheau'  layer  modes  may  be  attributable  to  the  stability  characteristics  of  the  evolving 
mean  velocity  profiles  (as  suggested  by  Cimbala,  Nagib  &  Roshko  1988,  for  lower 
Reynolds  number  flows,  far  downstream),  vortex  coalescence,  three-dimensionail  ef¬ 
fects,  or  some  combination  thereof. 

Behavior  similar  to  the  dual,  global,  and  local  modes  was  previously  observed 
by  Roberts  (1985)  and  Roberts  &  Roshko  (1985)  for  the  caise  of  a  forced  walce  behind 
a  splitter  plate,  and  Williamson  &  Roshko  (1988)  for  a  cylinder  executing  tremsverse 
oscillations  in  a  free  stream,  for  example. 
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Fig.  3.5  Global  locking;  Global  wake  structure  is  synchronized  with  the  forcing 
frequency,  (a)  5/  =  0.3,  (b)  5/  =  0.5,  (c)  5/  =  0.7,  (d)  5/  =  0.9. 
(n,  =  8,  fie  =  1.5  X  10^). 

3.3.2  Speed  and  spacing  of  the  vortical  structures 


Flow  visualization  techniques  usually  considered  only  for  qualitative  anal¬ 
ysis  of  flow  fields.  In  the  present  experiments,  however,  estimation  of  the  normalized 
speed  (celerity)  of  the  vortical  structures  relative  to  the  cylinder,  Uy  =  fix  ,  divided 
by  the  free-stream  velocity. 


A, 


(3.3.1) 


Fig.  3.6  Local  locking:  Local  wake  structure  is  synchronized  with  the  forcing  fre¬ 
quency.  (a)  5/  =  1.1,  (b)  5/  =  1.5,  (fti  =8,  Re  =  1.5  x  10^). 


Fig.  3.7  Sheau-  layer  forcing:  (a)  5/  =  2.0,  (b)  5/  =  3.3.  (fli  =  8,  Re  = 
1.5  X  10^). 

was  straightforward,  since  5/  was  known  o  priori  and  the  ratios  of  the  vortex  spac¬ 
ing  to  the  cylinder  diameter,  fi/2a  and  fj,/2a,  could  be  obtained  from  photographs 
like  those  in  Figs.  3.5  to  3.7.  See  Fig.  3.9.  This  analysis  assumes  that  the  dye  marks 
the  centers  of  the  vortical  structures.  The  normalized  celerity,  A,  ,  is  also  estimated 
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Fig.  3.8  Natural  (unforced)  shedding.  Flow  visuzJized  with  dye  issuing  as  in 
Fig.  3.2  -  Fig.  3.8. 


in  Sec.  4.3,  p.  4.18,  for  a  particular  forced  case,  using  LDV  measurements.  While 
A*  is  the  normalized  celerity  of  the  vortical  structures,  averaged  over  one  forcing 
cycle,  it  is  also  the  spacing,  normalized  by  the  distance  moved  by  the  flow  in  one 
forcing  cycle.  The  normalized  transverse  spacing  is  analogously  defined  as 


Fig.  3.9  Diagram  of  (fr)o,  (^s,)o,and 
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Note  that  a  negative  value  of  is  unphysical  sis  it  indicates  vortical  struc¬ 
tures  moving  ahead  of  the  cylinder.  Similarly,  a  value  of  A^  >  1  is  also  unlikely 
for  a  wake,  since  it  corresponds  to  vortical  structures  moving  downstream  faster 
than  the  free-stream  velocity  ahead  of  the  cyhnder.  In  a  finite  test  section,  how¬ 
ever,  blockage  effects  result  in  a  local  free-stresim  velocity,  Uq  ,  mestsured  above  smd 
below  the  wake,  that  is  somewhat  higher  than  Uoo  ■  Hence,  it  is  not  impossible  for 
A I  to  be  somewhat  greater  than  1. 

Figures  3.10  and  3.11  shows  the  variation  of  Aj  and  Ay  with  5/,  just  behind 
the  cylinder.  The  spacing  for  the  first  two  forcing  periods,  nearest  the  cylinder,  are 
shown.  The  squares  in  these  figures  denote  the  spacing  for  the  first  pair  of  vortical 
structures,  directly  behind  the  cylinder.  The  circles  denote  the  second  pair.  The 
vertical  lines  maxk  the  range.  The  dashed  vertical  lines  indicate  where  there  was  a 
rapid  contraction  of  the  spacing,  because  the  vortical  structures  were  being  “rolled” 
into  larger  lower  spatial  frequency  structures.  This  rapid  contraction  is  not  evident 
at  the  highest  values  of  5/ ,  because  the  roll  up  does  not  occur  until  several  more 
forcing  periods  downstream  of  the  cylinder. 

Note  the  decrease  in  A^  around  S/  =  1.  This  change  is  reminiscent  of  Trit- 
ton’s  (1959)  high-speed  and  low-speed  vortex  shedding  modes.  In  the  global  mode 
the  vortex  structures  move  downstream  quickly  (c/.  Tritton’s  high-speed  mode)  and 
in  the  local  mode  the  structures  move  downstreeim  more  slowly  (c/.  Tritton’s  low- 
speed  mode).  It  should  be  noted,  however,  that  because  of  the  finite  test  section 
height,  a  larger  wake  displacement  thickness,  especially  at  the  lower  values  of  5/, 
can  exaggerate  the  speed  of  the  vortical  structures  relative  to  the  cylinder.  It  is 
tempting,  anyway,  to  associate  the  instability  in  the  wake  structure  observed  above 
5/  ~  1  with  the  decreased  initial  speed  of  the  vortices. 
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Fig.  3.10  Variation  of  initial  A,  with  5/  (Qi  w  8).  Squares:  A^  for  the  first  pair 
of  structures,  circles,  for  second  pair. 


Figure  3.12  shows  the  variation  of  the  spacing  ratio. 


with  Sf.  The  squares  in  Fig.  3.12  denote  the  initial  spacing  ratio,  (Ay/Ai)^,  the 
circles,  the  spacing  ratio  one  period  farther  downstream,  (Ay/Ai)j.  Weihs  (1972) 
suggested  that,  to  first  order  behavior,  the  weike  would  become  narrower,  or  wider, 
depending  on  the  initial  spacing  ratio.  For  the  data  in  Fig.  3.12,  the  “neutral” 
value  of  (Ay/Ai)o ,  i.e.,  when  (Ay/Ax)^  =  (Aj,/Ar)j ,  which  marks  the  switch  from  a 
“widening”  wadce,  {^y/^z)o  <  (Ay/Ai)j ,  to  “narrowing”  one,  (Aj/Ar)^  >  (Ay/Ai)j , 
occurs  around  {Xy/Xz)^  =  (Ay/A^),  =  0.41  (5/  a  0.6).  Perhaps  related  to  the 
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Fig.  3.11  Variation  of  initial  Ay  with  5/  (fii  w8).  Squares:  A^  for  the  first  pair 
of  structures,  circles,  for  second  pair. 

switch  in  (Ay/Ax)g  at  5/  a  0.6  is  the  coincidence  of  the  value  of  the  wake  displace¬ 
ment  thickness  for  the  forced  case  and  the  unforced  case.  See  Fig.  3.17,  at  =  8. 
Around  (Ay/Ai)g  =  0.85  (5/  ss  1.2),  the  data  indicate  a  switch  back  from  “narrow¬ 
ing”  to  “widening.”  This  switch  marks  the  center  of  the  local  mode,  the  division 
between  the  global  and  the  shear  layer  modes.  Note  that  the  terms  “widening”  and 
“narrowing”  refer  to  the  ratio  of  the  transverse  and  streaimwise  spacing  of  the  vorti¬ 
cal  structures,  not  the  absolute  spacing.  The  spacing  of  the  vortical  structures  was 
also  measured  at  a  fixed  x/a  =  9,  downstream  of  the  cylinder.  Figure  3.13  shows 
the  v-ariation  of  the  spacing  of  the  structures,  with  5/,  for  =8,  and  Fig.  3.14, 
for  fli  =  4.  The  spacing  of  the  forced  vortical  structures  is  not  shown  in  the  figure, 
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Fig.  3.12  Vziriation  of  initial  ratio,  Ay/A,  with  5/  (fli  w8).  Squsu-es:  Ay/A^  for 
first  pair  of  structures,  auid,  circles,  for  the  second  pair. 

for  the  higher  values  of  5/ ,  because  they  were  not  discernible  at  i/a  ~  9 . 

Note  that  within  the  spread  of  the  data.  A,,  Ay,  and  Ay/A^,  at  xja  ~  9, 
seem  to  be  constant  over  most  of  the  range  shown.  They  do  seem  to  be  a  function  of 
fl ,  though.  This  suggests  that,  for  a  given  Qi ,  the  geometry  of  the  wake  structure  is 
similar  over  that  range.  Note  that  the  vertical  spacing  Ay  is  depressed  at  the  lower 
values  of  5/.  This  is  probably  because  of  the  finite  test  section  height.  In  addition, 
the  scatter  in  the  data,  especially  in  the  fli  =  4  case,  may  be  attributable  to  the 
subjective  method  of  locating  the  centers  of  the  vortical  structures.  For  =  4, 
the  average  values  of  A^ ,  Ay,  and  Ay/Ai  were  0.96,  0.28,  and  0.3  respectively.  For 
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Fig.  3.13  Variation  of  A^,  Ay,  zmd  Ay/A,  with  5/  at  xfa  =  9  with  fli 
Triangles:  Af ,  squares:  Ay ,  and  diamonds  Ay/ A^ . 


fli  =  8,  they  were,  0.89,  0.37,  and  0.41.  Note  that,  at  the  same  value  of  flj  = 
average  value  of  Ay /A*  at  i/a  =  9  matches  the  “neutral”  value  of  (Ay/Ai)^ 
found  earlier. 
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Fig.  3.14  Variation  of  A*,  A^,  and  Ay/A,  with  5/  at  xja  =  9  with  fii  =  4. 
Triangles:  A,  .squares:  Ay ,  and  diamonds  Ay/ A,. 

3.3.3  A  rough  estimate  of  the  vortex  strength 

A  measure  of  the  strength  of  the  vortical  structures  in  the  wake  can  be 
obtained  from  their  speed  and  spacing.  A  vortex  street  will  have  a  self-induced 
velocity  tx, ,  which  will  depend  on  the  strength,  Fv  ,  and  the  streamwise  and  lateral 
spacing  of  vortical  structures,  and  £y ,  and  also  on  factors  such  as  the  presence 
of  the  cylinder,  finite  test  section  height,  finite  core  size  of  the  vortical  structures, 
and  viscous  effects.  See  Goldstein  (1938,  §242-246)  and  Roshko  (1954c)  for  similar 
estimates  of  the  strength  of  vortical  structures  in  a  bluff  body  wadce. 
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Using  the  potential  flow  model  of  a  periodic  vortex  street,*  of  the  strength 
of  the  vortical  structures  in  the  wake,  it  can  be  seen  {e.g.,  Goldstein  1938)  that  the 
strength  per  unit  strecimwise  spacing  of  the  vortical  structures  in  the  wake  can  be 


written  as 


(3.3.3) 


This  expression  is  not  valid  at  £y  =  0,  because  when  the  structures  are  lined  up, 
u,  =  0  and  Fy  can  be  any  value.  In  addition,  u,  can  be  related  to  the  velocity  of 
the  'wtical  structures,  Uy,  from  Sec.  3.3.2,  i.e.. 


Uj  —  Uo  Uy  ^  Uo  ^  £x  * 


(3.3.4) 


Uo  varies  with  downstream  position,  and  is  not  known  a  priori.  Using  continuity, 
however,  and  neglecting  the  core  size  of  the  vortical  structures,  Uo  can  be  written 
in  terms  of,  Fy/^*,  Uooi  and  £y/h,  i.e., 

%  =  I'oo  +  ^  ^  ,  (3.3.5) 

where  h  is  the  test  section  height.  A  little  algebra,  and  the  definitions  for  and 


A, ,  then  yield 


_  2  _ ^  ~ 

Uoo^z  tanhf^^^  - 


(3.3.6) 


^0°^*  tanh  ^ 

Notf  that  the  singularity  in  this  expression  is  directly  related  to  the  fact  that  the 
spe<.  d  of  the  vortex  street  is  zero  for  the  czise  £y  =  0 . 


Figure  3.15  shows  the  variation  of  TvfUoo^x  with  Sj.  The  dashed  line  in 
the  figure  is  a  plot  of  Eq.  3.3.6  assuming  constant  A^  and  Aj, .  The  rise  in  Fy/Uoo^r 
at  the  lowest  values  of  S/  could  be  an  artifact  of  finite  test  section  height  effects, 
represented  by  the  the  iy/h  term  in  Eq.  3.3.6.  Similar  plots  for  the  Qi  =  4  case 
were  not  interpretable,  because  the  spread  in  the  data  was  too  large. 

•  Using  such  a  model  can  only  yield  a  relative  measure  of  the  strength  of  the  vortices  in  the  wake. 

See  Sec.  4.3  for  a  better  measure  of  the  strength  of  the  vortical  structures,  for  a  particular  forced 


case. 
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Sf 

Fig.  3.15  V^fUooiz  vs.  Sj.  (fti  =8). 


Similar  to  Fy/l/oo^i?  a  measure  of  the  strength  of  an  individued  vortical 
structures  can  be  written  as 


Tv 

Uooa 


_  . 

\  ^/  /  tanh  ^ 


(3.3.7) 


Figure  3.16  shows  the  variation  of  Fv/t^ooU  with  Sf.  The  dashed  line  in  the  figure 
is  Eq.  3.3.7,  zissuming  constant  Xz  and  Ay. 


For  £y/h  1,  Eqs.  3.3.6  and  3.3.7  become 


Fv 

Uootz 


2  (  1  -  A,  )  coth  ' 


(3.3.8) 


$ 


Though  not  really  applicable  to  this  idealized  case,  it  was  suggested  by 
Roshko  (1954c)  that,  for  a  real  wake,  the  mean  local  velocity  along  the  line  of 
vortex  centers,  u*,  bt  sed  to  calculate  the  velocity  of  the  vortex  street  relative  to 
the  fluid,  where 

=  2  (3.3.10) 

for  the  potential  flow  model  discussed  here.  Using  u*  instead  of  Uj  in  the  preceding 
analysis  is  equivalent  to  setting  tanh(7r£y/£,)  =  1  in  the  Eqs.  3.3.6  to  3.3.9.  This 
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leads  to  simpler  estimates, 


and 


Tv 

Uooi. 


2  ( 1  -  A,  )  , 


Tv 

Uooa 


4{l-  . 

•>/ 


(3.3.11) 


(3.3.12) 


3.3.4  Ejection  of  circulation  into  the  flow 

The  rate  at  which  circulation  is  introduced  into  a  flow  about  a  body  can  be 
written  {e.g.,  Leonard  1987)  as 

^1  =  ~2>lbody^body  ,  (3.3.13) 

at  I flow 

where  Abody  is  the  cross-sectional  area  of  the  body  and  ^body  is  the  rotational 
acceleration  of  the  body.  In  words,  net  circulation  is  introduced  into  a  flow  from  a 
body  only  when  it  experiences  rotary  acceleration.** 


Integrating  Eq.  3.3.13  with  respect  to  time  yields  for  the  case  of  a  circular 
cylinder  of  radius  a , 

^Lfiow  =  —  27ra^A^body  ,  (3.3.14) 

where  A  denotes  “change  in.”  Hence,  when  a  rotating  cylinder  goes  from  one 
rotation  rate  to  another,  a  net  circulation  will  be  ejected  from  the  surface.  In 
dimensionless  form  this  is 

=  -27rAn  (3.3.15) 

UooO 


*•  It  is  important  to  note  that  “flow”  here  refers  to  both  the  boundary  layer  and  outer  flows. 
To  accentuate  this  point  recall  that  an  impulsively  started  and  lifting  airfoil  has  no  rotational 
acceleration  but  ejects  a  starting  vortex  into  the  outer  flow.  There  is  no  paradox  because  an 
amount  of  circulation  equal  to  that  shed  into  the  outer  flow  is  trapped  in  the  boundary  layer 
of  the  airfoil. 
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If  all  of  the  circulation  of  a  particular  sign  created  during  a  forcing  period  and 
due  to  the  angular  acceleration  of  the  cylinder  where  put  into  alternating  vortices 
Eq.  3.3  .15  would  indicate  a  vortex  strength  of  up  to 

=  AirQi  .  (3.3.16) 

Uoo^ 

For  the  present  data  fli  =4  and  8.  This  yields  values  of 

7^  «  50;  100  (3.3.17) 

t/oofl 

which  are  much  larger  tham  those  presented  in  Fig.  3.16.  That  the  strength  of 
the  vortical  structures  in  Fig.  3.16  is  only  rough  does  not  account  for  the  large 
discrepancy  between  the  values  in  Fig.  3.16  and  Eq.  3.3.17.  A  plausible  explanation 
for  this  discrepancy  is  that  most  of  the  vorticity  generated  at  the  cylinder  surface  is 
reabsorbed,  or  is  cancelled,  combining  with  vorticity  of  the  opposite  sign,  before  the 
vortical  structure  is  shed  into  the  wake.  This  is  consistent  with  the  observation  that 
the  vortical  structures  for'T  very  close  tc  the  cylinder  and  hence  have  time  to  absorb 
both  the  positive  and  the  negative  circulation  generated  by  the  cylinder’s  rotational 
acceleration.  See  Sec.  4.3  for  a  description  of  the  flow  close  to  the  cylinder. 


It  is  likely  that  the  discrepancy  will  decrease  at  the  lowest  values  of  5/, 
because  the  trend  is  for 


Fv 

Uooa 


(3.3.18) 


cf.  Eq.  3.3.7.  At  the  very  lowest  S/,  the  shedding  of  vortices  appear  much  like 
alternating  “starting  vortices.” 
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3.4  Displacement  thickness  and  drag  coefficient  estimation 


In  order  to  provide  a  more  quantitative  measure  of  the  effects  of  forcing 
on  the  cylinder  and  the  resulting  wake,  the  cylinder  drag  and  wake  displacement 
thickness  were  estimated  from  wake  streamwise  mean  and  rms  velocity  profiles.  For 
flow  that  is  two-dimensional,  the  displacement  thickness  is  normzdly  defined  as 

£1  = 

h  -  Jo  [  Uo 

In  this  expression,  rj  =  yfh  is  the  normadized  vertical  (cross-stream)  position,  u{r)) 
is  the  mean  streamwise  velocity,  y  and  h  are  the  vertical  position  and  water  channel 
depth  respectively,  and  Uq  is  the  velocity  in  the  (free-stream)  region  outside  the 
cylinder  wake  and  water  channel  boundary  layers.  Recall  Uq  and  Uoo  in  Fig.  1.1. 
By  analogy,  the  “displaced  area”  for  a  three-dimensional  flow,  can  be  defined  as 


dtj  . 


(3.4.1) 


(3.4.2) 


where  A  is  the  test  section  cross-sectional  sirea,  and  C  =  2/6  is  the  normalized 
spanwise  coordinate,  with  z  the  spanwise  coordinate  and  b  the  water  charnel  span. 
Continuity  then  yields  for  Uoo,  the  free-stream  velocity  far  ahead  of  the  cylinder. 


Uoo  =  J  «(t?,C)  drjdC  . 

Using  Eq.  3.4.3  to  simplify  3.4.2  then  yields 

^  =  1  - 
A  Uo  ■ 


(3.4.3) 


(3.4.4) 


This  expression  is  useful  in  that  it  provides  information  about  the  flow  over  an 
entire  cross-section  of  the  water  channel,  while  requiring  only  the  measurement  of 
Uo  and  Uoo  (see  Fig.  1.1). 
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In  order  to  remove  the  contribution  of  the  initial  boundary  layer  (in  the 
absence  of  the  cylinder),  to  leading  order,  it  is  useful  to  compare  Uq  (with  the 
cylinder  in  place)  with  the  undisturbed  flow  velocity  in  the  absence  of  the  cylinder, 
Uq  ,  at  the  same  streamwise  location.  This  yields  an  expression  for  the  approximate 
displaced  area  by  the  wake;  i.e.. 


A  Uo  ■ 


(3.4.5) 


Figure  3.17  depicts  the  variation  of  A”  f  A  with  5/  and  fli .  The  dashed  curves  in 
this  figure  denote  A* /A  levels  equal  to,  or  less  than,  those  corresponding  to  the 
unforced  case.  The  parameter  space  denoted  by  the  solid  lines  represents  an  increase 
in  A*  J A  over  the  unforced  case.  The  displacement  area  can  be  made  substantially 
larger,  or  smaller,  depending  on  the  frequency  ai  .-  amplitude  of  oscillation.  The 
minimum  for  A* /A  over  the  range  shown  occurs  in  the  neighborhood  of  5/  a  1, 
and  fli  ss  3,  and  is  roughly  half  that  of  the  unforced  case.  Measurements  for  A‘/A 
were  taken  in  the  range  0.17  <  5/  <  3.3  in  steps  of  about  0.17,  and  0  <  flj  < 
16  in  steps  of  1.  The  data  presented  in  Figs.  3.17  and  3.18  W3is  averaged  over 
2  minutes  (approximately  40  natural  unforced  shedding  cycles).  Wake  midspan 
velocity  profiles  were  measured  for  5/  in  the  range  of  0.15  to  1.4,  holding  fixed 
at  f2i  =  2.  Figure  3.18  compau-es  the  wake  profiles  for  several  values  of  5/,  at 
fli  =  2. 


For  a  finite  test  section  and  flow  that  is  two-dimensional  in  the  mean  over  a 
portion  of  the  the  span,  it  can  be  shown  that  the  sectional  drag  coefficient  can  be 
estimated  by  the  expression  (Dimotakis  1978): 


^  (1  -  6-/hr  [Jo  Uo  \  Uo  )  Jo  Ul  2  i  h  j  • 

(3.4.6) 

In  this  expression,  u'  and  v'  are  the  streamwise  and  cross-stream  rms  velocities,  and 
6*  is  the  displacement  thickness  defined  in  Eq.  3.4.1.  The  displacement  thickness 
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Fig.  3.17  Variation  of  A* /A  with  Sf  and  fii .  Dashed  lines  mark  the  region  where 
the  wake  displacement  thickness  is  less  than  in  the  unforced  case  (note 
the  minimum  around  5/  =  1  and  =3).  Solid  lines  mark  the  region 
where  it  is  greater  than  or  equal  to  the  unforced  case. 


emd  drag  coefficient  estimates  summzurized  in  Figs.  3.19  and  3.20  were  calculated 
from  midspan  mean  velocity  profiles,  using  Eqs.  3.4.1  and  3.4.6,  neglecting  the 
fluctuating  terms,  which  are  not  only  small,  but  very  nearly  cancel  (Dimotakis 
1978).  As  can  be  seen,  there  is  a  broad  minimum  in  Cd  and  S*/h  around  S/  =  1 . 
The  estimated  Cq  for  the  unforced  case  is  a  factor  of  six  greater  than  for  this  forced 
case.  Similarly,  there  is  a  factor  of  five  reduction  in  8*/h.  Note  that  the  value  of 
8*  fh  for  a  particular  Sj  and  flj  is  generally  less  than  the  corresponding  value 
of  A* /A,  indicating  an  increase  in  8* /h  outside  the  midspan  region.  Since  Cd 
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Fig.  3.18  Cylinder  wake  mean  velocity  profiles,  measured  at  midspan  (fli  =2). 
Squares:  unforced,  circles:  5/  =5  0.2,  emd  triangles:  S/  ~  1. 

increases  with  8*  [h  it  should  also  increase  outside  this  midspan  region.  A  spanwise 
variation  in  8*  Jh  may  be  attributed  to  the  interaction  of  the  sidewalls  with  the  shed 
vortices.  See  Roberts  (1985),  Roberts  &  Roshko  (1985),  Koochesfahani  (1987),  and 
Kurosaka  et  al.  (1988). 

Figure  3.21  shows  the  correlation  between  the  calculated  sectional  drag  coef¬ 
ficients  and  sectional  wake  displacement  thicknesses.  The  dashed  line  in  that  figure 
denotes  the  approximate  relationship  (e.g.,  Batchelor  1970,  §5.12), 

s* 

Cd  ~  —  1  (3.4.7) 

a 


6Vh 
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Fig.  3.19  Cylinder  wake  displacement  thickness  variation  with  5/  (fli  =2). 
for  wake  profiles  taken  far  downstream  of  the  cylinder. 

Another  loczJ  minimum  exists  near  5/  w  0.2.  The  flow  visualization  in 
Fig.  3.4  suggests  that  this  corresponds  to  the  duu  mode  in  which  the  cylinder 
releases  two  vortices  of  opposite  sign  during  each  half-cycle.  When  one  of  these  two 
vortices  shed  from  the  bottom  appears  on  top,  a  dip  is  observed  in  the  corresponding 
velocity  profile.  See  circles  in  Fig.  3.18.  This  phenomenon  was  only  observed  in  a 
small  region  of  the  ( S/,  Jl]  )-parameter  space. 
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Fig.  3.20  Drag  coefficient  variation  with  5/  (fii  =2). 

3.5  A  note  on  the  interpretation  of  linear  stability  analyses 

This  section  was  originally  written  in  response  to  a  reviewer  of  the  paper 
“Rotary  oscillation  control  of  a  cylinder  wake,”  by  Tokumaru  Dimotakis  (1991). 
It  is  included  here  as  a  point  of  interest. 

There  has  been  much  interest  in  the  use  of  linezu"  stability  analyses  to  describe 
unforced  bluff  body  flows.  See  Koch  (1985),  Triantafyllou  et  al.  (1986),  Monkewitz 
Nguyen  (1987),  Provansal  et  al.  (1987),  Chomaz  et  al.  (1988),  Monkewitz  (1988), 
and  Karniadakis  &  Triantafyllou  (1989).  Using  linear  theory,  it  has  been  noted  that 
the  response  to  periodic  forcing  of  the  absolutely  unstable  near  wake  region  of  a 
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Fig.  3.20  Drag  coefficient  variation  with  5/  (fii  =2). 
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5*/h 


Fig.  3.21  Co  vs.  S*/h.  Circles:  forced  (fti  =2),  square:  unforced  (fti  =0). 

bluff  body  will  be  overwhelmed  by  the  “intrinsic”  vortex  shedding  mode  observed 
in  the  absence  of  external  forcing.  See  Huerre  k  Monkewitz  (1985),  Monkewitz  k 
Nguyen  (1987),  Chomaz  et  al.  (1988),  aind  Karniadeikis  k  Triantafyllou  (1989).  This 
argument  is  based  on  a  linearized  model  and  suggests  that  the  periodic  response 
will  be  overwhelmed  because  the  treinsient  response  grows  exponentially  in  time 
while  the  forcing  is  only  periodic  in  time.  Since  the  system  cannot  be  considered 
linear  for  large  excursions  from  the  “base”  state,  however,  this  argument  can  only  be 
relied  upon  for  infinitesimal  forcing  amplitudes  and  short  times.  It  is  not  surprising, 
perhaps,  that  this  interpretation  of  the  results  of  linear  theory  is  not  in  accord  with 
experiments  and  simulations  which  show  that  finite  amplitude  forcing  of  the  wake 
of  a  bluff  body  can  result  in  a  wake  structure  that  is  synchronous  with  the  forcing. 
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over  a  range  of  forcing  frequencies  and  amplitudes;  e.g.,  Koopman  (1967),  Provansal 
et  al.  (1987),  and  KarniadzJds  &  Triantafyllou  (1989). 

Provansal  et  al.  (1987),  examined  the  (nonlinear)  Landau  stability  model 
(Landau  1944,  Stuart  1958,  1960)  and  showed  how  it  could  be  used  to  describe 
qualitatively  the  conditions  necessary  for  synchronization  at  small  but  finite  forcing 
amplitudes,  at  Reynolds  numbers  in  the  vicinity  of  the  onset  of  vortex  shedding. 
See  also  Landau  &  Lifshitz  (1987,  §26  &  §30).  Evidently,  while  the  linear  stability 
analysis,  and  hence  the  notion  of  the  bluff  body  wake  as  a  globally  vmstable  flow  with 
a  region  of  absolute  instability,  is  useful  in  describing  the  system  dynamics,  given  an 
a  priori  knowledge  of  the  final  averaged  flow  (Karniadakis  A:  Triantafyllou,  1989), 
prediction  of  the  final  dynamic  state,  especially  in  the  presence  of  finite  external 
forcing,  requires  a  different  approach. 


CHAPTER  4 


A  particular  forced  case 


In  this  chapter,  a  particular  case  of  a  cylinder  executing  forced  rotairy  oscilla¬ 
tions  in  a  uniform  stream  is  examined.  An  estimate  of  the  unsteady,  phase-ave.  ged 
streaunfunction  was  made.  The  speed  and  strength  of  the  shed  vortical  structures 
was  also  estimated. 


4.1  Experimental  setup 

The  experiments  documented  here  were  performed  in  the  20"  x  20"  Free 
Surface  Water  Tunnel  at  GALCIT.  The  cylinder  used  for  the  near  wake  velocity 
meaisurements  is  an  amodized  amd  machined  aluminum  tube,  4"  in  diameter,  and 
mounted  and  driven  by  the  apparatus  described  in  Ch.  2.  The  command  signail  for 
the  cylinder  velocity  was  generated  by  a  computer  controlled  function  generator. 
The  position  of  the  cylinder  output  by  the  BEI  13-bit  absolute  position  encoder  was 
read  using  the  computer’s  parallel  interface. 

Velocities  were  meaisured  using  a  one-chamnel,  Bragg  cell  frequency- shifted 
laiser  Doppler  velocimeter  with  a  trainslatable  (and  rotatable)  focad  volume.  A  Lexel 
Model  85  Argon-ion  laser  operating  in  single  line  mode  ( 514.5  nm)  at  about  200  mW 
was  used  for  the  LDV.  The  laiser  beam  for  the  LDV  velocity  measurements  was  split 
into  two  beams  of  nearly  the  same  intensity  using  a  cube  beam  splitter.  The  beams 
were  then  Bragg  cell  frequency-shifted  to  allow  for  flow  reversals.  The  Bragg  cell 
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offset  between  the  two  beams  weis  200  kHz  for  the  data  presented  in  this  chapter.  A 
Dove  prism  was  introduced  in  line  after  the  Bragg  cells  for  the  purpose  of  rotating 
the  orientation  of  the  measurement  volume.  An  achromat  was  used  to  focus  the 
beams  in  a  focal  (measurement)  volume  located  in  the  midspan  plane  of  the  water 
tunnel.  The  light  scattered  from  the  focal  volume  was  gathered  by  an  achromat 
and  focused  on  a  photodiode  with  an  integral  low  noise  pre-amplifier  designed  by 
Dan  Lang.  The  signal  from  the  pre-amplifier  was  then  band-pass  filtered  above  and 
below  the  frequency  corresponding  to  the  mean  flow  velocity.  The  band-pass  range 
was  chosen  to  pass  instantaneous  velocity  fluctuations.  A  tracking  phase-locked 
loop  designed  by  Dan  Lang  and  Paul  Dimotakis  was  then  used  to  lock  a  TTL 
square  wave  to  the  dominant  frequency.  The  TTL  signal  was  read  by  a  counter¬ 
timer  board  on  a  data  acquisition  computer  and  stored  on  disk  for  later  processing. 
The  transmitting  and  receiving  hardware  was  mounted  on  a  motor  driven  x-y-z 
traverse  (positioner). 


4.2  Measurement  details 

For  adl  the  data  presented  in  this  chapter,  5/  =  0.94,  =  2 ,  and  Re  =  2.1  x 

10^ .  This  Reynolds  number  was  chosen  to  be  higher  than  in  Ch.  3  to  shorten  the 
time  scale  of  the  flow  dynamics,  and  hence  also  of  the  experiment.  The  tangential 
velocity,  V0{r,  6)  was  measured  along  eleven  B  =  constant  rays  in  the  upper  half 
plane,  at  15  degree  intervals.  Along  most  of  these  rays,  ve  w«is  measured  in  steps 
6r/a  =  0.01 ,  beginning  near  the  cylinder  surface,  for  1.01  <  t  fa  <  1.15,  where  r/a 
is  the  normalized  radial  position.  Then  in  steps  of  6r/a  =  0.025  for  1.15  <  r/a  < 
1.25,  steps  of  8r/a  =  0.05  for  1.25  <  rfa  <  1.5  and  finally  in  steps  of  Sr/a  =  0.25 
for  1.5  <  r/a  <  2.5.  Figure  4.1  shows  the  B  =  constant  measurement  rays.  The 
velocities  at  each  location  were  meaisured  at  f  6t  =  0.02  time  intervals.  At  each 
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location,  41  forcing  periods  were  averaged  at  constant  forcing  phase,  yielding  50 
realizations  of  the  flow  averaged  at  constant  forcing  phase.  The  velocity  field  in  the 
lower  half  plane  was  taken  to  be  the  same  as  that  at  the  top,  but  a  half  forcing 
period  out  of  phase. 

Before  the  velocity  data  were  averaged  at  constant  phase,  they  were 
smoothed  in  time  using  a  gaussian  filter  with  a  3  dB  cutoff  around  five  times  the 
forcing  frequency.  No  other  explicit  smoothing  was  performed,  but,  because  of  the 
wide  spacing  of  the  6  =  constant  rays,  there  was  sm  implicit  smoothing  of  data 
interpolated  in  between. 


Fig.  4.1  Measurement  locations. 
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4.2.1  Streamfunction 

Calculation  of  the  stream  function  was  straightforward.  Given  that  the 
stream  function  'J'  is  known  at  some  point  in  the  flow  xo ,  and  the  velocity  normal 
to  any  path  C  connecting  the  two  points  x  and  Xo ,  ««  =  u  •  n,  is  also  known,  the 
(dimensionless)  stream  fimction  V’(x)  can  be  calculated,  for  two-dimensional  flows, 
by  the  expression 

V.  (x/<.)  =  ^=  /  + 

U oo^  Jc  O, 

Further,  since  the  stream  function  is  constant  along  a  solid  boundary  surface,  it  is 
convenient  to  begin  integrating  from  there,  and  to  choose  some  constsint  V’(xo/a)  = 
00  on  that  surface.  For  the  present  experiments,  Eq.  4.2.1  was  written  as 

=  r''  */»  ■  ('‘•2.2) 

J\  ^oo 


It  should  be  noted  that  for  a  number  of  the  measurements,  the  boundary 
layer  neeu’  the  wzdl  could  not  be  resolved  with  the  present  setup,  because  of  large 
gradients  in  the  velocity,  and  difficulties  in  positioning  the  measurement  volume. 
This  did  not  cause  large  errors  in  the  calculation  of  0,  however,  because  the  velocity 
at  the  wall  was  known  o  priori.  For  the  integration  of  Eq.  4.2.2,  a  (straight)  line 
was  fit  between  the  velocity  at  the  wall,  and  the  nearest  measured  velocity. 

The  forced  shedding  process  for  5/  w  1  and  fli  ss  2  is  illustrated  in  Figs. 
4.3  to  4.12.  The  illustrations  on  the  top  denote  fluid  marked  with  dye  introduced 
just  ahead  of  the  ten  o’clock  position  on  the  cylinder.  They  were  traced  from  a 
video  tape  of  the  flow.  The  circles  with  the  Q  in  them  indicate  the  direction  and 
magnitude  (arbitrary  scale)  of  the  rotation  rate.  The  illustrations  on  the  bottom 
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show  the  stresLin  lines  calculated  from  data  averaged  at  constant  forcing  phase.  Brief 
notes  on  the  averaged  constant-phase  streamline  data  are  contained  in  the  captions 
of  the  Figs.  4.2  -  4.12*.  Figure  4.2  illustrates  the  average,  over  all  forcing  phases, 
of  the  stream  function. 


The  stream  function  data  were  first  estimated  by  integrating  along  6  = 
constant  rays.  Between  the  measurement  points,  the  streamfunction  was  estimated 
using  spline  interpolation,  first  along  the  rays,  and  then  between  them.  The  stream¬ 
line  plots  were  generated  using  a  contour  plotting  routine  with  points  interpolated 
on  a  300  x  400  rectangular  grid. 
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Fig.  4.2  Streamlines  averaged  over  all  phases.  Note  the  delayed  average  separation 
point  and  the  small  averaged  “recirculation”  region  realized  in  this  forcing 
regime.  This  looks  very  much  like  the  (unforced)  time  averaged  flow 
around  the  cylinder  at  higher  Reynolds  numbers  ( >  5  x  10^ ),  which  also 
has  a  similar  drag  coefficient  of  about  0.2  to  0.3. 


In  Figs.  4.3  -  4.12,  the  flow  visualization  data  are  illustrated  on  the  top  and 
the  phase  averaged  stream  line  pattern  calculated  from  Idv  data  are  on  the  bottom. 


The  entire  set  of  50  phase  averaged  realizations  of  the  flow  is  included  in  .\ppendix  B. 
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The  dot  in  the  figures  mark  the  approximate  cylinder  position,  the  airrows,  the 
direction  of  rotation,  and  the  shaded  regions,  the  range  of  motion. 

It  appears  that  the  forced  shedding  mechanism  outlined  here  replaces  the 
natural  shedding  mechanism  with  one  that  gathers  vorticity  close  to  the  cylinder 
surface,  and  then  releases  a  vortical  structure  into  the  flow.  This  would  suggest 
that  the  formation  of  vortical  structures  in  the  forced  case  is  less  affected  by  the 
stability  characteristics  of  the  flow  downstream  of  the  cylinder.  This  is  in  contrast 
to  observations  of  unforced  vortex  shedding  which  show  that  vorticity  is  introduced 
into  the  flow  from  more-or-less  steady  sepziration  points,  with  the  resulting  vortical 
structures  forming  in  the  wake  of  the  cylinder. 

When  comparing  the  streamline  patterns  in  this  chapter  and  Appendix  B 
with  simulations,  several  factors  should  be  considered.  As  can  be  seen  from  Fig.  4.1, 
the  spatial  resolution  of  the  data  for  r/a  >  1.5  corresponds  to  only  a  handful  of 
measurement  locations  per  vortical  structure  in  that  region.  The  continuity  of  the 
velocity  field,  however,  diminishes  this  problem.  As  stated  in  Sec.  4.2,  the  data  in 
the  lower  half  plane  was  taken  to  be  the  szune  as  in  the  upper  half  plane,  but  a 
hadf  forcing  cycle  out  of  phase.  Because  of  asymmetries  in  the  flow,  this  match  (at 
y  =  0 )  is  not  perfect  and  causes  small  “glitches”  in  the  streamline  pattern  around 
the  joint  between  the  upper  amd  lower  half  planes.  Compauing  streamline  patterns 
a  half  period  apau't  gives  an  indication  of  the  asymmetry  in  the  data.  Because 
the  present  data  wais  measured  at  midspam,  there  is  no  out-of-plane  velocity,  in 
the  mean.  This  does  not,  however,  remove  three-dimensional  effects  like  vortex 
stretching  and  instamtaneous  out-of-plane  flow. 

Because  of  factors  such  as  these,  comparisons  of  the  present  data  with  sim¬ 
ulations  will  not  be  perfect.  It  should  be  emphasized  that  the  basic  mechamisms 


-4.7- 

gathering  vorticity  close  to  the  cylinder  (allowing  it  to  mix  and  cancel)  and  then 
ejecting  it  into  the  flow,  should  not  be  strongly  affected  by  these  factors.  A  simu¬ 
lation  that  captures  the  behavior  of  the  unsteady  boundary  layer  in  this  flow  will 
likely  have  results  similar  to  the  present  data. 


Fig.  4.3  ft  =  0.0.  At  this  point,  the  cylinder  has  just  finished  shedding  a  vortex 
(at  i/a  =  1.5),  and  the  kernel  (bump)  of  the  next  vortex  to  be  shed  can 
be  seen  at  about  45  degrees  on  the  cylinder  surface. 
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Fig.  4.4  ft  =  0.1.  The  bump  has  grown  larger,  note  also  the  beginnings  of  a 
second  bump  on  the  cylinder  surface  at  about  90  degrees. 
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Fig.  4.6  ft  =  0.3.  The  bumps  continue  to  grow.  Note  that  the  main  bump  seems 
to  be  held  in  place  relative  to  the  cylinder  surface,  possibly  by  the  flow 
induced  by  the  previously  shed  vortices.  Note  that  the  cylinder  motion 
is  in  the  clockwise  direction  which,  intuitively,  would  tend  to  have  the 
opposite  effect. 
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Fig.  4.7  ft  =  0.4.  Again,  the  bumps  continue  to  grow,  but  now  are  becoming 
more  rounded.  The  large  bump  almost  looks  as  if  it  is  being  “pinched” 
into  the  flow. 
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Fig.  4.8  ft  =  0.5.  The  presence  of  the  secondary  bump  is  less  noticeable;  most 
likely  it  has  been  merged  with  the  primary  structure. 
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Fig.  4.9  ft  =  0.6.  The  cylinder  begins  to  move  in  the  opposite  direction  to  the 
flow  on  that  side.  This  appears  to  be  when  the  vortex  begins  to  be 
released  into  the  flow. 


Fig.  4.11  ft  =  0.8.  The  main  vortex  begins  to  move  off  to  the  right.  The  next 
vortex  to  be  shed  continues  to  be  held  in  place  both  by  the  vortices 
which  have  been  shed,  and  the  rotation  of  the  cylinder. 
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4.3  Strength  of  shed  vortices 

Section  3.3.3  proposed  a  measure  of  the  strength  of  the  vortical  structures 
shed  into  the  wake.  In  this  section,  another  estimate  for  the  circulation  of  the  vor¬ 
tical  structures  is  made  by  measuring  the  circulation  of  the  flow  about  the  cylinder, 


rc(r)  =  r  f  vg(r,0)dO  , 
Jc(r) 


for  V8urying  phase,  /t,  and  radius  of  integration,  r.  The  associated  dimensionless 
quantity  is 

rc(r) 

Uooa  ‘ 

Figure  4.15  shows  the  variation  of  Tc/UooO.  with  ft  and  r/a. 


The  circulation  in  an  annulus  between  two  circles  of  radii  rj  and  r2  can  be 


written  as 

^  ^  .  (4.3.1) 

t/ooflln  Uooa  UooO. 

Hence,  an  estimate  for  the  strength  of  a  vorticaJ  structure  (residing  in  an  annulus) 


_  rc(rn,4x)  I  _  rc(rniin)  (4  3  2) 

UooO’  UooO  ImM  UooO  min 

where  r^ax  and  Tmin  are  chosen  to  maximize  Eq.  4.3.2  in  the  wake.  The  data 


in  Fig.  4.13  is  typical  of  the  variation  of  TdUooO  with  radius  of  integration  r/a 
at  a  constant  phase.  Figure  4.14  is  the  strezunline  pattern  related  to  the  data  in 
Fig.  4.13.  From  data  like  Fig.  4.13,  the  estimated  strength  of  the  vortical  structures 


(4.3.3) 


In  addition,  from  the  slope  of  the  contours  in  Fig.  4.15, 

d(r/a)  _  Ax 

d{ft)  Sf  ' 
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Fig.  4.13  Variation  of  the  normalized  circulation,  Fc/I/ooO,  taken  round  contours 
of  varying  radius,  rja  {ft  =  0.54).  The  rise  in  the  curve  between 
rja  —  1.2  and  1.7  marks  a  patch  of  positive  (clockwise)  circulation 
and  the  drop  between  r/a  ~  1.7  and  2.3  a  negative  (counter-clockwise) 
patch,  cf.  Fig.  4.14. 

it  is  found  that  the  average  (normalized)  convection  speed  of  the  vortices  in  the 
near  wake,  based  on  the  slope  of  the  T^/UooO^  =  0  contour,  is 

A,  5S5  0.6  ,  for  -  >  1.2  .  (4.3.4) 

a 

The  transverse  spacing  of  the  vortices  was  difficult  to  estimate  but  appezu-ed  to  be 
in  the  range 


0.2  <  Ay  <  0.3  , 


(4.3.5) 
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Fig.  4.14  Streamline  pattern  associated  with  the  data  in  Fig.  4.13  {ft  =  0.54). 

Interestingly,  using  the  method  of  Sec.  3.3.3,  p.  3.17,  for  the  periodic  vortex  street 
yielded  estimates  close  to  Eq.  4.3.3.  Using  Eq.  3.3.7,  amd  and  Ay  from  Eqs. 
4.3.4  and  4.3.5,  yielded  an  estimate  for  the  strength  of  the  shed  vortical  structures, 

1.1  <  <  1.4,  (4.3.6) 

C  ocU 


and  similarly  for  the  simplified  estimate.  Eq.  3.3.12,  which  takes  no  account  of  the 
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transverse  spacing,  a  value 

K  1.0.  (4.3.7) 

U<x>o. 

What  is  remarkable  is  that  very  little  of  the  circulation,  calculated  at  the 
cylinder  surface,  actually  finds  its  way  into  the  flow.  If  all  of  the  circulation  of  a 
particular  sign  created  during  a  forcing  period  and  due  to  the  angular  acceleration 
of  the  cylinder  where  put  into  alternating  vortices  Eq.  3.3.16  would  indicate  a  vortex 
strength  of  up  to 

=  4:rn,  a  25  (4.3.8) 

Uooa 

for  the  data  in  this  chapter.  For  the  forced  case  presented  in  this  chapter,  the 
strength  of  the  vortical  structures  shed  into  the  wake  is  only  about  4%  of  this 
value. 


Fig.  4.15  Fc/t^oofl  va.  r/a  and  ft.  The  daished  lines  denote  net  counter-clockwise 
rotation;  the  solid  lines,  clockwise.  The  value  TdUooa  at  contour  level 
‘1’  =  0.6,  ‘2’  =  0.5,  ‘3’  =  0.4  ... ,  down  to,  ‘13’  =  -0.6.  Contour  level 
‘7’,  solid  line,  marks  the  zero  contour. 


CHAPTER  5 


Lifting  cylinder 


5.1  Preface 

In  this  chapter,  a  novel  method  for  estimating  the  lift  on  the  cylinder  is  pre¬ 
sented.  The  data  calculated  using  this  method  compare  favorably  with  similzu-  data 
published  for  steadily  rotating  cylinders.  In  addition,  the  effect  of  superimposing 
forced  rotary  oscillations  on  the  steady  rotation  of  the  cylinder  is  also  presented; 
t.e.,  a  cylinder  rotary  motion  of, 

ff(t)  =  £2^  [Ho  +  a,  sin(27r/t)]  ,  (5.1.1) 

a 

where  is  the  angular  velocity  of  the  cylinder,  Qo  and  Qj  are  the  normalized 
amplitudes  of  the  steady  amd  oscillating  components  of  the  rotation,  and  /  is  the 
forcing  frequency. 
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5.2  Experimental  Setup 

The  experiments  documented  in  this  chapter  were  performed  in  the  GALCIT 
20"  X  20"  Free  Surface  Water  Timnel  (cf.  the  LSWC  used  for  the  wake  measurements 
in  Ch.  3). 

The  cylinder  used  for  the  near  wake  velocity  measurements  is  machined  Plex¬ 
iglas  cylinder,  1"  in  diameter  (  AR  =  18.7  ),  and  mounted  and  driven  by  the  appara¬ 
tus  described  in  Ch.  2.  The  command  signal  for  the  cylinder  velocity  was  generated 
by  a  computer-controlled  function  generator.  The  position  of  the  cylinder  output 
by  the  BEI  13-bit  absolute  position  encoder  was  read  using  the  computer’s  parallel 
interface. 

Velocities  were  measured  using  a  two-channel  (u,  u),  laser  Doppler  velocime- 
ter.  The  optics,  traverse  and  analog/digital  LDV  processor,  outlined  below,  are 
from  Lang’s  (1985)  thesis.  The  u -channel  was  only  used  to  measure  the  free  stream 
velocity.  The  v -channel  velocities  were  used  to  estimate  the  lift  coefficient.  The 
v -channel  was  Bragg  cell  frequency-shifted.  A  Lexel  Model  85  Argon-ion  laser  op¬ 
erating  in  single  line  mode  ( 514.5  nm)  at  about  200  mW  was  used  for  the  LDV. 
The  laser  beam  for  the  LDV  velocity  measurements  was  split  into  four  beams  of 
nearly  equal  intensity  using  two  cube  be£im  splitter  in  series  mounted  in  series  at 
-1-45“  and  —45“  respectively.  The  beams  for  the  u -channel  were  then  Bragg  cell 
frequency-shifted  to  allow  for  flow  reversals.  The  Bragg  cell  offset  between  the  two 
beams  was  300  kHz  for  the  data  presented  in  this  chapter.  Each  of  the  four  beams 
was  then  aligned  (to  insure  that  they  intersect  in  a  single  focal  volume)  by  pass¬ 
ing  them  through  wedge  prisms.  An  achromat  was  used  to  focus  the  beams  in  a 
focal  (measurement)  volume  located  in  the  midsptin  plane  of  the  FSWT.  The  light 
scattered  from  the  focal  volume  was  gathered  by  an  achromat  and  focused  on  a  pho- 
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todiode  with  an  integral  low  noise  pre-amplifier  designed  by  Dan  Lang.  The  signal 
from  the  pre-amplifier  was  then  band-pass  filtered  above  and  below  the  frequency 
corresponding  to  the  mean  flow  velocity  to  isolate  the  u -velocity  and  around  the 
Bragg  offset  to  isolate  the  u -velocity.  The  band-pass  ranges  were  chosen  to  pass 
instantaneous  velocity  fluctuations  and  were  non-overlapping. 

For  the  data  in  Figs.  5.2,  5.3,  and  5.4,  the  u-  and  t; -channel  outputs  from 
the  band  pass  filters  were  then  passed  thorough  a  pair  of  transimpedance  amplifiers 
(Lang  1985)  before  being  sent  to  the  analog/digital  LDV  processor  (Lang  1985). 
The  signals  were  then  processed  by  an  analog  and  digital  LDV  processor  and  the 
velocity  data  was  stored  to  disk  using  data  acquisition  software  written  by  Dan 
Lang.  See  Lang  (1985)  for  a  complete  description  of  the  analog/ digital  signal  pro¬ 
cessing  hardware.  Each  data  point  in  Fig.  5.2  is  an  average  over  about  10,000 
(instantaneous)  velocity  measurements. 

For  the  data  in  Figs.  5.5,  5.7,  5.8,  and  5.9,  a  pmr  of  tracking  phase-locked 
loops  designed  by  Dan  Lamg  and  Paul  Dimotedds  were  used  to  lock  TTL  square 
waves  to  the  dominant  u-  and  v -channel  frequencies  from  the  band-pass  filters. 
The  TTL  signals  were  then  read  by  a  (multiple)  counter-timer  bo2u-d  on  a  data 
acquisition  computer  amd  stored  on  disk  for  later  processing. 

Flow  visualization  was  accomplished  by  introducing  red  dye  into  the  flow 
ahead  of  the  cylinder  through  a  hypodermic  tube.  A  blue  filter  was  used  to  darken 
the  red  dye  marker  relative  to  the  white  background.  Photographs  were  taken  using 
a  35mm  camera.  Illumination  was  provided  from  behind  by  a  large  (backlit)  white 
sheet  of  paper  suspended,  outside  the  water  tunnel,  behind  the  model. 
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5.3  Lift  estimation  by  the  virtual  vortex  method 


Since  the  cylinder  in  the  present  experiments  experiences  no  net  rotational 
acceleration  it  cannot  shed  any  (net)  vorticity  into  the  wake,  in  the  mean.  It  follows 
that  the  net  circulation  in  the  wake  cannot  contribute  to  F.  As  a  consequence,  the 
vorticity  contributing  to  F  must  be  restricted  to  a  relatively  small  neighborhood 
around  the  cylinder.  This  suggests  that  the  effect  on  the  flow,  some  distance  from 
the  lifting  body,  can  be  approximated,  in  a  multipole  exp2msion  sense,  by  a  single 
virtual  vortex  of  strength  F.  It  is  also  reasonable  to  expect  that  this  virtual  vortex 
resides  somewhere  within  the  cylinder  boundaries. 


Ideally,  in  an  infinite  domain,  the  velocity  u(i,  y  =  0) ,  a  distance  x  upstreaim 
of  a  vortex  center  is 


F 

27rx 


(5.3.1) 


Using  Eqs.  1.2.5  and  5.3.1,  v  can  be  written  in  terms  of  the  lift  coefficient.  Cl. 
For  the  infinite  domain  case  this  can  be  written  as 


v(x)  _  Cl  X 
Uoo  2ir  a 


(5.3.2) 


Because  of  the  finite  height  of  the  test  section,  a  proper  account  should 
include  spatially  periodic  image  vortices  located  above  and  below  the  cylinder.  See 
Fig.  5.1.  Equation  5.3.1  is  then  replaced  by 


V 


(5.3.3) 


where  h  is  the  test  section  height  and  iq  is  the  streamwise  position  of  the  virtual 


Fig.  5.1  Diagram  of  periodic  image  vortices.  Only  the  image  vortices  immediately 
above  and  below  are  pictured.  The  shaded  circle  denotes  the  cylinder  (not 
to  scale),  and  the  smaller  circles  the  virtual  vortex  and  its  images. 


vortex.  •  Rearranging  Eq.  5.3.3  yields, 


(5.3.4) 


A  transverse  displacement  of  the  virtual  vortex  from  the  cylinder  axis,  j/q  ,  is  a  higher  order 
effect. 
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Fig.  5.2  Normalized  transverse  velocity  ahead  of  a  cylinder  executing  steady  ro¬ 
tation  (fli  =  0):  2hvlT  vs.  (i  —  xo)/a,  Rt  —  3.8  x  10^.  Qq  =  0.5  (  □ 
),  1.5  (  O  ),  2  (  A  ),  2.5  (  -h  ),  3  (  ^  ),  4  (  o  ),  5  (  *  ),  6  (  A  ),  8 
(  1  ),  10  (  —  ),  and  Hi  =  0.  The  dashed  line  is  csch(7r(z  —  xo)/h). 

or,  using  Eq.  1.2.5, 


V 


a 

2* 


Cl  csch 


(5.3.5) 


The  transverse  velocity  was  meatsured  at  several  locations  sdiead  of  the  cylin¬ 
der  axis,  for  a  range  of  0.5  <  fto  ^  10-  The  validity  of  Eq.  5.3.3  as  a  model  for 
the  decrease  in  the  velocity  ahead  of  the  cylinder  is  demonstrated  in  Fig.  5.2,  a 
plot  of  2hv/T  vs.  (x  —  xo)/a.  The  parameters  h,  a,  and  i  were  known  a  priori 
and  vfUoo  was  measured.  F  and  xq  were  determined  using  Eq.  5.3.5  and  a  version 
of  the  Levenberg-Marquardt  nonlinear  least-squares  fit  routine  from  Press  et  al. 
(1986).  Figure  5.3  compares  the  data  of  Reid  (1924)  with  Cl  determined  from  the 
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f’o 


Fig.  5.3  Ci  based  on  data  fit  in  Fig.  5.2  i2e  =  3.8  x  lO’,  fli  =  0,  AR  =  18.7  (□). 

Data  of  Reid  (1924)  AR  =  13.3,  Re  =  4.0  x  10^  (  Q  )>  -Re  =  5.6  x  10^ 
(  A  ),  Rc  =  7.9  X  10^  (  +  ),  Re  =  1.2  X  10®  (  — ♦  )•  Data  of  Prandtl 
=4.7,  Re  =  5.2  x  10^  ( - ). 

nonlinear  least  squares  curve-fit. 

The  curve-fit  data  also  show  that  the  centroid  of  the  virtual  vortex,  iq  <  is 
slightly  ahead  of  the  cylinder  axis.  See  Fig.  5.4.  The  centroid  appears  to  be  slightly 
ahead  of  the  cylinder  body,  «.e.,  xo/a>  1,  at  the  lowest  value  of  Do  (=0.5).  Ideally, 
this  should  not  occur.  The  velocities  being  measured  in  that  case,  however,  are 
small  (0.004  <  v/Uoo  <  0.02)  and  errors  of  this  order  cause  anomalies  in  the  fit 
parameters. 


Fig.  5.4  Centroid  of  the  virtual  vortex,  xo,  based  on  the  data  fit  in  Fig.  5.2. 

In  the  following  sections,  Cl  is  estimated  using  a  single  velocity  ahead  of  the 
cylinder,  and  taking  the  centroid  of  the  virtual  vortex  to  be  on  the  cylinder  axis. 
In  addition,  the  other  parameters  were  set  at  x/a  =  10,  /i/a  =  40,  x/h  =  0.25  and 
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5.4  Steady  rotation:  f}i=0. 

Figure  5.5  shows  the  lift  calculated  using  Eq.  5.3.5  with  iq  =  0  and  v/Uoo 
measured  at  x/a  =  10.  The  published  data  of  Pramdtl  (1925)  and  Reid  (1924)  are 
shown  for  reference.  Note  that  the  present  data  seem  to  overestimate  Cl  at  the 
lowest  values  of  flo  •  This  may  partially  be  attributable  to  the  fact  that  the  origin 
of  the  virtual  vortex  is  not  on  the  cylinder  axis.  However,  the  fact  that  the  data 
of  Prandtl  (1925)  and  Reid  (1924)  show  a  negative  Cl  »  at  the  lowest  values  of  Slo » 
suggests  that  the  lower  Reynolds  number  in  the  present  experiments  may  also  be 
a  factor.  Note  that  Tritton  (1977)  reports  that  a  sphere  also  experiences  a  lift  in 
the  “wrong”  direction  at  low  rotation  rates  and  higher  Re  because  of  a  turbulent 
separation  on  the  side  of  the  sphere  moving  opposite  to  the  flow.  As  can  be  seen 
from  the  data,  Fig.  5.5,  larger  cylinder  aspect  ratios  ( AR )  yield  larger  maximum 
lift  coefficients  (  Cl„„  )•  Compare  the  data  in  the  present  experiments  employing  a 
cylinder  with  a  larger  AR  =  18.7  with  the  data  of  Reid  (  AR  =  13.3)  and  Prandtl 
(Afi  =4.7). 

Recall  that  Goldstein  (1938),  in  ein  interpretation  of  an  intuitive  argument 
given  by  Prandtl  (1925),  suggested  that  Cl„„  =  4  7r  w  12.6.  The  present  data 
which  show  the  estimated  Cl  exceeding  Az  casts  doubt  on  this  proposed  upper 
limit.  Note  that  the  value  of  Cl  measured  at  the  highest  value  of  flo  exceed 
“Cl„„”  by  almost  25%.  In  addition,  the  trend  suggests  that  the  estimated  Cl 
can  be  made  even  larger  by  increasing  flo  beyond  12. 

It  should  be  noted  that  in  real  flow,  diffusion  and  unsteady  flow  processes 
can  transport  vorticity  away  from  the  cylinder  at  start-up,  weakening  Prandtl’s 
proposed  Cl„„  (Eq.  1.2.7).  A  more  plausible  explanation  is  that  three-dimensional 
(end)  effects  will  tend  to  reduce  the  mesm  spanwise  lift  measured  at  a  particular 
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Fig.  5.5  Comparison  with  published  data  of  Ct  va.  Qq  for  steadily  rotating  cylin¬ 
ders  whose  ends  are  flush  with  the  test  section  walls.  Present  data  (circles) 
Re  =  4.2  X  10^  ,  (squares)  Re  =  6.8  x  10^ ,  AR  =  18.7 ,  in  water.  Data  of 
Reid,  (dots)  Re  =  4.0  x  10^  <  Re  <  1.2  x  10®,  AR  =  13.3,  in  air.  Data 
of  Prandtl  (dashes).  Re  =  5.2  x  10^ ,  AR  =  4.7,  in  air. 

rotation  rate  to  below  that  which  could  be  attained  in  a  two-dimensional  flow. 
Thom  (1926)  showed  that  the  sectional  lift  coefficient  decreased  toward  the  ends  of 
a  rotating  cylinder.  It  is  interesting  that  Goldstein  also  relates  how  Prandtl,  in  one 
of  his  experiments,  associated  the  limiting  of  (=  4)  with  a  separation  of  the 

flow  from  the  side  walls  of  the  test  section.  To  remedy  this  situation,  Prandtl  added 
co-rotating  end  disks  to  the  cylinder  and  then  observed  an  increased  (=  10)- 

No  further  note  was  made  of  the  end  effects  for  that  configuration. 
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5.5  Rotation  with  forced!  oscillations:  fli  ^  0,  5/  =  0.7. 

For  the  data  presented  in  this  section,  the  forcing  Strouhal  number, 

2a/ 

Uoo  ' 

was  fixed  at  a  value  of  0.7,  while  the  steady  and  fiuctuating  components  of  the 
rotary  motion  were  varied. 


0  2  4  6  8  10  12 


Fig.  5.7  Comparison  of  Cl  vs.  Qq  data  for  several  amplitudes  of  oscillations, 
with  constant  forcing  frequency  and  Reynolds  number  (5/  =  0.7,  Re  = 
6.8  X  10^).  Squares;  flj  =  0,  crosses;  1.2,  and  diamonds:  2.3. 
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Fig.  5.6  Flow  visualization  for  various  mean  rotation  rates  ( flo  )■  With  oscillations 
(left.  =  2.1,  5/  =  0.7),  and  without  (right,  fli  =0),  cf.  stars  and 
circles  in  Fig.  5.8.  From  top  to  bottom,  the  mean  (normalized)  rotation 
rate  is  Uq  =  0,  1,  2,  3,  4.  and  5.  Re  =  4.1  x  10^.  Note  that  the  picture 
on  the  top  right  corner  is  the  unforced  (non-rotating)  case. 
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Fig.  5.8  Comparison  of  Ci  vs.  f2o  data  for  the  steady  case  and  an  oscillating 
case,  at  two  Reynolds  numbers.  Steady  rotation:  (circles)  Re  =  4.2  x  10* , 
fli  =  0,  (squares)  Re  =  6.8  x  10*,  Uj  =  0.  Net  rotation  with  oscillations, 
(stars)  Re  =  4.2  x  10*,  =2.1,  5/  =  0.7,  (diamonds)  Re  =  6.8  x  10*, 

=  2.3,  Sf  =  0.7. 

It  W21S  found  that  forced  rotary  oscillations  increased  Cl  for  0  <  fio  <  2.5 
and  decreased  it  for  2.5  <  flo  <  4.5.  See  Fig.  5.7.  It  is  no  coincidence  that  for 
the  steadily  rotating  case,  flj  =  0,  the  cylinder  has  noticeable  vortex  shedding  for 
flo  <  2.5,  and  none  for  flo  >  2.5.  Flow  visualization  (Fig.  5.6)  shows  that  for 
flo  <  2.5,  forced  oscillations  of  the  cylinder  help  close  the  wake,  creating  a  flow 
that,  on  average,  is  closer  to  potential.  In  contrast,  for  2.5  <  flo  <  4.5,  where 
the  wake  would  normally  close  with  steady  rotation  alone,  oscillations  hinder  this 
effect.  The  data  for  the  three  cases  presented  in  Fig.  5.7  coincide  for  flo  >  4.5,  i.e.. 
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the  oscillations  of  the  cylinder  have  little  effect  on  the  lift  beyond  this  point. 

There  is  a  cross-over  region  of  the  Ci  vs.  flo  curve  at  fio  between  2  and 
4.  The  data  in  Figs.  5.7  and  Fig.  5.8  suggest  that  Cl  in  that  region  is  not  a 
strong  function  of  the  forcing  amplitude  or  the  Reynolds  number.  While  there  is  an 
indication  in  Fig.  5.6  (c)  that  the  forced  oscillation  frequency  ( 5/ )  and  the  natural 
shedding  frequency  ( S„  )  are  close  below  the  cross-over  region  ( fto  <  2.5  ),  estimates 
of  Cl  varied  little  with  5/  in  that  region.  See  Fig.  5.9  (squares). 

As  can  be  seen  in  Fig.  5.8,  the  data  for  Re  =  4.2  x  10* ,  and  Re  =  6.8  x  10* 
are  close,  except  for  a  slightly  delayed  break  in  the  lift  curve  around  Cl  =  4.5  for 
Re  =  4.2  X  10* .  This  is  possibly  due  to  a  delayed  boundary  layer  separation  from 
the  test  section  walls  in  the  latter  case. 


CHAPTER  6 


Conclusions 


The  type  of  active  program  control  of  the  flow  described  in  the  present  ex¬ 
periments  relies  on  the  forces  between  the  cylinder  and  the  surrounding  flow.  It 
is  this  interaction  that  regulates  the  ejection  of  circulation  into  the  separated  flow 
region  and  actuates  the  entire  flow  fleld.  This  is  exemplified  by  the  control  of  the 
wake  structure  demonstrated  in  Chs.  3  -  5  as  well  as  the  up-stream  influence  char¬ 
acterized  in  Ch.  5.  In  addition,  the  way  this  type  of  control  exploits  the  natural 
tendency  to  have  vortical  structures  in  the  flow  is  in  contrast  to  the  more  traditional 
control  philosophy  that  unsteadiness  is  something  that  should  be  canceled,  e.^.,  the 
feed-forward  control  employed  by  Liepmann,  Brown  &  Nosenchuck  (1982)  and  Liep- 
mann  &  Nosenchuck  (1982),  auid  the  feedback  control  employed  by  Ffowcs- Williams 
Zhao  (1989). 

In  Ch.  3,  the  efficacy  of  oscillatory  cylinder  rotation  as  an  actuation  mech¬ 
anism  for  actively  controlling  the  cylinder  wake  was  examined.  It  wzis  found  that 
considerable  control  could  be  exercised  over  the  structure  in  the  wake  with  such 
forcing.  In  pzirticular,  a  large  increase,  or  decreaise,  in  the  resulting  displacement 
thickness,  estimated  cylinder  drag,  and  associated  mixing  with  the  free  stream  can 
be  achieved,  depending  on  the  frequency  and  amplitude  of  oscillation.  Not  sur¬ 
prisingly,  the  present  results  show  that  working  in  a  control  domain  in  which  the 
structures  shed  axe  synchronous  with  the  forcing  provides  the  greatest  control  au¬ 
thority  over  the  wake  structure.  In  addition,  while  these  results  were  obtained  for 
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moderate  Reynolds  numbers,  i.  e.  Re  ~  10^ ,  there  are  indications  that  the  descrip¬ 
tion  of  the  flow  phenomena  presented  here  is  qualitatively  the  ssune  over  a  large 
range  of  Reynolds  numbers.  The  streamline  data  in  Ch.  4  show  that  under  control 
(forcing),  the  vortical  structures  are  formed  very  close  to  the  cylinder.  This  is  an 
indication  that  the  formation  of  the  wake  structure,  under  control,  has  more  to  do 
with  the  controlled  ejection  of  circulation,  and  the  associated  forces  on  the  body, 
than  with  the  stability  characteristics  of  some  generic  unforced  wake.  This  forced 
shedding  is  in  contrast  to  the  natural  (unforced)  shedding  of  vortices,  where  the 
structures  form  only  after  the  vorticity  has  been  shed  into  the  wake,  from  more-or- 
less  steady  separation  points.  Hence,  it  is  possible  that  the  stability  characteristics 
of  the  wake  can  play  a  more  significant  role  in  the  unforced  case  than  in  the  forced. 
In  addition,  it  is  shown  that  under  control,  very  little  of  the  circulation  introduced 
at  the  cylinder  wall  appears  in  the  wake.  Evidently,  most  of  the  positive  and  nega¬ 
tive  vorticity  introduced  into  the  boundary  layer  fluid  combines  before  the  fluid  is 
allowed  to  enter  the  wake. 

The  virtual  vortex  method  explored  in  Ch.  5  was  found  to  be  successful  in 
characterizing  flow  ahead  of  the  cylinder.  In  addition,  the  method  yields  estimates 
for  Cl  which  which  agree  with  the  data  published  by  Thom  (1924)  for  lift  coeffi¬ 
cients  greater  than  2,  i.e.,  flo  >  1-  While  the  method  seems  to  overestimate  Cl 
at  the  smellier  values  of  flo  <  1,  it  is  not  certain  whether  this  is  due  to  a  fail¬ 
ure  in  the  method,  a  difference  in  Re  between  the  experiments,  or  a  pronounced 
variation  of  the  sectional  lift  coefficient  along  the  cylinder  span  (Thom  1926).  The 
Cl  calculated  at  the  larger  values  of  Jlo  were  also  found  to  be  greater  than  in  the 
published  data,  but  this  was  probably  due  to  the  larger  aspect  ratio  of  the  cylinder 
in  the  present  experiments.  In  addition,  the  maximum  (obtainable)  lift  coefficient 
proposed  by  Prandtl  (1925),  was  exceeded  in  the  present  experiments,  possibly  be¬ 
cause  Prandtl’s  arguments  neglected  unsteady  effects.  Finally,  it  was  found  that  the 
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addition  of  forced  rotary  oscillations  to  the  steady  rotation  of  the  cylinder  increased 
Cl  in  cases  where  the  wake  would  normally  be  separated  (in  the  steadily  rotating 
case),  but  decreased  it  otherwise. 


APPENDIX  A 


Forces  on  a  body 


This  appendix  was  extracted  form  notes  collected  over  the  course  of  the 
present  investigations.  This  particular  segment  of  the  notes  was  chosen  because  it 
lends  support  to  the  “virtual  vortex  method”  used  in  Ch.  5.  In  addition,  it  shows 
how  the  unsteady  circulation  around  the  body  is  related  to  the  unsteady  forces. 


A.l  Forces  in  three-dimensional  flow 


The  ’unsteady  force  on  a  body  moving  through  an  incompressible  flow  (For 
steady  case  see  Milne-Thomson  1968,  p.  672),  can  be  written  as 

Fb  =  ~P^J  +  J  [-pn- pu{u-n)  +  fiu;  X  n]dS  ,  (A. 1.1) 

where,  V  is  a  volume  surrounding,  but  not  including,  the  body,  cr  is  the  surface 
enclosing  V  and  the  body,  p  is  the  fluid  density,  u  is  the  velocity,  p  is  the  pressure, 
p  is  the  viscosity,  u)  =  V  x  u  is  the  vorticity,  and  n  is  the  outward  normal  to  the 
surface  a.  See  Fig.  A.l. 


If  0 ,  the  Bernoulli  constant,  is  defined  as 


0  =  (Uoo  •  Uoo)  , 

Eq.  A.  1.1  cam  be  rewritten  ais 

Fs  =  -p—J  udV  +  y  |-/3n -I- i  p  (u  •  u)n  -  pu(u  •  n)  +  pw  X  n 


(A.1.2) 

dS  . 
A.1.3) 
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Fig.  A.l  3-D  Geometry. 

Note  that  the  nonlinear  terms  in  Eq.  A.1.3  can  be  written  (Milne- Thomson  1968, 
p.68)  as 

J  u  (u  •  n)  <i5  —  J  ^  {u  -  u)  ndS  =  J  wxudV  -h  J  u(V  •  u)  dV  .  (A. 1.4) 

Then,  noting  also  that  the  laist  term  in  Eq.  A.1.4  is  zero  in  incompressible  flow,  we 
have 

Fb  =  — p-T  f  ^dV  4-  /  [— /3n  4- X  n]  rfS  +  P  f  uxufdV  .  (A. 1.5) 
dt  Jy  Jtr  Jv 


In  order  to  extract  the  “classic”  Lift  =  pUooT  from  the  momentum  equation, 
it  is  useful  to  subtract  the  free-stream  velocity,  Uo© ,  from  the  velocity  field.  Let 


OO  • 


Uo 


u  -  U 


(A.1.6) 
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Noting  that 

u;  =  Vxu  =  Vx  (uo  +  Uoo)  =  V  X  uo 

and  also  that  the  vector  circulation  may  be  defined  as 

r  =  f  udV  =  /  V  X  Uo  dV  , 

Jv  Jv 

the  last  term  in  Eq.  A. 1.5  can  be  written  as 

I  u  X  uj  dV  =  I  (Uoo  +  Uo)  xuj  dV 
Jv  Jv 

=  Uoo  X  /  u>  dV  +  I  Uo  X  u>  dV 
Jv  Jv 

=  Uoo  X  r  +  /  Uo  X  (V  X  Uo)  dV 

Jv 


(A.1.7) 


(A.1.8) 


(A.1.9) 


It  is  desirable  to  write  the  force  in  terms  of  surface  integrals,  so  note  that 
re-applying  the  identity  (Eq.  A. 1.4)  to  uo ,  with  V •  uo  =  0,  yields  for  the  last  term 
in  Eq.  A.1.9 

J  Uo  X  (V  X  Uo)  dV  =  -  y  Uo  (uo  •  n)  dS  -f  J  ^  ■  uo)  ndS  .  (A. 1.10) 

For  ease  of  notation  define  the  surface  integral 

E  ^  ^  _  I  /  (A.1.11) 

y<r  Uoo  \Uoo  J  ha,  2  \Uoo  UooJ  ha 

such  that  Eq.  A.  1.10  can  be  written  as 

Uo  X  (V  X  Uo)  dV  =  —pUl^baJl  (A. 1.12) 


L 


where  2  6  a  *  is  the  frontal  area  of  the  body.  This  then  yields  for  the  force  on  a 
body  in  incompressible  flow, 

Fb  =  /  uodV  +  /  [-/3n  + /iu>  X  n]  dS  -I-  pU*,  x  F  +  pUlJ)a'L 

dt  Jy  J„ 

(A.1.13) 


*  For  the  case  of  a  circular  cylinder,  a  is  the  radius  and  6  is  the  span. 


Fig.  A.2  2-D  Geometry. 
A.2  Forces  in  two-dimensional  flow 


Writing  this  for  two-dimensionzd  flow, 


r  — *  r  6  e  j  ,  W  — *  U}  Bx  ,  Uqo  ^  ®r  1 


(A.2.1) 


yields  the  force  per  unit  span, 

Fb /6  =  -p-T  f  udA-  f  0ndl  -  pUooT By  ^  (  uBx  x  ndl  +  pUl^a S2_d  , 
dt  Ja  Jc  Jc 

(A.2.2) 

where  A  is  the  two-dimensional  cross-sectional  area  of  the  volume  V,  C  is  the 
boundary  of  A  (See  Fig.  A.2),  and 

_  /■  Uo  /  uo  \  dl  1  /■  /  Uo  uo  N  dl 

"  JcUZK^  ’')  ^  2  ic  Uoo  '  "  ■ 


n 


a 


(A.2.3) 
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In  terms  of  the  force  coefficient, 

Fb 


Cfj.D  —  1 


\pUl,2ab 


a  d  r  u  dA  f 

~~  UZ'dtjy,uZ~^  ~  Jc 


C  pUl,a^  a 


r  2  f  ua  dl  ^ 

Tr — +  '5"/  77-exXn—  +  S2-D 
C/oofl  Re.  Jc  Uoo  a 


A. 3  Lift  coefficient 


Calculating  the  lift  coefficient, 


Cl  =  Cf,. 


D  ^y 


a  d 

uZ 
r 

U6oa 


d  f  V  dA  f 
dt  Jc  Uoo  0}  Jc 

,  2  f  ua  dl  By 

+  -^  /  yz - -  2^2-D 

Re  Jc  Uoo  a 


0  ,  sdl 

-—5 — (n-ey)  — 

C  pUlo  a  a 


-y  » 


and  following  G.  I.  Taylor  (1925),  the  /?  integral  is  eliminated  by 
wake  vertically,  so  that  (n  •  e,)  =  0  (recall  that,  by  construction, 
the  wake).  Then 

a  d 


Cl  =  - 


where 


U^dt 


f  V  dA  r  2  f  loa  dy  ^ 

Jc'vZ'^  ~  UZi  ^  !cVZ~^  ~ 


fie  = 


is  the  Reynolds  number. 


Averaging  Eq.  A. 3.2  yields 


Cl  =  - 


U^a 


T  f 

Re  y*.ke  \UooJ  a 


(A.2.4) 


(A.3.1) 

traversing  the 
(3  =  0  outside 

)  •  ey  ,  (A.3.2) 


(A.3.3) 
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Because  F  tends  to  a  constant  outside  a  small  neighborhood  of  the  body  {e.g., 
Thom  1931)  and  Cl  is  a  property  of  the  body,  it  is  apparent  that 


2 

Rt 


dy 


—  S2-D  •  ®y 


const.  , 


and  since  S  approaches  zero  for  large  contours,  it  is  proposed  that 

2  f  /  ujd\  dy  2  Au 

/wake  ~d  ~  ~  — — 

for  any  contour  about  the  body  that  traverses  the  wake  vertically,  where  Au  is  the 
velocity  jump  across  the  wake  at  a  large  distance  from  the  body.  Equation  A.3.3 
can  then  be  written  as 


Cl  =  - 


r 

Uooa 


(A.3.4) 


For  even  moderate  Re,  the  second  term  in  Eq.  A.3.4  may  be  neglected,  since 
Au/Uoo  will  also  tend  to  be  small. 


APPENDIX  B 


Streamlines  averaged  at  constant  forcing  phase. 


Forcing  function: 

where  fii  =2,  5/  =  0.94,  and  i?c  =  2.1  x  10^  for  the  present  experiments.  In  the 
figures,  a  negative  rotation  is  in  the  clockwise  direction.  See  Ch.  4  for  a  description 
of  these  data. 


X 

Fig.  B.l  Time- averaged  streamlines. 


constant  p 

78. 
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A  method  for  estimating  the  mean  lift  of  a  rotating  cylinder  is  presented. 
Estimates  based  on  this  method  compare  favorably  with  similar  data  published  for 
steadily  rotating  cylinders.  It  was  also  found  that  the  addition  of  forced  rotary 
oscillations  to  the  steady  rotation  of  the  cylinder  increases  the  lift  coefficient  in  the 
cases  where  the  wake  would  normally  be  separated  in  the  steadily  rotating  case,  but 
decreases  it  otherwise. 


1.  Introduction 

A  rotating  cylinder  moving  in  a  uniform  stream  experiences  a  force  normal 
to  the  direction  of  motion.  Goldstein  (1938)  makes  several  historical  references 
for  both  rotating  spheres  and  cylinders,  crediting  Magnus  (1853)  with  the  first 
laboratory  experiments  examining  the  lift  on  a  rotating  cylinder.  Early  in  this 
century,  experiments  on  a  circular  cylinder  rotating  about  its  axis  in  a  uniform  flow 
were  performed  by  Reid  (1924),  Prandtl  (1925),  Thom  (1926,  1931),  for  example. 
Their  measurements  were  for  Reynolds  numbers  in  the  range, 

4  X  10^  <  Re  =  <  1.2  X  10®  ,  (1) 

where  a  is  the  cylinder  radius,  Uoo  is  the  free  stream  velocity,  and  v  is  the  kine¬ 
matic  viscosity.  More  recently,  experiments  and  simulations  have  been  performed 
at  low  Reynolds  numbers  by  Badar,  et  al.  (1989),  Ingham  &  Tang  (1990),  and 
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Tang  Sz  Ingham  (1991),  for  the  steady  flow  past  a  rotating  cylinder,  and  at  higher 
Reynolds  numbers  of  10^  and  10'*  by  Badar,  et  al.  (1990)  for  the  flow  past  a  cylinder 
impulsively  started  in  both  rotation  and  translation. 


As  is  well  known,  the  mean  lift  coefficient  of  a  cyi'.nder  can  be  written  as 

Cl  =  —.r2  =  -  77 —  ^  (2) 

pU^Q,  UqqO, 

where  p  is  the  fluid  density,  L  is  the  lift  per  unit  span,  and  F  is  the  mean 
circulation  taken  around  a  contour  enclosing  the  lifting  body.  See  for  example  Taylor 
(1925),  Thwaites  (1960,  V.9).  That  this  circulation  could  be  measured  around 
contours  close  to  the  cylinder  was  shown  experimentally  by  Thom  (1931). 


In  the  potential  flow  prototype  of  a  rotating  cylinder  in  a  uniform  free  stream, 
the  rotation  of  the  cylinder  is  modeled  by  a  point  vortex  of  strength  F  at  its  center. 
When  modeling  a  flow,  F  should  be  considered  as  a  parameter  to  be  adjusted  so  that 
the  flow  is  properly  represented;  e.</.,  to  satisfy  the  Kutta  Condition.  Proposing 
that  Cl  can  be  made  arbitrarily  large  by  arbitrarily  increasing  F  is  not  helpful;  F 
is  the  mean  lift  (c/.  Eq.  2). 


Prandtl  (1925)  argued  that  the  maximum  circulation,  Fn,ax »  which  could  be 
realized  about  the  rotating  cylinder  was  equal  to  the  circulation  at  which  the  up¬ 
stream  and  downstream  stagnation  points  join  on  the  bottom  of  the  cylinder.  From 
the  potential  flow  analog  of  that  case,  that  corresponds  to  a  value  of 

Finax  ~  .  (3) 

See  also  the  discussion  in  Goldstein  (1938  I,  pp.  81-84).  Beyond  this  point,  Prandtl 
argued  that  no  circulation  could  be  shed  by  the  cylinder  to  infinity,  as  it  is  ramped 
up  from  rest  to  a  constant  rotation  rate,  fixing  the  total  circulation  at  infinity,  and 
the  lift  of  the  rotating  cylinder  thereby.  If  this  argument  were  correct,  the  maximum 
steady-state  lift  coefficient  that  can  be  realized  would  be  given  by, 

=  47r  «  12.6  .  (4) 

Note  that  Prandtl  relates  the  real  and  potential  flow  case  of  the  rotating  cylinder  not 
by  matching  the  circulation  calculated  by  taking  a  contour  around  the  cylinder  at 
the  surface,  as  one  would  exjiect.  but  by  matching  the  circumferential  velocities 
in  the  two  cases. 
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The  present  experiments  examine  the  effect  of  cylinder  rotation  on  the  flow 
ahead  of  the  cylinder.  In  order  to  assess  the  effects  of  forcing,  the  flow  ahead  of  the 
cylinder  was  described  in  terms  of  the  strength  and  position  of  a  virtual  vortex.  The 
strength  of  the  virtual  vortex  was  then  related  to  the  lift  coefficient  of  the  cylinder. 
Results  obtained  for  a  steadily  rotating  cylinder  (no  forced  oscillations)  compare 
favorably  with  similar  data  published  in  the  literature.  It  was  also  found  that  the 
addition  of  forced  rotary  oscillations  to  the  steady  rotation  of  the  cylinder  helped 
increase  the  lift  coefficient,  in  cases  where  the  wake  would  normally  be  separated 
in  the  steadily  rotating  case,  and  decrease  it  otherwise.  It  is  believed  that  this  is 
the  first  investigation  of  the  flow  past  a  cylinder  with  both  net  rotation  and  forced 
oscillations. 


2.  Experimental  setup 

The  present  experiments  were  performed  in  the  GALCIT  20"  x  20"  Free 
Surface  Water  Tunnel*  (see  Ward  1976).  A  machined  Plexiglas  cylinder,  1"  in 
diameter  and  an  aspect  ratio  of  18.7  was  supported  10"  above  the  bottom  of  the 
water  channel,  between  0.5"  thick  Plexiglas  fairings  placed  flush  to  the  side  walls 
of  the  channel.  Power  was  transferred  from  the  motor  to  the  cylinder  using  steel 
reinforced  neoprene  timing  belts.  These  were  entirely  enclosed  in  the  fairings  and 
did  not  interfere  with  the  flow.  The  cylinder  was  driven  by  a  high  performance 
DC  motor  from  PMI  Motion  Technologies  (JR16M4CH  Servodisk^^ ),  capable  of 
tracking  an  arbitrary  command  signal  within  the  bandwidth  and  slewing  rates  that 
were  investigated.  The  command  signal  for  the  cylinder  angular  velocity  was  gen¬ 
erated  by  a  computer-controlled  function  generator.  The  angular  position  of  the 
cylinder  was  measured  using  a  BEI  Motion  Systems  Co.  13-bit  absolute  position 
encoder  connected  to  the  motor  driveshaft. 

Velocities  were  measured  using  a  two-channel  (u,  u),  laser  Doppler  velocime- 
ter.  The  optics,  traverse  and  analog/digital  LDV  processor,  outlined  below,  are 


The  setup  of  the  present  experimeuts  is  similar  to  that  used  in  Tokumaru  &  Dimotakis  (1991), 
and  Tokumaru  (1991). 
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as  were  employed  by  Lang  (1985).  The  u-cliannel  was  only  used  to  measure  the 
free  stream  velocity.  The  u-channcl  velocities  were  used  to  estimate  the  lift  coeffi¬ 
cient.  A  Lexel  Model  85  Argon-ion  laser  operating  in  single-line  mode  (514.5  nm) 
at  about  200  mW  was  used  for  the  LDV.  The  laser  beam  for  the  LDV  velocity  mea¬ 
surements  was  split  into  four  beams  of  nearly  equal  intensity  using  two  cube  beam 
splitters  in  series,  at  -f  45°  and  —45° ,  respectively.  The  beams  for  the  u -channel 
were  then  Bragg-cell  frequency-shifted  to  accommodate  flow  reversals.  The  Bragg¬ 
cell  frequency  offset  between  the  two  beams  was  300  kHz  for  the  present  data.  For 
alignment,  to  insure  that  the  beams  intersected  in  a  single  probe  volume,  each  of 
the  four  beams  was  passed  through  wedge  prisms.  An  achromat  was  used  to  focus 
the  beams  in  a  focal  (measurement)  volume  located  in  the  midspan  plane  of  the 
FSWT.  Light  scattered  from  the  focal  volume  was  gathered  by  two  achromats,  in 
series,  and  focused  on  a  photodiode  with  an  integral  low  noise  pre-amplifier.  The 
signal  from  the  pre-amplifier  was  then  band-pass  filtered  above  and  below  the  fre¬ 
quency  corresponding  to  the  mean  flow  velocity  to  isolate  the  u -velocity  and  around 
the  Bragg  offset  to  isolate  the  u -velocity.  The  band-pass  ranges  were  chosen  to  pass 
instantaneous  velocity  fluctuations  and  were  non-overlapping. 

For  the  data  in  Figs.  2-4,  the  signals  from  the  band-pass  filters  were  processed 
by  an  LDV  signal  processor  (Lang  1985)  and  the  velocity  data  were  stored  on  disk. 
Each  data  point  in  Fig.  2  is  an  average  over  about  10,000  (instantaneous)  velocity 
measurements.  For  the  data  in  Fig.  6.  a  pair  of  tracking  phase-locked  loops  was 
used  to  lock  TTL  square  waves  to  the  dominant  u-  and  u-channel  frequencies  from 
the  band-pass  filters.  The  TTL  signals  were  then  read  by  a  (multiple)  counter- timer 
board  on  a  data  acquisition  comi)uter  and  stored  on  disk  for  later  processing. 

Flow  visualization  was  accomplished  by  introducing  red  dye  into  the  flow 
ahead  of  the  cylinder  through  a  hypodermic  tube.  A  blue  filter  was  used  to  darken 
the  red  dye  marker  relative  to  the  white  background.  Photographs  were  taken  using 
a  35  mm  camera.  Illumination  was  provided  from  behind  by  a  large  (backlit)  white 
sheet  of  paper  suspended,  outside  the  water  tunnel,  behind  the  model. 
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3.  Cylinder  motion 


The  angular  motion  of  the  cylinder  in  the  present  experiments  can  be  de¬ 
scribed  by  a  single  equation;  i.e.. 


=  Q.0  +  ill  sin(27r/t)  , 


where 


Q  = 


6  a 

uZ 


(5) 


(6) 


is  the  normalized  rotation  rate  of  the  cylinder,  a  is  the  cylinder  radius,  Ucx>  is  the 
free-stream  velocity,  t  is  the  time,  6  is  the  angular  velocity  of  the  cylinder,  /  is 
the  forcing  frequency,  and,  Q.q  and  Qi  are  amplitudes  of  the  steady  and  harmonic 
components  of  the  cylinder  motion.  The  normalized  forcing  freqtiency  is 

2a/ 


5,  = 


(7) 


i.e.,  the  forcing  Strouhal  number.  The  Reynolds  numbers  in  the  present  experiments 
are  Re  =  3.8  x  10^  and  6.8  x  10^ . 


4.  Lift  estimation  by  the  virtual  vortex  method 

The  cylinder  in  the  present  experiments  experiences  no  net  rotational  acceler¬ 
ation  after  it’s  initial  start-up.  Hence,  once  steady  state  flow  has  been  established, 
it  can  not  shed  any  additional  vorticity  into  the  wake,  in  the  mean.  It  follows 
that,  some  time  after  start-up,  the  circulation  in  the  wake  cannot  contribute  to  the 
mean  circulation,  F,  taken  around  the  cylinder.  As  a  consequence,  the  vorticity 
contributing  to  F  must  be  restricted  to  a  relatively  small  neighborhood  around  the 
cylinder.  This  suggests  that  the  effect  on  the  flow,  some  distance  from  the  lift¬ 
ing  body,  can  be  approximated,  in  a  multipole  expansion  sense,  by  a  single  virtual 
vortex  of  strength  F . 

Ideally,  in  an  infinite  domain,  the  transverse  velocity  v{x,y  =  0),  a  distance 
X  upstream  of  a  vortex  center  is 

F 


v{x,y  =  0)  =  - 


2'kx 


(8) 
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Fig.  1  Diagram  of  periodic  image  vortices.  Only  the  image  vortices  immediately 
above  and  below  are  pictured.  The  shaded  circle  denotes  the  cylinder  (not 
to  scale),  and  the  smaller  circles  the  virtual  vortex  and  its  images. 

Using  Eqs.  2  and  8,  v  can  be  written  in  terms  of  the  lift  coefficient,  Ci-  For  the 
infinite  domain  case  this  can  be  written  as 

=  0)  _  /gx 

Uoo  27r  a 

Because  of  the  finite  height  of  the  test  section,  however,  a  proper  account  should 
include  spatially  periodic  image  vortices  located  above  and  below  the  cylinder.  See 
Fig.  1.  Equation  8  is  then  replaced  by 


v[x,y  =  0)  =  -  ^  esch  ^TT 


X  —  Xq 

h 


where  h  is  the  test  section  height  and  Xq  is  the  streamwise  position  of  the  virtual 
vortex.  Rearranging  Eq.  10  yields, 

2hv{x,y  =  0)  ^  ,,,, 


=  —  csch  TT 


or.  using  Eq.  2. 


=  0)  =  ±Cr  csch 

Uoo  2h 
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Fig.  2  Normalized  transverse  velocity  ahead  of  a  cylinder  executing  steady  rotation 
(f2i  =  0);  2/iu/r  vs.  {x  —  xo)fa.  Re  =  3.8  x  10^ .  Qq  =  0.5  (  □  ),  1.5  ( 
O  ),  2  (  A  ),  2.5  (  +  ),  3  (  -.  ),  4  (  o  ),  5  (  O,  6  (  A  ),  8  (  I  ),  10  ( 
—  ),  and  fli  =  0.  The  dashed  line  is  csch(7r(x  —  xo)fh). 

A  more  general  form  of  Eq.  12  which  allows  for  both  a  streamwise  {  xq  )  and  trans¬ 
verse  ( yo )  displacement  of  the  virtual  vortex  is 


v{x,y  =  0)  ^  ^  cos  (tt  sinh  (tt 

C/qo  2/i  ^  sin^  (’’’  ^)  +  sinh^  (tt  ) 

Note  that  when  mecisuring  v  along  the  x  ( y  =  0 )  axis,  a  transverse  displacement 
of  the  virtual  vortex  from  that  axis,  yo  i  is  a  higher  order  effect  than  a  streamwise 
displacement,  xq,  t.e.. 


(14) 
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Fig.  3  Cl  based  on  data  fit  in  Fig.  2  ile  =  3.8  x  10^,  f2i  =0,  >!l  =  18.7  (□). 
Data  of  Reid  (1924)  M  =  13.3,  Re  =  4.0  x  lOM  O  ),  -Re  =  5.6  x  10“  ( 
A  ),  /2e  =  7.9  X  10“  (  +  ),  Re  =  1.2  x  10^  (  — »  ).  Data  of  Prandtl  (1925) 
>R  =  4.7,  Re  =  5.2  x  10“  ( - ). 


5.  Steady  rotation:  Di  =0. 

To  test  the  utility  of  Eqs.  10  -  12,  the  transverse  velocity  was  measured  at  six 
locations  ahead  of  the  cylinder  axis,  for  a  range  of  0.5  <  Dq  ^  10  ( i?e  =  3.8  x  10^ ). 
The  validity  of  this  model  for  the  transverse  ( v )  velocity  ahead  of  the  cylinder  is 
demonstrated  in  Fig.  2,  a  plot  of  2/ii»/r  vs.  (x  —  XQ)la.  The  parameters  /i,  a, 
and  X  were  known  a  priori  {/i/a  =  40),  and  v/Uoo  was  measured.  F  and  xq  were 
determined  using  Eq.  11  and  a  version  of  the  Levenberg-Marquardt  nonlinear  least- 
squares  fit  routine  from  Press  et  al.  (1986).  Figure  3  compares  the  results  of  Reid 
(1924)  and  Prandtl  (1925)  with  Cl  determined  from  this  nonlinear  least  squares 
curve-fit. 
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Fig.  4  Centroid  of  the  virtual  vortex,  xq,  based  on  the  data  fit  in  Fig.  2. 


It  may  appear  that  the  present  results  measured  at  Re  =  3.8  x  10^  overes¬ 
timate  Cl  at  the  lowest  values  of  flo-  In  fact,  the  results  of  Prandtl  (1925)  and 
Reid  (1924),  at  Re  >  4  x  10^,  show  a  negative  Cl  at  the  lowest  values  of  flo-  Ex¬ 
periments  at  much  lower  Re  {Re  <  100,  Badar,  et  al.  1989,  Tang  &  Ingham  1990) 
yield  a  Cl  somewhat  higher  tliaxi  the  present  experiments.  Tritton  (1977)  reports 
that  a  sphere  also  experiences  a  lift  in  the  “wrong”  direction  at  low  rotation  rates 
and  higher  Re  because  of  a  turbulent  separation  on  the  side  of  the  sphere  moving 
opposite  to  the  flow.  This  suggests  that  the  discrepancy  at  low  Qq  may  well  be 
attributable  to  Reynolds  number  effects. 

As  can  be  seen  from  Fig.  3,  larger  cylinder  aspect  ratios  (Al)  yield  larger 
maximum  lift  coefficients  (Cl^^^  )•  Compare  the  results  in  the  present  experiments 
employing  a  cylinder  with  a  larger  Al  =  18.7  with  those  of  Reid  ( Al  =  13.3)  and 
Prandtl  (Al  =  4.7). 


The  curve-fit  results  also  show  that  the  centroid  of  the  virtual  vortex,  xq  ^  is 
slightly  ahead  of  the  cylinder  axis.  See  Fig.  4.  In  fact,  this  centroid  appears  to 
be  slightly  upstream  of  the  cylinder  body,  i.e.,  xo/a  >  1,  at  the  lowest  value  of 
flo  (=0.5).  This  may  counter  intuition,  as  very  little  of  the  vorticity  in  the  flow  is 
situated  upstream  of  the  cylinder.  As  can  be  seen  from  Fig.  5  and  Eq.  14,  allowing 
for  a  virtual  origin  away  from  the  cylinder  axis  is  a  higher  order  correction. 


0  2  4  6  8  10 
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Fig.  5  Lift  estimates  holding  the  virtual  origin  coincident  with  the  cylinder  axis, 
A’s,  and  allowing  it  to  vary,  *’s. 

Recall  that  Goldstein  (1938).  in  an  interpretation  of  an  intuitive  argument 
given  by  Prandtl  (1925),  suggested  that  =  4  7r  ss  12.6.  The  present  results, 

that  show  the  estimated  Cl  exceeding  47r,  cast  doubt  on  this  proposed  upper  limit. 
Note  that  the  value  of  Cl  measured  at  the  highest  value  of  flo  exceeds  ” 

by  more  than  20%.  In  addition,  the  trend  suggests  that  the  estimated  Cl  can  be 
made  larger  by  further  increasing  flo-  or  Al.  or  both. 
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Diffusion  and  unsteady  flow  processes  can  transport  vorticity  away  from  the 
cylinder  at  start-up,  weakening  Prandtl’s  proposed  (Eq.  4).  A  more  plausible 

explanation  is  that  three-dimensional  (end)  effects  will  tend  to  reduce  the  mean 
spanwise  lift  measured  at  a  particular  rotation  rate  to  below  that  which  could  be 
attained  in  purely  two-dimensional  flow.  Thom  (1926)  showed  that  the  sectional 
lift  coefficient  decreased  toward  the  ends  of  a  rotating  cylinder.  It  is  interesting 
that  Goldstein  also  relates  how  Prandtl,  in  one  of  his  experiments,  associated  the 
limiting  of  (=  4)  with  a  separation  of  the  flow  from  the  side  walls  of  the 

test  section.  To  remedy  this  situation,  Prandtl  added  co-rotating  end-disks  to  the 
cylinder  and  then  observed  an  increased  (=  10).  No  further  discussion  of 

end-effects  for  that  configuration  was  included. 


6.  Rotation  with  forced  oscillations:  #  0,  Sf  =  0.7. 

For  the  data  presented  in  this  section,  the  forcing  Strouhal  number,  5/,  is 
fixed  at  0.7,  while  the  steady  and  fluctuating  components  of  the  rotary  motion  are 
varied.  The  Reynolds  number  is  somewhat  higher  than  in  the  previous  section  at 
6.8  X  10^ .  Cl  is  estimated  using  Eq.  12  and  the  mean  velocity  at  a  single  point 
ahead  of  the  cylinder,  v(x/a  =  10).  Because  xq/x  •C  1,  the  centroid  of  the  virtual 
vortex  is  taken  to  be  on  the  cylinder  axis,  i.e.,  xo/h  w  0.  See  Eq.  14.  In  addition, 
the  ratio  of  the  cylinder  radius  to  the  water  tunnel  depth  is  a/h  =  0.025. 

It  was  found  that  forced  rotary  oscillations  increased  Cl  for  0  <  Hq  <  2.5 
and  decreased  it  for  2.5  <  flo  <  4.5.  See  Fig.  6.  For  the  steadily  rotating  case, 
fli  =  0,  the  cylinder  has  noticeable  vortex  shedding  for  Dq  <  2.5,  and  none  for 
Do  >  2.5. 


Flow  visualization  (c/.  Figs.  7-9)  reveals  that,  for  Do  <  2.5,  forced  oscilla¬ 
tions  of  the  cylinder  help  close  the  wake,  creating  a  flow  that,  on  average,  is  closer 
to  potential  (the  forced  oscillations  may  also  decrease  spanwise  variations  in  the 
flow).  In  contrast,  for  2.5  <  Do  <  4.5,  where  the  wake  would  normally  close  with 
steady  rotation  alone,  oscillations  have  the  opposite  effect.  The  data  for  the  three 
cases  presented  in  Fig.  6  coincide  for  Do  >  4.5,  i.e..  the  oscillations  of  the  cylinder 


12 


16 

14 

- 1 - 1 - 1 - 1 - 

12 

- 

10 

nAO 

- 

cy  8 

O 

l<JO 

CKO 

□ 

_ 

. 

6 

6^° 

4 

O 

2 

OA  ° 

_  OA  O 

01 

^ _ . 

0  2  4  6  8  10 


Oo 


Fig.  6  Comparison  of  Ci  vs.  Qq  data  for  several  amplitudes  of  oscillations,  with 
constant  forcing  frequency  and  Reynolds  number  (5/  =  0.7,  Re  =  G.S  x 
lO^).  □:  =0,  A:  1.2.  and  Q  :  2.3. 

have  little  effect  on  the  lift  beyond  this  point.  There  is  also  a  cross-over  region  of 
the  Cl  vs.  flo  curve  at  Qq  between  2  and  4,  where  Cl  does  not  appear  to  be  a 
strong  function  of  the  forcing  amplitude. 


7.  Streak  iSow  visualization  and  streamline  patterns 

For  the  photographs  in  Figs.  7-9.  dye  was  introduced  upstream  of  the  cylin¬ 
der  near  the  mean  “stagnation”  streamline.  The  overlays  on  the  right  side  of  the 
figures  depict  the  streamline  patterns  calculated  from  a  potential  flow  model  (with 
the  cylinder  bounded  from  above  and  below  as  in  the  photographs)  using  the  lift 
coefficients  calculated  from  the  present  experiments.  As  can  be  seen,  the  stagnation 
streakline  in  the  photographs  can  be  approximated  by  the  stagnation  stream  line 
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Fig.  7  Top:  stationary  cylinder,  =  (0,0,0),  Cl  =■  0.0.  Bottom: 

oscillation  with  zero  mean  rotation  rate,  Sf)  —  (0,2.1, 0.7),  Cl  = 

0.0. 


Fig.  8  Top:  steady  rotation,  =  (0.5, 0,0),  Cl  =  0.6.  Bottom:  os¬ 
cillation  with  same  mean  rotation  rate  as  “Top:”,  =  (0.5, 2. 1,0. 7), 

Cl  =  2.0. 


in  the  potential  flow  case. 
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Fig.  9  Top:  steady  rotation,  (flo.lli,5/)  =  (4,0,0),  Ci  =  10.5.  Bottom:  oscil¬ 
lation  with  same  mean  rotation  rate  as  “Top:’’,  {Qo,Cli,Sf)  =  (4, 2. 1,0. 7), 
Cl  =  9.2. 

8.  Conclusions 

These  experiments  have  shown  that  higher  cylinder  aspect  ratios  yield  higher 
maximnm  lift  coefficients.  The  maximnm  lift  coefficients  in  the  present  experiments 
exceed  that  proposed  by  Prandtl  (1925),  possibly  because  Prandtl’s  arguments  ne¬ 
glected  unsteady  effects.  It  was  also  found  that  the  addition  of  forced  rotary  os¬ 
cillations  to  the  steady  rotation  of  the  cylinder  increased  Cl  in  cases  where  the 
wake  would  normally  be  separated  (in  the  steadily  rotating  case),  but  decreased  it 
otherwise.  Finally,  the  virtual  vortex  method  was  found  to  be  successful  in  charac¬ 
terizing  the  flow  ahead  of  the  cylinder.  The  method  yielded  estimates  for  Cl  which 
agree  with  the  data  published  by  Thom  (1924)  for  lift  coefficients  greater  than  2, 
t.e.,  flo  >  1-  At  the  smaller  values  of  Oq  <  1,  the  measurements  yielded  values 
of  Cl,  that  are  higher.  This  discrepancy  is  attribnted  to  a  difference  in  the  Re 
between  the  experiments. 

This  research  was  sponsored  by  the  Air  Force  Office  of  Scientific  Research, 
URI  AFOSR  Grant  No.  F49620-86-C-0134. 
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Streamline  patterns  in  the  neighborhood  of 
a  cylinder  executing  rotary  oscillations 
in  a  uniform  flow 


streamline  patterns  in  the  neighborhood  of 
a  cylinder  executing  rotary  oscillations 
in  a  uniform  flow 
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Graduate  Aeronautical  Laboratories,  California  Institute  of  Technology  301-40, 

Pasadena,  CA  91125  USA 

A  particular  case  of  a  circular  cylinder  executing  forced  rotary  oscillations  in  a 
uniform  stream  is  examined.  Estimates  of  the  unsteady  streamfunction,  averaged  at 
constant  phase,  are  presented.  These  data  show  that  vortical  structures  are  formed 
very  close  to  the  cylinder.  The  estimated  strength  of  these  vortical  structures  is 
small  when  compared  with  the  fluctuation  of  the  vorticity,  in  the  boundary  layer, 
caused  by  the  cylinder  motion.  A  po.ssible  (natural)  feedback  mechanism  limiting 
the  strength  of  the  vortices  in  the  wake  is  also  outlined. 


1.  Introduction 

The  flow  past  a  cylinder  executing  harmonic  rotary  oscillations  about  its  axis 
in  a  free  stream  has  been  approached  by  several  investigators,  using  flow  visualiza¬ 
tion,  force  measurements,  and  simulations. 

Okajima  et  al.  examined  this  flow  using  simulations  and  force  measurements 
for  Reynolds  numbers  in  the  range, 

40  <  /?e  =  <  6  X  10" 

and  normalized  amplitudes  (peak  rotation  raUi).  in  the  range 

0  a 

0.2  <  =  —  <  1.0  , 

U  oc 
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and  forcing  Strouhal  numbers  in  the  range 

2a  f 

0.05  <  Sf  =  <  0.3  , 

^QO 

where  a  is  the  cylinder  radius,  Uoo  tbe  free-stream  velocity,  i/  the  kinematic  vis¬ 
cosity,  /  the  forcing  frequency,  and  6  the  angular  velocity  of  the  cylinder.  They 
noted  a  “phase  locking”  of  the  vortex  shedding  (generation)  to  the  cylinder  motion 
over  a  range  of  forcing  frequencies  and  amplitudes.  Flows  about  a  cylinder  execut¬ 
ing  rotary  motions  have  also  been  investigated  by  Taneda  (1978),  Wu  et  al.  (1989), 
and  Filler  et  aZ.(1991).  The  present  data  is  an  outgrowth  of  the  work  presented  in 
Tokumaru  &  Dimotakis  (1991),  and  Tokumaru  (1991). 

In  the  work  to  be  discussed  here,  we  examine  a  particular  case  of  phase- 
locked  vortex  shedding  from  a  cylinder  executing  harmonic  rotary  oscillations  in 
a  free  stream.  The  effect  of  the  cylinder  oscillation  is  examined  using  velocity 
measurements  in  the  vicinity  of  the  cylinder,  averaged  at  constant  forcing  phases. 
While  not  perfect,  these  experimental  data  demonstrate  the  forced  shedding  process 
in  more  detail  than  flow  visualizations  and  single  point  measurements  that  can  be 
found  in  the  literature. 


2.  Experimental  facility 

The  present  experiments  were  performed  in  the  GALCIT  20"  x  20"  Free 
Surface  Water  Tunnel*  (see  Ward  197G).  A  machined  anodized  aluminum  cylinder, 
4"  in  diameter  with  an  aspect  ratio  of  about  4.2  was  supported  10"  above  the 
bottom  of  the  water  channel,  between  0.5"  thick  Plexiglas  fairings.  Power  was 
transferred  from  the  motor  to  the  cylinder  using  steel  reinforced  neoprene  timing 
belts.  These  were  entirely  enclosed  in  the  fairings  and  did  not  interfere  with  the 
flow.  The  cyhnder  was  driven  by  a  high  performance  DC  motor  from  PMI  Motion 
Technologies  ( JR16M4CH  Servodisk"^*^  ),  capable  of  tracking  an  arbitrary  command 
signal  within  the  bandwidth  and  slewing  rates  that  were  investigated.  The  command 

*  The  setup  of  the  present  expcriiiiciits  is  similar  to  that  used  in  Tokumaru  &  Dimotakis  (1991), 
and  Tokumaru  (1991). 
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signal  for  the  cylinder  angular  velocity  was  generated  by  a  computer-controlled 
function  generator.  The  angular  position  of  the  cylinder  was  measured  using  a 
BEI  Motion  Systems  Co.  13-bit  absolute  position  encoder  connected  to  the  motor 
driveshaft. 

Velocities  were  measured  using  a  Bragg  cell  frequency-shifted  laser  Doppler 
velocimeter  with  a  translatable  (and  rotatable)  focal  volume.  A  Lexel  Model  85 
Argon-ion  laser  operating  in  single  line  mode  ( 514.5  nm)  at  about  200  mW  was 
used  for  the  LDV.  The  laser  beam  for  the  LDV  velocity  measurements  was  split 
into  two  beams  of  nearly  the  same  intensity  using  a  cube  beam  splitter.  The  beams 
were  then  Bragg  cell  frequency-shifted  to  allow  for  flow  reversals.  The  Bragg  cell 
offset  between  the  two  beams  was  200  kHz.  A  Dove  prism  was  introduced  in  line 
after  the  Bragg  cells  for  the  purpose  of  rotating  the  orientation  of  the  measurement 
volume.  An  achromat  was  used  to  focus  the  beams  in  a  focal  (measurement)  volume 
located  in  the  midspan  plane  of  the  water  tunnel.  The  light  scattered  from  the  focal 
volume  was  gathered  by  a  pair  of  achromats  onto  a  photodiode  with  an  integral 
low  noise  pre-amplifier  designed  by  Dan  Lang.  The  signal  from  the  pre-amplifier 
was  then  band-pass  filtered  above  and  below  the  frequency  corresponding  to  the 
mean  flow  velocity.  The  band-pass  range  was  chosen  to  pass  instantaneous  velocity 
fluctuations.  A  tracking  phase-locked  loop  was  then  used  to  lock  a  TTL  square  wave 
to  the  dominant  frequency.  The  TTL  signal  was  read  by  a  counter-timer  board  on  a 
data  acqiiisition  computer  and  stored  on  disk  for  later  processing.  The  transmitting 
and  receiving  hardware  was  mounted  on  a  motor  driven  x-y-z  traverse  (positioner). 


3.  Results  and  discussion 

The  cylinder  was  rotated  sinusoidally  in  time,  i.e.,  with  a  rotation  rate,  9{t) , 
of 

0(t)  = -01  sin(27r/(),  (1) 

where  /  denotes  the  forcing  frequency,  a  the  cylinder  radius,  and  0i  the  peak 
rotation  rate  of  the  cylinder.  For  the  present  exi)eriments  the  forcing  Strouhal 
number  5/,  normalized  peak  rotation  .  and  the  Reynolds  number  Re,  were  all 
fixed  {Sf  =  0.9,  Di  =2,  and  Re  =  2.1  x  10'* ). 
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3.1  Velocity  measurements 

The  angular  velocity  component  in  the  flow,  vg(r,6),  was  measured  along 
eleven  9  =  constant  rays  in  the  upper  half  plane,  at  15  degree  intervals.  The  points 
in  Fig.  1  denote  the  measurement  volumes  used  in  the  present  experiment.  Each 
forcing  period  was  divided  up  into  fifty  (f6t  =  0.02)  time  slices  (phases).  The 
velocities  at  each  location  were  averaged  at  constant  phase  {i.e.,  at  each  of  the  fifty 
forcing  phases),  over  approximately  fourty-one  forcing  cycles. 

The  velocity  data  were  smoothed  in  time  using  a  Gaussian  filter  with  a  3  dB 
cutoff  around  five  times  the  forcing  frequency.  No  other  explicit  smoothing  was 
performed,  but,  because  of  the  wide  spacing  of  the  9  =  constant  rays,  there  was  an 
implicit  smoothing  of  data  interpolated  in  between. 


3.2  Streamfunction 


Calculation  of  a  stream  function  from  velocity  data  was  relatively  straightfor¬ 
ward.  It  was  calculated  at  each  forcing  phase  in  order  to  visualize  the  development 
of  the  flow  with  time.  In  particular,  given  that  the  stream  function  'F  is  known 
at  some  point  in  the  flow  xq  ,  and  the  normal  velocity,  =  u  •  n ,  to  any  path  C 
connecting  two  points  x  and  Xq,  is  also  known,  the  (dimensionless)  stream  function 
•0(x)  can  be  calculated,  for  two-dimensional  flows,  by  the  expression 


i/’{x/a) 


^(x) 


f  ^  —  +  (Xo/a)  . 

Jc  Uoo  a 


(2) 


Further,  since  the  stream  function  is  constant  along  a  solid  surface,  it  is  convenient 
to  begin  integrating  from  there,  and  to  choose  some  '0o(xo/a)  =  constant  on  that 
surface.  For  the  present  experiments,  V’o  was  chosen  to  be  zero  on  the  cylinder 
surface,  so  Eq.  2  could  be  written  «is 


V>(r/a) 


veir/a) 


dr  fa  , 


(3) 


where  r  is  the  radial  coordinate  with  r  =  0  on  the  cylinder  axis. 
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The  stream  function  data  were  first  estimated  along  each  B  =  constant  ray 
using  Eq.  3  and  a  simple  trapezoidal  integration  scheme.  Between  the  measurement 
points,  the  streamfunction  was  estimated  using  spline  interpolation,  first  along  the 
rays,  and  then  between  them.  The  velocities  in  the  lower  half-plane  were  approxi¬ 
mated  using  the  measurements  from  the  upper  half-plane,  but  a  half  forcing  cycle 
out  of  phase.  The  streamline  plots  were  generated  using  a  contour  plotting  routine 
with  points  interpolated  on  a  300  x  400  rectangular  grid. 

Figures  2  to  11,  show  selected  stream  lines  calculated  from  the  velocity  data 
averaged  at  constant  forcing  phase.  The  dot  in  the  figures  inside  the  cylinder 
mark  the  approximate  cylinder  position,  the  arrows,  the  direction  of  rotation,  and 
the  shaded  regions,  the  range  of  motion.  Notes  on  the  averaged  constant-phase 
streamline  data  are  contained  in  the  captions  of  the  Figs.  2  -  11.** 

In  Fig.  2  (/t  =  0.0),  the  cylinder  has  just  finished  shedding  a  vortex  (at 
xfa  =  1.5),  and  the  kernel  (bump)  of  the  next  vortex  to  be  shed  can  be  seen  at 
about  45  degrees  on  the  cylinder  surface.  The  bump  has  grown  larger  in  Fig.  3 
(/t  =  0.1),  note  also  the  beginnings  of  a  second  bump  on  the  cylinder  surface  at 
about  90  degrees.  In  Fig. 4  {ft  —  0.2),  both  the  main  bump  and  secondary  bump 
have  grown  larger,  and  moved  around  the  cylinder  in  the  clockwise  direction.  There 
is  also  a  large  vortical  structure  being  shed  off  the  lower  surface  at  the  same  time. 
In  Fig.  5  [ft  =  0.3),  the  bumps  continue  to  grow.  The  main  bump  seems  to  be 
held  in  place  relative  to  the  cyliiuler  smface.  possibly  by  the  flow  induced  by  the 
previously  shed  vortices.  The  cylinder  motion  is  in  the  clockwise  direction  which, 
intuitively,  would  tend  to  have  the  opposite  effect.  In  Fig.  6  {ft  =  0.4),  again,  the 
bumps  continue  to  grow,  but  they  now  appear  to  be  more  rounded.  The  larger 
bump  looks  as  if  it  is  being  “pinched”  into  the  flow.  In  Fig.  7  {ft  =  0.5),  the 
presence  of  the  secondary  bump  is  less  noticeable;  most  likely  it  has  been  merged 
with  the  primary  structure.  In  Fig.  8  {ft  =  0.6),  the  cylinder  begins  to  move  in 
the  opposite  direction  to  the  flow  on  that  side.  This  appears  to  be  when  the  vortex 
begins  to  be  released  into  the  flow.  Note  in  Fig.  9  {ft  =  0.7),  the  beginnings  of 
the  next  vortex  to  be  shed.  Then,  in  Fig.  10  (/t  =  0.8),  the  main  vortex  begins  to 

**  PostScript^  plots  and  averaged  velocity  data  for  the  entire  set  of  50  phases  of  the  flow  can  be 
obtained  in  electronic  form  by  c-inailiiig  a  rcfpicst  to  arcliive@hydraO.caltcch.edu. 
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move  off  to  the  right.  The  next  vortex  to  be  shed  continues  to  be  held  in  place  both 
by  the  vortices  which  have  been  shed,  and  the  rotation  of  the  cylinder.  Finally,  in 
Fig.  11  {ft  =  0.9),  the  main  vortex  continues  to  move  off  to  the  right.  Figure  12 
illustrates  the  stream  function  averaged  over  all  forcing  phases.  Note  the  small  size 
of  the  mean  recirculation  region. 

These  figures  show  that  the  forced  shedding  mechanism  outlined  in  these 
figures  replaces  the  natural  shedding  mechanism  with  one  that  gathers  vorticity 
close  to  the  cylinder  surface,  and  then  releases  it  as  a  vortical  structure  into  the 
flow.  This  is  in  contrast  to  observations  of  unforced  vortex  shedding  which  show 
that  vorticity  is  introduced  into  the  flow  from  more-or-less  steady  separation  points 
with  the  resulting  vortical  structures  forming  in  the  wake  of  the  cylinder. 


3,3  Strength  and  speed  of  shed  vortical  structures 

The  angular  velocity  data  in  the  present  experiment  is  also  well  suited  to 
calculating  the  circulation  Fc  around  the  cylinder, 

^  ,  r  /  (4) 

Uqc^  <Z  Jc{r)  b'oo 

where  r  is  the  radius  of  integration  around  a  circular  contour  C{r)  concentric  with 
the  cylinder. 


The  data  in  Fig.  13  is  typical  of  the  variation  of  Fc/C/oo^  with  radius  of 
integration  r/a  at  a  constant  phase;.  Figure  14  is  the  streamline  pattern  related  to 
the  data  in  Fig.  13.  Figure  15  shows  the  vaxiation  of  Fc/t/oo^  with  ft  and  r/a. 


The  circulation  inside  an  annulus  between  two  circles  of  radii  ri  and  r2  can 
be  written  as 

r..  ir2  TJr.\  rJrA 

(5) 


Fc(r2)  FcCn) 


f  1 


Hence,  an  estimate  for  the  strength  of  a  vortical  structure  (residing  in  an  annulus) 
is 

F--  F(; ( 7'iiiax )  1  Fc{r,„i„)| 


U  on  O. 


Uora 


FcfT-niin) 


UnnO, 
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where  r^ax  aiid  r,iiiu  are  chosen  to  maximize  Eq.  6  in  the  wake.  Using  Eq.  5  and 
data  like  that  shown  in  Fig.  13  and  15  to  estimate  the  strength  of  the  vortical 
structures  yields 


(7) 


for  the  present  exi)eriment. 


Note,  however,  that  the  (peak-to-peak)  fluctuations  of  the  circulation  at  the 
cylinder  surface 


r  _ 2  r'c,iuax(^  —  O.) 

U^d  Ufx^Q, 


47rfli  ~  25  , 


for  the  present  exi>erinient,  is  much  larger  than  in  Eq.  7. 


(8) 


This  may  be  because  the  strength  of  the  vortices  in  a  wake  pattern  is  nec¬ 
essarily  limited  since  the  pattern  cannot  have  a  self-induced  velocity  greater  than 
the  free  stream,  i.e.,  the  wake  pattern  cannot  convect  upstream  faster  than  the 
cylinder.  This  suggest  that  a  natural  feedback  mechanism  is  at  work  limiting  the 
strength  of  the  vortices  well  below  values  like  Eq.  8.  For  example,  if  the  vortices  are 
too  strong  they  will  tend  to  move  downstream  more  slowly,  delay  the  shedding  of 
the  next  vortical  structure  into  the  wake.  The  oscillations  of  the  cylinder  generate 
both  positive  and  negative  vorticity  over  the  entire  surface  of  the  cylinder  during 
each  forcing  phase.  This  allows  the  vortex  to  gather  vorticity  of  both  signs  during 
a  forcing  cycle.  The  strength  of  the  vortices  then  depends  on  the  shedding  phase, 
where  the  shedding  phase  is  regtdated  by  previously  shed  vortices.  In  some  sense, 
this  feedback  can  be  likened  to  a  phase-lock  loop  circuit. 


Finally,  from  the  slope  of  the  contours  in  Fig.  15,  it  is  found  that  the  average 
(normalized)  convection  speed,  Uy/Uoo,  of  the  vortices  relative  to  the  cylinder, 
based  on  the  slope  of  the  Fc/f/c^n  =  0  contour,  is 


?/v  Sf  d(7-fa) 

uZ  "  T  d{ft)  ~ 


r 

-  >  1.2  . 
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for 
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3.4  A  rough  estimate  of  the  vortex  strength 


Using  the  Karman  vortex  street  as  a  model,  an  order  of  magnitude  estimate  of 
the  strength  per  unit  streamwisc  sjjacing  of  the  vortical  structures  can  be  derived. 
From  Goldstein  (1938,  §244), 


Fy 

Uoc  O. 


(10) 


where  £x  the  streamwise  spacing  of  the  vortical  structures.  In  addition,  £x/a  can 
be  related  to  Uy/f/oo  and  S/  by 


Lk 


(11) 


Note  that  by  construction  Fy  >0  ( Uy  <  1 )  for  a  wake,  and  a  value  of  Uy  <  0 
would  indicate  a  wake  pattern  that  is  traveling  upstream  faster  than  the  cylinder, 
so  for  steady  vortex  shedding  it  is  expected  that 


0  < 


tty 

uZ 


<  1 . 


Taking  Eqns.  10  and  11  and  applying  Eq.  12  then  yields 


(12) 
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Sf 


1.5  , 


(13) 


for  the  present  experiment.  This  is  consistent  with  the  value  of  Ty/UooO.  ~  1  in 
Eq.  7. 
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3.5  Caveats 

It  should  be  noted  that  there  are  several  limitations  to  the  present  results. 
For  a  few  measurements  near  the  cylinder  surface,  the  boundary  layer  could  not 
be  resolved  with  the  present  setup.  This  was  because  of  large  graxiients  in  the 
velocity  and  difficulties  in  positioning  the  measurement  volume.  This  did  not  cause 
large  errors  in  stream  function  integration,  however,  because  the  velocity  of  the 
cylinder  surface  was  known  a  piioii.  For  the  integration  at  the  cylinder  surface,  a 
(straight)  line  was  fit  between  the  velocity  at  the  wall,  and  the  nearest  measured 
velocity.  As  can  be  seen  from  Fig.  1,  the  spatial  resolution  of  the  data  for  rja  >  1.5 
corresponds  to  only  a  handful  of  measurement  locations  per  vortical  structure  in 
that  region.  Continuity  of  the  velocity  field,  however,  diminishes  this  problem.  In 
addition  the  data  in  the  lower  half  plane  was  taken  to  be  the  same  as  in  the  upper 
half  plane,  but  a  half  forcing  cycle  out  of  phase.  Because  of  asymmetries  in  the  flow, 
this  match  (at  y  =  0)  is  not  perfect  and  causes  small  “glitches”  in  the  streamline 
pattern  around  the  joint  between  the  upper  and  lower  half  planes.  Comparing 
streamline  patterns  a  half  period  apart  gives  an  indication  of  the  asymmetry  in  the 
data.  Because  the  present  data  was  measured  at  midspan,  there  is  no  out-of-plane 
velocity,  in  the  mean.  This  does  not,  however,  remove  three-dimensional  effects  like 
vortex  stretching  and  instantaneous  out-of-plane  flow.  Finally,  because  the  data 
was  averaged  at  constant  phase,  fluctuations  in  the  flow  that  are  not  in  phase  with 
the  forcing  will  be  averaged  out.  Because  of  factors  such  as  these,  comparisons  of 
the  present  data  with  two-dimensional  simulations  will  not  be  exact.  It  should  be 
emphasized,  however,  that  the  basic  mechanisms  gathering  vorticity  close  to  the 
cylinder  (allowing  it  to  mix  and  cancel)  and  then  ejecting  it  into  the  flow,  should 
not  be  strongly  affected  by  these  factors.  A  good  simulation  (that  captures  the 
behavior  of  the  unsteady  boundary  layer)  should  yield  results  similar  to  the  present 
data. 
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4.  Conclusions 

Estimates  of  the  unsteady  streamfunction,  averaged  at  constant  phase,  were 
made.  They  revealed  that  the  vortical  structures  are  formed  very  close  to  the  cylin¬ 
der  and  then  released  into  the  wake.  This  is  in  contrast  to  the  natural  (unforced) 
shedding  of  vortices,  where  the  structures  form  after  the  vorticity  has  been  shed 
into  the  wake,  from  more-or-less  steady  separation  points. 

In  addition  is  was  found  that  the  strength  of  these  structures  was  small, 
by  more  than  an  order  of  magnitude,  when  compared  with  the  fluctuation  of  the 
vorticity  in  the  boundary  layer  caused  by  the  cylinder  rotational  acceleration.  It 
is  proposed  that  this  is  because  the  strength  of  the  vortices  in  a  wake  pattern  is 
necessarily  limited  since  the  pattern  cannot  have  a  self-induced  velocity  greater 
than  the  free  stream  (i.c.,  the  wake  pattern  cannot  convect  upstream  faster  than 
the  cylinder).  Evidently,  most  of  this  vorticity,  introduced  into  the  boundary  layer 
fluid  combines  before  the  fluid  is  allowed  to  enter  the  wake. 

This  research  was  sponsored  by  the  Air  Force  Office  of  Scientific  Research, 
URI  AFOSR  Grant  No.  F49620-8G-C-0134. 
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Fig.  1  Measurement  locations. 
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Fig.  11  ft  =  0.9. 
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Fig.  12  Streamlines  averaged  over  all  phases.  Note  the  delayed  average  separation 
point  and  the  small  averaged  “recirculation”  region  realized  in  this  forcing 
regime.  This  looks  very  much  like  the  (unforced)  time  averaged  flow  around 
t’le  cylinder  at  higher  Reynolds  numbers  (>  5  x  10^),  which  also  has  a 
similar  drag  coefficient  of  about  0.2  to  0.3. 
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Fig.  13  Variation  of  the  normalized  circulation,  rc/f/cx:«i  t^'-ken  round  contours  of 
varying  radius,  rfa  {ft  =  0.54).  The  rise  in  the  curve  between  r/a  =  1.2 
and  1.7  marks  a  patch  of  positive  (clockwise)  circulation  and  the  drop 
between  r/a  =  1.7  and  2.3  a  negative  (counter-clockwise)  patch,  cf 
Fig.  14. 


16  April  1992 


Fig.  15  Fc/UooO.  vs.  r/a  and  ft.  The  dashed  lines  denote  net  counter-clockwise 
rotation;  the  solid  lines,  clockwise.  The  value  Tc/UocO-  at  contour  level 
T’  =  0.6,  ‘2’  =  0.5,  ‘3’  =0.4  ...,  down  to,  T3’  =  —0.6.  Contour  level 
‘7’,  solid  line,  marks  the  zero  contour. 
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